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In  the  three  year  period  of  this  contract  we  have  built 
and  installed  a  long  baseline  stellar  interferometer  on  Mt. 
Wilson.  A  number  of  papers  describing  the  instrument  and 
results  from  that  instrument  have  appeared  [Shao  et  al. 
Astronomy  and  Astrophysics  193 (1988) : 357 ;  Shao  et  al. 
Astrophysical  Journal  327 (1988) :905;  Colavita  et  al.  Applied 
Optics  26(1987) :4106;  Colavita  et  al.  Applied  Optics 
26(1987) ;4113;  Mozurkewich  et  al.  Astronomical  Journal 
95 (1988) : 1269] .  The  Mark  III  interferometer  is  the  most 
astronomically  productive  interferometer  to  date.  The 
success  of  the  Mark  III  has  led  the  USNO  to  decide  to  build  a 
dedicated  four  element  interferometer  array  for  fundamental 
astrometry  and  NRL/ONR  to  build  a  six  element  imaging  array. ^ 

To  summarize  the  last  three  years,  we  list  the  major  / 
milestones. 


9/86  first  fringes  with  12M  N-S  baseline 


11/86  automated  operation  (165  obs/night) 
11/86  first  precise  diameter  measurements 


5/87  8.6M  E-S  baseline  in  operation 

7/87  star  tracker  using  photon  camera  and  correlation 
centroid  estimator  put  into  operation 

8/87  modifications  to  enable  simultaneous  astrometric 
and  diajp^ter  measurements 

9/87  install  laser  metrology  system  | 

5/88  E-S  baseline 


5/88  E-S  baseline 

5/88  /  design  and  build  variable  baseline  (3-30  meters) 
/  components 


6/?8  design  and  build  group  delay  correlator 


ziKc 


In  conclusion,  we  have  successfully  completed  the  task 
of  building  a  long  baseline  interferometric  observatory. 
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Summary.  The  Mark  III  interferometer  is  an  operational  long 
baseline  stellar  interferometer  on  Mt.  Wilson  with  four  pt>ssible 
baseline  configurations  from  9m  NE-SW  to  20m  N-S.  The 
interferometer  was  designed  to  be  a  highly  automated  astronom¬ 
ical  instrument  to  measure  stellar  positions  and  diameters  to  a 
magnitude  limit  of  seven.  Initial  fringe  observations  were  made 
in  September  1986  with  a  12-m  N-S  baseline.  In  the  following 
months,  semi-automated  astrometric  and  stellar  diameter  mea¬ 
surements  were  also  made.  This  paper  describes  the  hardware  and 
software  components  of  the  instrument  and  its  operational 
characteristics. 

The  interferometer  has  several  novel  features.  One  is  the  use  of 
optimal  estimation  and  control  algorithms  (e.  g.  Kalman  filters)  in 
the  control  loops.  Another  is  the  ability  to  operate  both  as  a 
closed-loop  phased  interferometer  and  eventually  as  an  open-Untp 
or  absolute  coherent  interferometer.  High  thermal  stability  and 
mechanical  accuracy  should  permit  the  instrument  to  psiint  blind 
at  an  astronomical  object  and  maintain  optical  path  equality  to 
within  the  limits  set  by  the  atmosphere.  In  this  absolute  inter¬ 
ferometric  mode  of  operation,  it  should  be  possible  to  observe 
faint  astronomical  objects  that  are  too  dim  for  phase  tracking.  In 
theory,  measurements  of  amplitude,  group  delay,  and  closure 
phase  will  be  possible  to  14  mag. 

Key  words:  stellar  interferometry  -  astrometry  -  stellar  diameters 
optical  array  -  proper  motion 


I.  Introduction 

The  Mark  HI  interferometer  is  the  third  of  a  scries  of  stellar 
interferometers  at  Mt.  Wilson  Observatory.  The  first  instrument 
demonstrated  white-light  fringe  tracking  in  1979  (Shao  and 
Staelin.  1980).  The  second  instrument  observed  fringes  in  1982 
and  was  used  as  a  technology  test  bed  for  astrometric  measure¬ 
ments  until  1984  (Shao  et  al..  1987).  Several  devices  critical  to  the 
operation  of  a  stellar  interferometer  were  developed  for  the  Mark 
II  interferometer.  One  such  device  was  a  high  speed,  ultra  high 
accuracy,  laser  controlled  optical  delay  line.  An  operational 
procedure  for  rapid  switching  between  stars  was  also  developed. 


Traditionally,  the  most  difficult  task  of  a  stellar  interferometer  has 
been  finding  and  tracking  the  central  fringe.  In  rapid  switching, 
this  procedure  is  totally  automated.  Thus,  operation  of  the  Mark 
II  helped  identify  the  technical  problems  of  making  astrometric 
measurements  with  a  long  baseline  interferometer. 

Like  the  previous  instrument,  the  Mark  III  interferometer  was 
also  designed  to  lie  a  test  bed  for  new  techniques  and  technologies 
relevant  to  optical  interferometry  and  synthetic  optical  aperture 
synthesis.  However,  the  principal  purpose  of  the  Mark  III 
interferometer  was  to  initiate  a  test  program  of  fundamental 
astrometric  measurements  to  demonstrate  a  significant  improve¬ 
ment  in  accuracy.  Another  goal  was  to  initiate  a  program  of 
accurate  stellar  diameter  measurements.  Most  important  of  all. 
the  instrument  was  designed  to  be  astronomically  productive. 

Long  baseline  stellar  interferometers  have  traditionally  been 
considered  one  of  the  most  technologically  challenging  astronom¬ 
ical  instruments.  In  the  Mark  III  interferometer,  the  goal  has  been 
to  build  an  instrument  that  is  reliable,  easy  to  operate,  and  capable 
of  extremely  accurate  astronomical  measurements.  In  order  to 
achieve  these  goals,  a  number  of  active  subsystems  were  in¬ 
corporated  into  the  design  of  the  instrument.  In  implementing 
these  subsystems,  the  design  philosophy  involved  keeping  the 
mechanical  and  optical  systems  as  simple  as  possible  and  shifting 
the  bulk  of  the  complexity  into  software. 

Although  the  Mark  III  was  originally  designed  as  a  closed  loop 
phased  interferometer,  operational  experience  has  shown  that  it 
has  sufficient  accuracy  and  stability  to  operate  as  an  absolute 
interferometer.  An  absolute  interferometer  is  one  whose  mechan¬ 
ical  accuracy  and  stability  in  open  loop  operation,  without  fringe 
tracking,  is  sufficiently  high  that  optical  path  errors  are  dominated 
by  atmospheric  turbulence.  This  level  of  accuracy  and  stability, 
usually  associated  with  the  primary  mirror  of  a  large  astronomical 
telescope,  is  the  level  of  accuracy  exhibited  by  the  Mark  III 
interferometer  with  the  1 2-m  baseline.  Operation  with  the  20-m 
baseline  is  expected  to  exhibit  the  same  level  of  accuracy  and 
stability. 

The  ability  to  operate  as  an  absolute  interferometer  makes 
possible  three  new  modes  of  interferometric  observations  for  long 
baseline  interferometers.  They  are  photon  starved  amplitude 
measurements,  photon  starved  group  delay  astrometry,  and 
photon  starved  closure  phase  measurements.  Photon  starved 
amplitude  measurements  have  long  been  a  part  of  speckle 
interferometry  with  single  apertures.  The  importance  of  photon 
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Fig.  1.  The  Mark  III  interrcrometer  on  Ml.  Wilson 


Starved  observations  is  that  they  make  possible  astrometry  and 
aperture  synthesis  of  faint  extragalactic  objects.  These  three  types 
of  measurements  are  discussed  after  the  description  of  the 
instrument. 


2.  The  instrument 

The  interferometer  is  located  in  California  on  Mt.Wilstm  ap¬ 
proximately  80  m  east  of  the  1 00-inch  telescope,  near  the  site  of  the 
Mark  I[  interferometer.  Figure  I  shows  a  view  of  the  site.  Six  and 
ten  meters  north,  south,  and  east  of  the  central  building  are  the 
siderostat  piers,  which  are  e.ich  constructed  of  steel-reinforced 
concrete  with  a  4  x  4  foot  base  and  a  height  of  8  feet,  4  feet  of  which 
is  above  ground.  Siderostats  are  currently  mounted  on  the 
innermost  piers.  Huts  provide  weather  protection  for  the  sidero¬ 
stats.  These  have  hinged  roofs  which  swing  back  exposing  the 
siderostats  to  the  sky.  Currently,  there  are  vacuum  pipes  between 
the  huts  and  the  central  building  to  eliminate  local  seeing  effects 
caused  by  heat  in  the  main  building.  Initial  operation  involved 
only  the  t2-m  baseline  connecting  the  north  and  .south  siderostats: 
however,  the  east  siderostat  can  be  used  in  place  of  the  north 
siderostat  for  a  9-m  NE-SW  baseline.  The  fourth  window  on  the 
central  building  is  for  future  expansion.  In  1988.  we  plan  to 
operate  on  the  outer  set  of  piers  shown  in  Fig.  I.  with  20-m  N-S 
and  16-m  NE-SW  baselines. 

Since  the  instrument's  principal  purpose  is  large  angle  astro¬ 
metry,  attention  was  paid  to  construction  practices  that  affect 
astrometric  accuracy.  The  central  building  is  air  conditioned 
during  the  day  so  that  at  night  there  will  be  a  minimum 
temperature  difference  between  inside  and  outside.  The  piers  in 
the  huts  that  hold  the  siderostats  go  down  to  bedrock  and  are  not 
directly  connected  to  the  concrete  pads  on  which  the  huts  are  built. 
The  piers  are  also  insulated  by  two  inches  of  foam.  The  optics 
inside  the  main  building  that  affect  astron.etric  measurements  are 
ail  located  on  a  super-invar  breadboard,  except  the  mirrors  in  the 
optical  delay  tine.  However,  the  position  of  the  optical  delay  line  is 
monitored  by  a  laser  interferometer  referenced  to  the  same  super- 
invar  breadtoard.  The  laser  is  frequency  stabili7ed  with  better 
than  10"’  long-term  stability. 

The  instrument  can  be  divided  into  five  subsystems:  I )  the  star 
tracker  or  autoguider,  2)  the  optical  delay  line.  3)  the  stellar  fringe 


acquisition  and  tracking  system,  4)  the  laser  metrology  system, 
and  S)  the  siderostat  pointing  and  observation  sequencing  system. 
A  description  of  the  optical  layout  of  the  instrument  is  followed  by 
a  more  detailed  explanation  of  the  five  subsystems. 

2.1.  Optical  layout 

Figure  2  is  an  optical  schematic  of  the  interferometer.  Light  from 
the  two  siderostats  is  directed  toward  the  central  building  by  fixed 
mirrors.  These  fixed  mirrors  are  necessary  to  keep  the  polarization 
vectors  from  the  two  interfering  beams  parallel  for  observation  of 
stars  located  over  all  parts  of  the  sky.  After  the  light  enters  the 
main  building,  the  beams  are  directed  by  piezoelectrically  con¬ 
trolled  mirrors  towards  the  vacuum  delay  lines.  After  the  light 
reflects  from  the  delay  line  retroreflectors.  the  beams  are  com¬ 
bined  at  the  beamsplitter  B.  After  the  beamsplitter,  part  of  the 
light  is  directed  into  optical  fibers  which  feed  GaAs  phototubes 
located  some  distance  from  the  central  breadboard,  and  part  of 
the  light  is  directed  toward  the  star  tracker.  Other  optical 
components  include  the  laser  interferometers,  and  alignment 
lasers  and  telescopes. 

2.2.  Star  tracker 

The  star  tracker  has  been  described  in  detail  previously  (Clark  et 
al..  1986:  Clark.  1986)  and  we  will  provide  only  a  brief  summary 
here.  The  star  tracker  system  is  designed  to  keep  the  two 
interfering  wavefronts  parallel  to  within  a  small  fraction  of  an 
arcsecond.  As  with  all  servo  systems,  the  star  tracker  is  composed 
of  an  error  sensor,  a  controller,  and  an  actuator.  For  the  Mark  III 
interferometer,  two  different  error  sensors  were  built  becau.se  a 
key  component  (a  photon  counting  camera)  for  the  original  design 
was  not  delivered  before  initial  observations  were  made  in 
September  1986.  Only  the  star  tracker  design  that  will  be 
permanently  attached  to  the  interferometer  is  described. 

An  optical  schematic  of  the  star  tracker  is  shown  in  Fig.  3. 
Light  from  the  star  in  a  3-inch  beam,  intercepted  by  the  two 
siderostats.  is  combined  at  a  beamsplitter.  Light  from  the  south 
siderostat  passes  through  an  annular  wedge  that  deflects  the  outer 
portion  of  the  beam  by  45":  the  inner  2-inch  part  of  the  beam  is 
undetlecied.  After  the  light  is  combined  at  the  beamsplitter,  an 
annular  mirror  at  one  output  of  the  beamsplitter  directs  the  light 
in  the  annulus  towards  a  telescope  and  imaging  detector.  The  light 
in  the  central  2-inch  core  from  both  sides  of  the  beamsplitter  is  fed 
to  photomultipliers  that  detect  the  fringes. 

When  the  two  components  of  the  light  beam  that  enter  the 
fringe  detection  phototubes  are  parallel  (as  they  must  be  for 
proper  operation  of  the  interferometer),  the  two  components  of 
the  annular  beam  entering  the  telescope  will  he  tilted  by  45' 
relative  to  each  other  because  of  the  45"  wedge  in  the  south  arm  of 
the  interferometer.  Thus  two  star  images  will  appear  on  the  face  of 
the  photon  counting  imaging  detector.  The  optical  arrangement  of 
this  star  tracker  is  slightly  different  from  that  used  for  some  other 
interferometers.  The  present  design  using  one  telescope  for  both 
star  trackers  eliminates  the  need  to  align  the  two  optical  systems  to 
each  other  with  subarcsecond  accuracy.  More  important,  the 
relative  alignment  of  two  separate  optical  systems  would  drift  with 
lime  and  temperature,  whereas  the  star  trackers  for  the  Mark  III 
interferometer  will  maintain  alignment  indeHnitely. 

The  photon  camera  used  in  the  star  tracker  is  a  commercial 
detector  based  on  the  PAPA  camera  (Papaliolioset  al..  1985)  with 
256  X  256  pixels.  The  detector  produces  a  pulse  with  an  .v,  r 
address  for  every  detected  photon.  The  digital  output  of  the 
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Fig.  2.  Optical  schematic  of  the  Mark  III  interferometer 


camera  is  interfaced  to  a  high  speed  signal  processor  which  iwates 
the  centroids  of  the  two  stellar  images.  Two  additional  optical, 
mechanical  components  are  shown  in  Fig.  ,V  The  variable  neutral 
density  filter  is  used  to  reduce  the  photon  flux  to  the  camera  to 
prevent  overloading  of  the  camera  electronics  when  observing 
bright  stars.  The  setting  of  the  filter  is  controlled  by  the  star  tracker 
computer.  There  are  also  two  shutters,  one  in  each  arm  of  the 
interferometer,  which  are  used  for  star  acquisition.  When  the 
siderostats  stop  slewing,  the  two  stellar  images  may  be  sutTicicntly 
far  from  their  correct  positions  that  proper  identification,  i.e.. 
which  star  came  from  which  siderostat.  may  not  be  possible.  The 
shutters  are  used  to  remove  this  ambiguity. 

The  camera  is  the  error  sensor.  The  actuator  is  a  PZT-driven 
gymbaled  mirror.  This  5-inch  mirror  is  driven  in  two  orthogonal 
directions  with  a  30"  peak-to-peak  motion.  The  bearings  are  cross 
flexure  pivots  and  the  mirror  mount  has  its  lowest  resonance  at 
200  Hz. 

The  heart  of  the  star  tracker  is  the  control  system;  it  makes 
extensive  use  of  optimal  estimation  techniques  and  only  a  brief 
outline  is  given  here.  A  detailed  description  is  given  by  Clark  et  al. 
(1986).  The  control  problem  is  separated  into  a  spatial  filtering 
part  and  a  temporal  filtering  part.  The  spatial  filter  takes  the 
photon  addresses  generated  by  the  photon  camera  and  derives  the 


centroid  of  the  star.  The  temporal  filter  takes  current  and  past 
centroid  estimates  and  generates  a  set  of  voltages  to  drive  the 
piezoelectric  actuators. 

The  image  of  the  star  on  the  camera  face  is  the  diffraction 
pattern  of  an  annular  telescope  aperture.  A  cross-correlation 
centroid  estimator  was  implemented.  The  estimated  v  coordinate 
of  the  centroid  is  given  by 

.\  =  JJ  /(V.  r)  /i, (V.  i)(/.vf/v,  (I) 

where  /  is  the  detected  image  and  /i,  is  the  correlation  function  lor 
the  v  direction.  The  noise  in  the  detected  image  is  photon  noise 
with  a  variance  given  by 

rTj=/(.v,.i).  (2) 


Since  the  theoretical  .shape  of  the  image  is  known,  there  is  a  closed 
form  solution  for  the  optimal  correlation  function,  i.e..  that 
function  /it.v.  i  )  which  gives  a  minimum  variance  for  .v  due  to 
photon  noise  in  the  measurement  of  /(.v.  r).  The  optimal 
correlation  function  is  given  by 


l-,7(x.^)/r.v] 

/(X.A) 


(3) 
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Fig.  X  Optical  schematic  of  star  tracker 


Table  I 


Correlation  function 

Single  photon  variance 
(arcsec)^ 

Optimal  [Eq.  (3)] 

0.303 

Quad  cell  [Eq.  (4)] 

3.09 

Linear  [Eq.  (S)] 

6.14 

This  type  of  centroid  estimator  has  limited  dynamic  range  but  very 
low  variance.  An  analysis  of  this  estimator  for  our  annular 
aperture  image  was  made  relative  to  centroid  estimators  using 
quad  cell  and  linear  correlation  functions; 

Quad  cell:  ^  j  ^  ^  ^ .  (4) 

linear:  h  =  x.  (5) 

The  results  are  shown  in  Table  I ,  Note  that  the  variance  for  the 
optimal  correlator  estimator  is  about  an  order  of  magnitude  better 
than  a  quad  cell. 

The  correlation  integral  (Eq.  (I)]  is  performed  by  a  custom 
digital  signal  processor.  The  processor  calculates  the  integral  on  a 
photon  by  photon  basis,  rather  than  a  pixel  by  pixel  basis.  This  is 
more  efficient  computationally  when  the  number  of  detected 
photons  is  much  leU  than  the  number  of  pixels  in  the  image. 
Currently,  centroid  estimates  are  calculated  every  1 0  ms.  The 
signal  processor  can  operate  at  photon  rates  well  above  the 
capabilities  of  the  camera.  The  optical  correlation  functions,  one 
in  X,  and  one  in  y  for  each  of  the  two  trackers  are  downloaded  to 
the  correlation  processor  from  the  main  control  computer.  The 


calculation  of  the  optimal  correlation  function  is  not  done  in  real 
time. 

The  second  part  of  the  control  system  is  the  temporal  filter. 
While  the  spatial  filter  is  implemented  mostly  in  digital  hardware 
and  microcode,  the  temporal  filter  is  implemented  entirely  in 
software  on  the  star  tracker  control  computer.  The  temporal  filter 
is  a  discrete  time  Kalman  filter.  The  Kalman  filter  assumes  a  state 
variable  formulation  of  the  tracking  problem  and  is  the  linear 
filter  which  minimires  the  mean  square  tracking  error.  This  filter 
takes  into  account  the  observation  noise  (photon  noise)  by 
adaptively  varying  the  tracker  bandwidth  according  to  the  noise 
variance.  The  output  of  the  filter  is  scaled  and  then  directed  to  a 
DIA  converter  and  high  voltage  amplifier  to  drive  the  PZTs. 

In  addition  lo  the  software  for  the  Kalman  filter,  there  is 
software  for  acquisition,  data  recording,  and  for  utilities  such  as 
moving  the  shutters  and  controlling  the  camera.  There  is  also 
software  for  complete  automatic  operation  and  for  off-line 
diagnostics  and  analysis  of  tracker  performance. 

2.3.  Optical  delay  line 

Two  optical  delay  lines  are  used  to  equalize  the  paths  from  the  star 
to  the  beamsplitter  via  the  two  siderostats  12m  apart.  The 
mechanical  and  optical  arrangement  is  shown  in  Fig.  4.  The 
mechanical  design  is  similar  to  that  used  in  delay  lines  for  high 
resolution  Fourier  spectroscopy  (Connes.  1975).  The  optical  delay 
is  changed  by  moving  cat's  eye  retroreflectors  along  tracks  which 
are  housed  in  two  2S-fi  vacuum  chambers.  The  peak-to-peak 
optical  path  difference  of  this  system  is  20  m,  and  the  delay  lines 
can  be  positioned  under  laser  control  to  an  accuracy  of  10  nm,  for 
a  dynamic  range  of  2  10’  resolution  elements  in  optical  path 
difference.  The  maximum  slew  rate  is  I  m  per  second  of  optical 
path  difference.  Twenty  seconds  is  required  to  go  from  one  end  to 
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Fig.  5.  Hetcri>dync  laser  inlerl'cronK’ter  Tor  delay  line 


the  other.  The  maximum  tracking  rate  is  2  mm  s  ‘ .  at  which  point 
the  rms  position  accuracy  has  degraded  to  20 nm. 

Figure  4  shows  the  optical  path  of  the  laser  beam  that 
measures  the  position  of  the  cat's  eye  retroreflector.  The  cat's  eye 
is  an  8-inch  f/3  parabola  with  a  1  -inch  flat  mirror  at  its  focus.  The 
1-inch  mirror  is  mounted  on  a  PZT  stack  that  has  a  5-pm  pcak-to- 
pcak  range.  The  cat’s  eye  is  suspended  on  a  flexure  bearing  over  a 
belt  driven  translation  stage.  A  voice  coil  positions  the  cat's  eye 
over  the  stage,  while  a  motor  at  one  end  of  the  25-ft  track  moves 
the  stage  via  a  timing  belt.  The  voice  coil  has  a  2-mm  range  and  the 
motor  has  a  9-pm  step  size.  The  three  actuators  PZT,  voice  coil, 
and  motor  -  have  considerable  overlap  in  terms  of  resolution  and 
dynamic  range. 

The  delay  line  sensor  is  a  heterodyne  laser  interferometer 
(Fig.  S).  The  light  from  a  stabilized  laser  is  split  into  two  beams  by 
a  polarizing  beamsplitter.  They  are  then  frequency  shifted  by  two 
acousto-optic  bragg  cells  with  a  2-MHz  frequency  difference. 
When  recombined,  the  beam  has  its  two  orthogonal  polarizations 
at  different  frequencies  with  a  2-MHz  frequency  difference.  The 
laser  is  a  commercial  transverse  Zeeman  stabilized  HeNe  laser 
with  a  one-day  stability  of  200  kHz  at  line  center.  474 THz.  The 


interferometer  is  a  standard  polarization  heterodyne  interfero¬ 
meter.  The  pha.se  of  the  2-MHz  signal  from  the  interferometer  is 
compared  to  the  2-MHz  reference  (Fig.  5)  and  a  one  wavelength 
change  in  the  monitored  optical  path  causes  a  2n  change  in  the 
relative  phases  of  the  interferometer  and  reference  2-MHz  signals. 
A  custom  digital  phase  comparator  circuit  counts  the  fringes  with 
A/M  resolution  at  a  maximum  slew  rate  of  2  10’’  fringes  fier 
second. 

Since  the  desired  optical  path  difl'erence  is  given  by  the 
difference  between  the  positions  of  the  two  delay  lines,  it  is  only 
necessary  to  actively  control  one  delay  line.  Thus,  vibrations  in 
one  delay  line  are  compensated  by  driving  the  other  in  exactly  the 
same  way. 

The  laser  interferometers  are  read  at  a  1-kHz  rate  by  the  fringe 
tracker  computer.  The  difference  between  the  two  interferometers 
is  the  measured  position,  which  is  subtracted  from  the  desired 
position  to  give  an  error  signal.  The  PZT  servo  is  a  simple  first- 
order  system,  i.e.,  the  PZT  drive  voltage  is  just  the  integral  of  the 
error  signal.  The  necessary  integral  is  performed  in  digital 
hardware,  and  the  velocity  of  the  PZT  mirror  is  updated  every 
millisecond  by  the  computer.  The  closed-loop  bandwidth  of  this 
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servo  is  approximately  100  Hz.  for  a  time  constant  of  1  to  2  ms. 
The  small-signal  bandwidth  of  the  entire  delay  line  is  the 
bandwidth  of  this  innermost  servo. 

The  voice  coil  servo  is  designed  to  keep  the  PZT  system  within 
its  S-pm  dynamic  range.  The  cat's  eye  flexure  suspension  has  a  1 .5- 
Hz  resonance.  Since  a  S-pm  accuracy  requires  a  response  signifi¬ 
cantly  faster  than  1 .5  Hz.  the  voice  coil  servo  has  to  be  at  least  a 
second-order  servo.  The  voice  coil  servo  is  implemented  entirely  in 
software.  The  error  signal  input  to  the  voice  coil  servo  is  the 
distance  the  PZT  mirror  is  from  the  center  of  its  range.  The  output 
is  a  current  to  the  voice  coil.  The  closed-loop  bandwidth  of  this 
servo  is  approximately  10  Hz. 

The  motor  servo  is  designed  to  keep  the  voice  coil  in  the  center 
of  its  range.  The  input  to  this  servo  is  the  position  of  the  voice  coil 
as  measured  with  a  lateral-cell  transducer.  The  motor  is  a 
microstep-driven  stepper  motor,  interfaced  to  the  computer 
through  a  programmable  pulse  generator.  In  the  tracking  mode, 
it  operates  as  a  first-order  servo,  with  a  small-signal  closed-loop 
bandwidth  of  about  1  Hz. 

The  delay  line  control  system  was  designed  to  interact  with  the 
rest  of  the  interferometer  through  a  simple  interface.  The  rest  of 
the  interferometer  does  not  see  three  separate  actuators  and 
control  loops  and  a  laser  interferometer.  Rather,  the  delay  line 
subsystem  is  programmed  by  setting  its  position  and  velocity.  At  a 
given  time,  the  delay  line  is  programmed  to  move  to  a  given 
position.  Thereafter,  the  optical  delay  is  changed  at  a  rate 
determined  by  the  programmed  velocity.  The  delay  line  will  move 
at  this  rate  until  the  programmed  position  or  velocity  is  changed. 

The  composite  three-stage  servo  system  has  a  moderately  high 
bandwidth  when  the  error  signal  is  small  (i.e.,  when  delay  line 
position  minus  target  position  is  less  than  about  4  pm).  However, 
when  switching  from  one  star  to  the  next,  the  programmed 
position  can  suddenly  change  by  15m.  Since  the  PZT  and  voice 
coil  have  small  dynamic  ranges,  and  in  addition  have  limits  on 
their  velocity  and  acceleration,  these  servo  systems  will  not  behave 
properly  in  response  to  such  a  large  step.  Left  uncorrected,  the 
resultant  oscillations  are  sufficiently  violent  to  misalign  optics. 
Therefore,  in  addition  to  the  linear  trucking  mode  of  operation 
described  above,  the  delay  line  has  several  other  modes  of 
operation.  A  significant  portion  of  the  delay  line  servo  code  detects 
input  conditions  which  would  require  a  change  in  the  servo  mode 
so  that  the  delay  line  remains  well  behaved  under  all  conditions. 

Table  2  shows  the  various  modes  of  operation  of  the  three 
delay-line  loops.  The  modes  of  the  various  servos  depend  on  the 
inputs  to  the  delay  line.  Rapid  changes  in  the  programmed 
position  or  velocity  is  the  usual  reason  for  changes  in  operating 
mode.  Hysteresis  in  the  transition  thresholds  prevents  oscillations 
caused  by  changes  in  the  servo  parameters  between  modes.  Inputs 
from  the  limit  switches  can  also  cause  changes  in  the  servo  modes. 
For  example,  two  limit  switches  are  used  at  each  end  of  the  track. 
The  first  informs  the  computer  that  the  delay  line  should  start 
decelerating  in  order  to  prevent  the  delay  line  from  moving 
beyond  the  second  limit  switch,  which  disables  the  motors.  Other 
interlocks  prevent  the  laser  fringe  counter  from  running  at  2 
million  fringes  per  second,  taking  the  delay  line  with  it.  if  the  laser 
becomes  misaligned.  For  a  more  detailed  description  of  the  delay 
line,  the  reader  is  referred  to  Holm  (1986). 


2.4.  Stellar  fringe  detector 

With  proper  operation  of  the  star  tracker  and  delay  line,  the  two 
beams  combined  at  the  beamsplitter  should  interfere.  This  section 


Table  2.  Inputs:  target  position,  target  velocity,  current  laser 
positions  (2),  limit  switches  (8),  voice  coil  position 


Servo 

Modes 

PZT 

Linear  track 

Off 

Voice  coil 

Linear  track 

Damp 

Off 

Motor 

Linear  track 

Accel./deoel. 

Constant  vel.  Off 

describes  the  procedure  by  which  the  amplitude  and  phase  of  the 
interference  pattern  are  measured. 

In  the  Mark  III  interferometer,  the  interfering  wavefronts  are 
parallel.  Hence,  the  two  outputs  of  the  beamsplitter  are  uniform 
over  the  aperture  and  complementary.  If  the  phase  difference  is  90 
degrees,  all  then  light  will  go  to  one  of  the  two  outputs;  if  the  phase 
difference  is  —90  degrees,  the  light  will  go  to  the  other  output. 

The  fringe  detection  algorithm  has  been  described  previously 
(Shao  and  Staelin.  1977).  The  basic  technique  uses  path  length 
modulation  and  synchronous  demodulation.  The  PZT  controlled 
mirror  in  one  of  the  vacuum  delay  lines  is  modulated  with  a  S(K)- 
Hz  triangle  wave.  If  the  path  difference  in  the  two  arms  of  the 
interferometer  is  within  the  coherence  length  of  the  light  (a  few 
wavelengths),  the  intensity  of  light  at  the  output  of  the  beamsplit¬ 
ter  will  be  as  shown  in  Fig.  6.  During  one  sweep  of  the  PZT.  the 
detected  photons  are  grouped  into  four  time  bins  denoted  by  /4.  B. 
C.  and  D.  The  amplitude  and  phase  of  the  fringe  are  given  by 
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In  the  above  equations,  the  C  -  /4  and  D  —  B  terms  are  the  sine  and 
cosine  components  of  the  fringe. 

The  measured  phase  is  unwrapped  so  that  if  the  atmosphere 
moves  the  fringe  position  beyond  ±  a.  the  correct  error  signal  will 


be  used  to  track  the  fringe.  The  unwrapped  phase  is  u.sed  as  the 
error  signal  in  a  first-order  servo.  As  mentioned  before,  the  delay 
line  is  programmed  by  position  and  velocity.  The  unwrapped 
phase  times  a  gain  term  is  added  to  the  previous  programmed 
position  of  the  delay  line  to  yield  the  new  programmed  position. 
The  programmed  velocity  is  the  sidereal  rate.  The  closed-loop 
bandwidth  of  this  servo  is  about  20  Hz.  A  similar  fringe  tracker 
was  used  on  the  Mark  I  and  Mark  II  interferometers  to  track 
stellar  fringes. 

The  fringe  tracking  servo,  implemented  entirely  in  software,  is 
the  simplest  part  of  the  stellar  fringe  detection  subsystem.  Other 
software  components  of  the  fringe  detection  subsystem  perform 
functions  such  as  stellar  fringe  search  and  initial  acquisition, 
fringe  recovery  after  tracker  failure,  central  fringe  identification, 
and  automatic  transition  between  operating  modes.  In  addition  to 
the  above  software  components,  other  system  components  that 
affect  the  operation  of  the  stellar  fringe  detector  are  the  PZT  stroke 
servo,  and  the  optics  for  two  color  fringe  measurements  and  for 
accurate  visibility  measurements. 

The  instrument  was  designed  to  make  accurate  astrometric 
measurements  and  accurate  stellar  diameter  measurements  on  as 
many  stars  as  possible.  The  design  of  a  stellar  fringe  tracker  is 
complicated  by  the  fact  that  the  requirements  of  high  tracker 
sensitivity,  low  systematic  astrometric  error,  and  low  visibility 
errors  are  not  entirely  consistent. 

The  stroke  of  the  PZT  can  be  varied  by  one  of  the  control 
computers.  Since  the  piezoelectric  constant  is  temperature  de¬ 
pendent.  the  peak-to-peak  amplitude  of  the  PZT  triangle  wave  is 
monitored  with  a  laser  interferometer  and  the  drive  voltage  is 
adjusted  to  maintain  the  proper  amplitude  to  within  0.5  nm.  For 
stellar  diameter  measurements,  the  stroke  is  set  at  the  effective 
wavelength  of  the  narrow  bandpass  filter.  For  astrometry,  the 
stroke  is  set  at  the  effective  wavelength  of  the  red  channel. 

When  the  interferometer  is  set  to  point  at  a  star,  the  delay  line 
is  programmed  to  go  to  the  expected  position  of  the  fringe  plus 
60  pm.  From  past  experience,  the  calculated  fringe  position  is 
never  more  than  30  pm  in  error  if  the  baseline  vector  used  in  the 
calculation  is  derived  from  stellar  fringe  data  less  than  a  week  old. 
The  delay  line  is  programmed  to  scan  for  the  fringe  at  approxi¬ 
mately  60pm  per  second.  The  fringe  amplitude  is  calculated  every 
atmospheric  coherence  time,  and  if  a  fringe  is  detected,  the  scan  is 
stopped  and  the  servo  turned  on. 

When  the  tracker  loses  the  fringe,  the  delay  line  is  moved  back 
Spm  and  a  search  is  restarted.  After  initial  fringe  acquisition,  the 
search  space  is  restricted  to  a  20-pm  interval  centered  about  the 
most  recent  position  where  the  fringe  was  found.  All  of  these 
search  parameters  are  superimposed  on  the  fringe  motion  due  to 
earth  rotation,  which  is  calculated  using  a  baseline  vector  typically 
derived  from  the  prei  ous  night's  data.  When  the  seeing  is 
exceptionally  good,  i.e..  sub-arc-second  seeing,  the  fringe  tracker 
will  stay  locked  on  one  fringe  for  many  minutes,  for  up  to  l(t' 
atmospheric  coherence  times. 

If  the  fringe  stays  locked  for  more  than  one  second,  a  central 
fringe  identification  algorithm  is  enabled.  The  visibility  during 
that  one-second  interval  is  compared  with  a  threshold.  If  the 
visibility  is  below  the  threshold,  it  is  assumed  that  the  tracker  is  on 
the  wrong  fringe  and  the  tracker  will  jump  one  fringe  to  the  right 
or  left.  The  direction  of  the  tracker  jumps  is  determined  by  the 
visibility  gradient.  Since  the  tracker  has  an  rms  phase  error  of 
about  1/6  of  a  fringe,  the  visibility  gradient  can  be  determined  by 
looking  at  the  distribution  of  visibility  versus  phase.  Since  the 
central  fringe  has  the  highest  visibility,  the  tracker  will  usually  hop 
in  that  direction. 
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The  operation  of  the  various  systems  which  comprise  the  star 
light  fringe  detector  requires  that  a  number  of  quantities:  phase, 
amplitude.  SNR.  visibility  gradient,  etc.  be  calculated  in  real  time. 
Although  the  system  is  rather  complicated,  the  resulting  tracker  is 
extremely  robust.  The  interferometer  can  be  commanded  to  point 
at  a  star  and  the  tracker  will  lock  onto  the  fringe  within  90s.  'Thus, 
the  software  complexity  was  introduced  to  enable  the  instrument 
to  operate  efficiently  under  a  wide  variety  of  conditions  with  a 
minimum  of  operator  intervention. 

2.5.  Laser  melralogy 

During  initial  operation  in  1986.  the  siderostat  laser  metrology 
system  was  not  operational.  However,  the  laser  system  is  impor¬ 
tant  for  astrometry  and  should  become  operational  in  the  summer 
of  1987.  The  basic  laser  interferometer  used  to  measure  changes  in 
instrument  geometry  was  described  in  the  delay  line  section  of  this 
paper.  Two  such  interferometers  measure  the  position  of  the  cat's 
eye  retroreflectors.  An  additional  12  interferometers  will  be  used 
to  measure  the  position  of  the  three  siderostat  mirrors. 

The  metrology  system  measures  the  motion  of  the  pivot  point 
of  the  siderostat  by  measuring  the  distance  between  four  reference 
corner  cubes  in  an  invar  plate,  which  is  bonded  to  the  concrete 
pedestal,  and  the  surface  of  a  zerodur  sphere  attached  to  the  mirror 
as  shown  in  Fig.  7.  The  sphere  is  made  accurate  to  /.  10  and  the 
center  of  the  sphere  is  coincident  with  the  mechanical  pivot  point 
to  20pm.  The  interferometers  have  5nm  resolution;  the  desired 
accuracy  of  the  pivot  point  monitor  is  100  nm  for  one-milliarcsec 
astrometry  with  a  20-m  baseline.  The  measurements  from  these 
lasers  will  be  used  to  correct  for  mechanical  errors  in  the  siderostat 
as  well  as  thermal  motion  of  the  siderostat.  as  described  below. 

The  basic  model  for  the  instrument  is  given  by 

Fringe  Pos  =  5  (8) 

S  is  the  unit  vector  to  the  star. 

B  is  the  baseline  vector, 

C  is  the  delay  offset. 

If  the  sidcrostats  were  perfect,  the  baseline  vector  would  be  the 
vector  that  joined  the  pivot  points  of  the  two  siderostats.  If  the 
delay  lines  were  commanded  to  move  to  their  zero  positions,  the 
delay  offset  C  Would  be  the  difference  in  the  optical  paths  from  the 
beamsplitter  to  the  pivot  points  of  the  two  siderostats.  Motion  of 
the  sidero.  .at  pivot  points  affects  both  B  and  C.  Motion  of  the 
internal  optics  only  affects  C. 

With  our  12-m  baseline  on  Mt.  Wilson,  the  rms  atmospheric 
fringe  motion  was  6pm  under  conditions  of  good  seeing.  For 
stellar  diameter  measurements  with  the  fringe  tracker  running, 
neither  instrument  drifts  nor  atmospheric  motion  is  important. 
However,  for  astrometric  measurements,  the  measured  fringe 
position  in  Eq.  t8)  is  made  relative  to  the  baseline  vector  and  delay 
offset.  Thus,  it  is  i.ecessary  that  both  0and  C  be  constant  over  a 
period  of  time  that  is  sufficiently  long  for  the  interferometer  to  be 
able  to  observe  the  fringe  position  of  a  number  of  stars. 

The  parameters  Band  C'can  change  for  two  reasons;  imperfect 
siderostat  bearings  and  thermal  drift.  Ideally,  the  siderostat  pivots 
about  a  point  on  the  surface  of  the  mirror.  In  practice,  the  mirror 
translates  as  it  rotates  in  a  rather  unpredictable  way  due  to 
irregularities  in  the  bearings  and  flexure  of  the  mirror  cell.  The 
laser  metrology  system  is  designed  to  measure  these  imperfections. 
The  result,  in  theory,  is  that  baseline  drifts  will  be  limited  to 
thermal  drift  of  the  top  of  the  concrete  pier  to  which  the  reference 
corner  cubes  in  the  laser  system  are  attached.  With  rapid  switching 
between  stars,  only  baseline  drift  with  a  time  scale  faster  than  the 
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Fi«-7.  Laser  metrology  system  lor  baseline  seetor  monitor 


SIDEROSTAT  WITH  4-LASER  PLATE 


switching  cycle  results  in  astrometric  errors  for  large  angle  relative 
astrometry. 

In  order  to  monitor  the  position  of  the  pivot  point,  a  minimum 
of  three  interferometers  is  necessary.  The  fourth  interferometer 
provides  the  redundancy  necessary  to  solve  for  the  position  of  the 
reference  corner  cubes.  The  interferometers  measure  the  position 
of  the  center  of  the  sphere.  Certain  parameters,  such  as  the 
distance  between  the  center  of  the  sphere  and  the  surface  of  the 
mirror  must  be  derived  from  stellar  data. 

The  mechanical  accuracy  of  the  instrument  without  the  laser 
metrology  system  is  approximately  10  IS  pm  if  the  baseline  was 


calculated  from  the  previous  night's  stellar  data.  The  .size  of  this 
error  implies  that  searching  for  the  fringe  would  take  only  a  few 
hundred  milliseconds,  an  insignificant  amount  of  time  compared 
to  the  .30  to  60  s  needed  for  the  siderostats  to  slew  to  a  new  star. 
The  mechanical  errors  are  roughly  equivalent  to  the  errors  due  to 
atmospheric  turbulence,  or  to  the  errors  in  the  primary  mirror  of  a 
high  quality  large  astronomical  telescope.  With  the  laser  system 
we  hope  to  improve  this  accuracy  by  two  orders  of  magnitude. 

While  this  higher  mechanical  accuracy  is  not  necessary  for 
stellar  diameter  measurcmeni.s.  it  is  crucial  for  ultra-high  accuracy 
astromet  ).  This  level  of  accuracy  is  also  necessary  for  future 
space-based  pha.se  coherent  long  baseline  interferometers. 


3.6.  Sulcrosuil  lonHol 

The  siderostat  control  subsystem  is  in  many  ways  similar  to  the 
computer  controlled  pointing  system  for  a  modern  astronomical 
telescope.  The  siderostat  itself  is  a  two  axis  motorized  gymbaled 
mirror  mount  <Fig.  8).  with  the  fixed  axis  tilted  75  from  vertical 
towards  the  West,  The  siderostats  u.se  conventional  precision  worm 
and  wheel  gears  and  microstepped  stepper  motors.  The  control 
computer  sets  programmable  pulse  generators  to  control  motor 
velocity.  The  average  step  size  is  0'.'02  to  provide  smooth  tracking, 
and  the  maximum  slew  rate  is  1 .8  per  second.  Typical  backlash  in 
the  gear  train  is  less  than  1". 

The  control  system  performs  several  tasks.  One  is  open  loop  or 
blind  pointing,  where  a  geometric  model  of  the  siderostat  is  used 
to  point  the  siderostat  and  track  the  star.  A  second  task  is  closed 
loop  pointing,  where  the  siderostat  is  moved  to  keep  the  high¬ 
speed  PZT  til!  mirrors  in  the  center  of  their  range.  The  third  task  is 
updating  the  parameters  of  the  geometric  model  when  a  star  is 
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TaMe  3.  Siderostat  model  parameters 
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Angle  between  fixed  and  movable  axis 
Angle  between  movable  axis  andmirror  normal 
Zero  point  of  motor  position  on  the  two  axes 


acquired.  In  addition,  the  siderostat  control  system  controls  the 
sequencing  of  observations  during  the  night. 

The  geometry  of  the  siderostat  described  earlier  was  that  of  an 
ideal  siderostat.  However,  the  stationary  axis  may  not  be  exactly 
1 5  degrees  from  horizontal;  it  may  not  be  pointed  due  west,  and  the 
two  axes  may  not  be  exactly  perpendicular.  Altogether,  there  are 
eight  parameters  in  the  geometric  model  of  the  siderostat.  With 
values  of  these  eight  parameters,  there  is  a  closed-form  relation 
between  the  angular  positions  of  the  motors  and  the  pointing  of 
the  instrument.  It  is  also  possible  to  invert  this  problem  i  closed 
form.  By  pointing  the  siderostat  at  a  number  (  >  4)  of  bright  stars 
whose  coordinates  are  known,  and  recording  the  positions  of  the 
motors  when  the  stars  are  centered  in  the  field  of  view  of  the 
interferometer,  the  eight  parameters  of  the  model  can  be  deter¬ 
mined  in  a  least  squares  fashion.  In  our  initial  setup,  we  used  the 
FK4  coordinates  of  a  half  dozen  stars  precessed  to  the  night  of 
observation. 

Since  the  siderostats  direct  the  starlight  into  the  central  beam¬ 
combining  building,  the  alignment  of  the  interferometer  optics 
inside  the  building  will  affect  the  pointing  of  the  siderostat.  The 
internal  alignment  of  the  interferometer  is  represented  by  two  of 
the  eight  parameters  describing  the  siderostat.  All  of  the  para¬ 
meters  are  listed  in  Table  3.  No  attempt  was  made  to  model  errors 
in  the  precision  worm  and  wheel  gear  or  to  model  flexure.  It  is 
most  likely  that  the  residual  errors  are  primarily  in  the  gears 
Residuals  from  a  least-squares  fit  for  the  siderostat  parameters  are 
of  order  5".  Typical  pointing  errors  for  stars  which  arc  not  used  in 
the  least-squares  fit  are  approximately  7"  to  10". 

However,  each  new  star  that  is  observed  provides  some 
information  on  the  eight  siderostat  parameters.  Hence,  rather 
than  an  off-line  least-squares  solution,  a  real  time  algorithm  could 
be  implemented  to  update  the  eight  parameters  every  time  the 
instrument  locked  onto  a  new  star.  In  the  Mark  III.  a  Kalman 
filter  is  used  to  update  the  parameter  vector.  In  operation,  the 
initial  values  of  the  parameters  at  the  beginning  of  the  night  are  the 
parameters  derived  from  the  previous  night's  observations.  These 
parameters  are  updated  for  every  new-  star  and.  under  normal 
conditions,  the  initial  values  would  be  essentially  lost  after 
approximately  6  stars.  Thereafter,  the  parameter  estimates  would 
reflect  thermal  drift  of  the  instrument  during  the  night. 

This  complex  siderostat  model  is  used  only  for  initial  acqui¬ 
sitions.  After  the  star  tracker  locks  onto  the  star,  the  positions  of 
the  PZT  tilt  mirrors  are  sent  to  the  siderostat  control  computer. 
The  siderostat  drive  rates  for  the  motors  are  then  modified  to  keep 
the  PZTs  in  the  center  of  their  range.  Since  the  PZT  tilt  mirrors  in 
effect  represent  realignment  of  the  internal  optics,  the  corrections 
to  the  sideral  rate  are  generated  by  changing  the  parameters 

Of,  0p . 

The  siderostat  control  computer  also  calculates  the  theoretical 
fringe  position  and  velocity  for  the  optical  delay  line.  These  values 


are  sent  to  the  delay  line  control  computer  through  a  shared 
multiport  memory  every  2  s.  The  theoretical  fringe  position  is 
calculated  from  a  baseline  vector  derived  from  the  previous  night's 
data,  in  the  future,  we  plan  to  update  the  baseline  vector  with  a 
Kalman  filter  in  the  same  way  the  siderostat  parameters  are 
updated.  On  the  first  night  of  interferometric  observations,  the 
baseline  vector  was  determined  by  conventional  surveying  techni¬ 
ques.  The  surveying  error  was  approximately  1  mm.  or  about  20" 
of  baseline  tilt. 

The  last  function  of  the  siderostat  control  computer  is  the 
sequencing  of  observations.  An  abbreviated  version  of  the  FK4 
catalog  is  kept  on-line  (In  1987.  the  FK  5  catalog  will  be  kept  on¬ 
line).  A  menu  driven  program  enables  the  operator  to  list  and  plot 
the  stars  visible  at  the  present  time  and  to  easily  define  the  stellar 
observation  sequence.  The  siderostat  control  computer  also 
communicates  with  the  other  two  control  computers  to  enable 
fully  automatic  operation  once  the  observing  sequence  is  pro¬ 
grammed.  The  siderostat  control  computer  signals  the  star  tracker 
computer  to  initiate  a  star  acquisition  sequence.  When  the  star 
tracker  has  locked  onto  a  star  it  signals  the  fringe  tracker 
computer  to  initiate  a  fringe  acquisition  sequence.  When  the 
allotted  observing  time  for  a  star  is  over,  the  siderostat  control 
computer  signals  the  other  two  computers  to  disengage  their  servo 
systems.  For  additional  information,  the  reader  is  referred  to 
Kelly  (1986). 


3.  Operational  characteristics 

The  instrument  was  built  and  debugged  over  a  33  month  period 
Fringes  were  observed  on  the  first  night  of  operation.  September 
12.  1986.  Over  the  following  months,  observational  techniques 
were  refined  and  automated.  Figure  9  is  a  plot  of  the  observations 
made  on  1986.  November  II.  The  vertical  axis  is  the  delay  line 
position  where  fringes  were  observed.  Each  dot  represents  a  l(K)-s 
observation.  A  total  of  165  observations  of  23  stars  were  made 
over  a  I0.5-h  period.  .Stars  over  a  declination  range  of  68  and  an 
hour  angle  range  of  8  h  were  observed.  All  these  observations  were 
of  stars  3.0  magnitude  or  brighter. 

The  initial  measurements  in  1986  w  ill  be  discussed  in  a  senes  of 
separate  papers.  Wide-angle  one-color  astrometry  show  s  night-to- 
night  repeatability  of  about  0'.’()4 1  Mozurkewich  et  al..  1987).  and 
the  formal  error  for  seven  nights  of  data  is  0"02  in  declination 
Significant  improvements  are  anticipated  w  hen  two-color  analv  sis 
is  applied  ar  J  when  the  laser  siderostat  monitor  is  installed  in  the 
summer  of  1987  Preliminary  measurements  on  the  atmospheric 
limitations  of  two  color  absolute  astrometric  measurements 
indicate  that  single  observations  of  ~  100  s  should  yield  accuracies 
of  about  ()?0I.  a  factor  of  5  improvement  over  one-color 
observations  (Colavity  et  al  .  1987).  Two-color  observations 
should  approach  0''003  for  a  single  night's  observation.  Last  of  all. 
we  have  made  preliminary  stellar  diameter  measurements  of 
pGem  and  xCMi.  Assuming  a  uniform  disk,  we  measure 
6.15  +  0.25  mas  for  zCMi  and  8.70  +  0.15  mas  for  ^Gem.  The 
uniform  disk  diameter  for  xCMi  is  not  in  agreement  with  prior 
measurements  (Hanbury  Brown  et  al..  1974)  in  part  due  to  the 
difference  in  observation  wavelength  and  uniform  disk  assimip- 
tion,  A  much  more  careful  discussion  of  accurate  stellar  diameter 
measurements  is  given  elsewhere  (Shao  et  al..  1988). 

The  instrument  was  never  operated  near  the  limiting  magni¬ 
tude  of  the  fringe  tracker  because  the  temporary  star  tracker  we 
were  using  could  not  track  stars  fainter  than  3  mag.  The  tracker 
using  the  photon  camera  should  be  capable  of  operation  on  faint 
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Observed  Oeloys  1  1  November  1986 


UT  (seconds) 

Position  of  delay  line  where  fringes  were  observed,  165  observations  of  23  stars 


sources,  to  12tnag  under  good  seeing  conditions  when  only 
mechanical  errors  need  to  be  tracked.  When  the  new  star  tracker  is 
installed  in  the  spring  of  1987.  the  limiting  magnitude  will  be  set  by 
the  fringe  detection  system.  We  plan  to  operate  the  interferometer 
in  two  distinct  modes  that  have  two  very  different  magnitude 
limits.  The  two  modes  represent  two  types  of  interferometers.  One 
is  a  closed  loop  phased  interferometer.  The  second  is  an  open  loop 
coherent  interferometer  or  absolute  interferometer.  We  plan  to 
test  this  second  mode  of  operation  in  the  next  two  years. 

3.].  Phase-coherent  interferometry 

The  Mark  III  was  designed  to  be  an  active  phased  interferometer. 
In  this  mode  of  operation,  the  fringe  tracker  compensates  for  ail 
optical  path  fluctuations,  including  atmospheric  turbulence  to 
~/l/6rms.  Fringe  phase  measurements  in  two  colors  are  used  for 
astrometric  purposes,  while  light  from  a  bandpass  filter,  <  400  A. 
is  used  for  stellar  diameter  measurements.  It  is  also  possible  to 
measure  the  diameter  of  a  star  at  emission  or  absorption  lines.  For 
diameter  measurements  in  narrow  spectral  lines,  we  make  use  of 
the  phase  coherent  property  of  an  active  interferometer.  For 
narrow  spectral  lines,  the  photon  flux  in  the  line  might  be  10'  ’  of 
the  total  flux  and  the  number  of  photons  per  atmospheric 
coherence  time  may  be  as  low  as  10'^.  However,  since  the  active 
optics  system  is  phase  coherent,  the  coherent  integration  time  is 
not  limited  to  the  atmospheric  coherence  time.  By  choosing  a 
longer  coherent  integration  time,  >  t  $.  a  fringe  visibility  measure¬ 
ment  can  be  made  in  10*  s  with  a  signal-to-noise  ratio  of  100. 

For  astrometry,  the  fringe  phase  and  the  position  of  the  optical 
delay  line  is  recorded  when  the  fringe  tracker  is  locked  on  the 


central  fringe.  In  this  case,  the  astrometric  error  due  to  photon 
noise  is  negligible  and  the  major  errors  are  instrumental.  A 
detailed  discussion  of  one-  and  two-color  astrometry  will  be 
presented  in  separate  papers  (Colavita  et  al.,  1987;  Mozurkewich 
etal.,  1988).  However,  the  astrometric  capability  of  the  instrument 
affects  all  other  modes  of  operation  of  the  instrument. 

With  the  known  coordinates  of  FK5  stars  it  is  possible  to 
derive  the  baseline  vector  and  delay  offset  by  measuring  the  fringe 
position  of  four  stars.  This  derived  baseline  is  as  accurate  as  the 
coordinates  of  the  FK  5  stars  used.  An  accurate  baseline  vector  in 
turn  enables  the  instrument  to  find  the  fringe  faster  and  to  increase 
the  number  of  observations  from  several  per  night  to  several 
hundred  per  night.  This  astrometric  capability  is  the  key  to 
operation  of  the  instrument  as  an  absolute  interferometer,  de¬ 
scribed  next. 

3.2.  Absolute  interferometer 

Within  the  next  year  we  plan  to  install  the  software  for  real  time 
baseline  calibration.  In  this  mode,  the  interferometer  will  rapidly 
switch  between  4-6  FK  5  stars  to  calibrate  the  baseline,  and  then 
point  blind  at  a  faint  astronomical  object  for  about  10^  s.  This 
cycle  will  be  repeated  to  yield  the  desired  integration  time  on  the 
faint  object.  Three  types  of  measurements  are  planned,  two  of 
them  in  the  near  future.  There  are  significant  advantages  to  phase 
coherent  operation  and  operation  in  the  absolute  mode  is  reserved 
for  objects  too  faint  for  the  fringe  tracker,  9mag  for  10cm  aper¬ 
tures.  All  three  types  of  measurements  planned  for  the  absolute 
interferometer  are  photon  starved  measurements.  Under  these 
conditions,  direct  phase  measurements  are  not  possible  and  the 
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only  quantities  that  can  be  measured  are  amplitude,  group  delay, 
and  closure  phase.  Photon  starved  conditions  occur  when  the 
average  number  of  photons  per  r<,  aperture,  per  tg  integration 
time,  is  less  than  1  in  the  coherent  optical  bandpass. 

In  open  loop  operation,  a  1 2-m  interferometer  will  exhibit  rms 
phase  fluctuations  of  approximately  20/.  at  0.6  pm  for  1"  seeing. 
The  fringe  detector  must  be  able  to  accomodate  this  lack  of 
temporal  coherence.  On  the  other  hand,  to  maximize  sensitivity, 
light  from  all  wavelengths  should  be  used  for  interferometric 
measurements.  The  basic  fringe  detector  (Fig.  10)  consists  of  a 
prism  spectrometer  with  an  imaging  photon  counting  detector. 
The  required  resolution  of  the  spectrometer  is  of  the  order  of 


where  X  is  the  wavelength  of  observation,  L  the  length  of  the 
baseline,  and  rgoeX*'*  is  equal  to  about  lOcm  at  0.6pm  for  I" 
seeing.  Conservatively,  SO  A  channels  are  assumed  in  the  calcu¬ 
lations  below.  In  the  following  sections  we  briefly  describe  the 
signal  processing  of  the  output  of  the  detector  for  the  three 
different  types  of  photon  starved  measurements. 

3.2.1.  Photon-starved  amplitude  measurements 

The  estimator  for  fringe  amplitude  must  be  modified  to  take  into 
account  a  bias  that  is  important  at  low  photon  rates.  An  unbiased 
estimator  for  visibility  squared  is  given  by 

2  L  <Ar-Dark>‘  J  '  ' 

where  Af^/t  +  B+C+D^  total  photon  flux.  Dark  =  dark 
count  and  <  >  denotes  time  average. 


The  signal-to-noise  ratio  of  the  above  visibility  estimator  has 
been  extensively  studied  both  for  amplitude  interferometry  and 
speckle  interferometry.  There  are  two  regions  of  interest,  depend¬ 
ing  on  photon  rate.  The  dependence  of  integration  time  needed  to 
achieve  a  given  signal-to-noise  ratio,  given  a  photon  flux  N  per  rg 
aperture,  per  Zg  integration  time,  per  spectral  channel  is  given  by 

Integration  time  x  Af '  ‘ ,  yv  >  1 ,  N>  Dark , 

xAf’M  >  Af>Dark.  (II) 

For  stellar  diameter  measurements,  the  numerator  and  de¬ 
nominator  of  the  expression  in  Eq.  (10)  must  be  evaluated  every 
to  ms  for  each  spectral  channel,  typically  60  channels  of  SO  A 
bandwidth  fora  12-m  baseline.  Ifcollecting  apertures  greater  than 
Zg  (10cm)  are  used,  there  must  be  a  spectrometer  for  each  rg 
subaperture.  In  photon  starved  opteration  there  is  no  hard 
magnitude  limit  as  there  is  in  active  phase  coherent  operation.  As  a 
practical  limit  we  require  a  signal-to-noise  ratio  of  30  in  I0*s  of 
integration.  Assuming  an  unresolved  star,  the  required  photon 
flux  is  1.6  10'^  photons  in  each  of  60  spectral  channels,  per  rg 
aperture  per  10  ms.  This  corresponds  to  a  magnitude  limit  of  14.3 
for  1 0-cm  apertures.  With  1  -m  apertures,  the  limiting  magnitude  is 
about  16.8  (see  Appendix).  Accurate  visibility  measurements  will 
require  much  higher  signal-to-noise  ratios  and  the  magnitude  limit 
will  degrade.  The  optical  efficiency  of  the  spectrometer  was 
assumed  to  be  100%  and  the  quantum  efficiency  of  the  imaging 
photon  counter  was  assumed  to  be  10%.  For  these  reasons,  these 
estimates  should  be  considered  optimistic. 

3.2.2.  Photon-starved  group  delay  astrometry 

The  signal  from  the  spectrometer  can  be  processed  in  a  different 
way  to  yield  astrometric  information  at  faint  magnitudes.  When 
the  optical  path  difference  in  the  two  arms  of  the  Interferometer  is 
different  by  many  wavelengths,  say  3  pm,  there  will  be  construc¬ 
tive  interference  at  A  =  0.5  pm  and  0.6  pm.  and  destructive  inter¬ 
ference  in  between.  More  precisely,  the  spectrometer  output  will 
exhibit  a  fringe  pattern  due  to  the  interference  of  light  from  the 
two  siderostats  as  given  by 

/(v)  = /g(v)  (1  +  sin 2n  vr) .  (12) 

where  /„  ( >■)  is  the  spectrum  of  the  star,  v  is  optical  frequency  ( 1 .  /.) 
and  T  is  the  optical  path  difference  in  the  two  arms  of  the 
interferometer,  i.e..  the  group  delay.  Equation  (12)  has  a  sine 
rather  than  a  cosine  because  the  Mark  HI  uses  a  dielectric 
beamsplitter  to  combine  the  beams.  The  Fourier  transform  of  the 
"fringed"  spectrum  /(v)  will  have  a  delta  function  at  r.  At  this 
delay,  the  Fourier  transform  has  both  an  amplitude  and  phase.  In 
the  photon  starved  mode,  where 

J/(v)r/v<l  pholon/rg/lOms,  (13) 

the  phase  is  not  measurable. 

The  basic  signal  processing  algorithm  is  to  average  the  power 
spectrum  of  /(r).  As  in  the  previous  section,  each  ig  subaperture 
must  be  separately  processed  if  the  collecting  apertures  are  large. 
Since  r  changes  due  to  turbulence,  a  time  average  of  the  power 
spectrum  is  a  gaussian  distribution  of  delta  functions.  For  a  1 2-m 
baseline,  the  standard  deviation  of  this  gaussian  is  '-12  pm  for 
r  seeing.  A  SNR  of  30  in  this  case  means  that  the  peak  in  the 
gaussian  is  approximately  30  times  the  rms  fluctuations  in  the  rest 
of  the  power  spectrum.  For  lO-cm  apertures,  the  magnitude  limit 
is  approximately  14.3  mag.  For  1-m  apertures,  the  corresponding 
limit  is  16.8  (see  Appendix). 
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3.2.3.  Photon-starved  closure  phase 

Closure  phase  measurements  were  first  made  with  radio  inter¬ 
ferometers  (Jennison,  19S8),  and  subsequently  suggested  for 
optical  interferometers  (Rogstad.  1968).  Recently,  closure  phase 
measurements  were  made  at  visible  wavelengths  by  masking  the 
aperture  of  a  large  telescope  (Baldwin.  1986;  Readhead.  1988)  at 
photon  rich  levels.  In  closure  phase  measurements,  three  inter¬ 
ferometers.  whose  baselines  form  a  closed  triangle,  measure  three 
complex  visibilities  corrupted  by  turbulence.  However,  the  closure 
phase,  the  sum  of  the  phases  of  the  three  complex  visibilities,  is 
independent  of  the  atmosphere.  In  photon  starved  measurements, 
raw  phases  cannot  be  measured.  Instead,  the  product  of  the 
complex  visibilities  is  averaged,  i.e.. 

triple  product  =  <,-1,  e'*'' /f.e'"’  /tje"').  (14) 

where  /4,  e'”'  is  the  complex  visibility  from  interferometer  1.  The 
importance  of  averaging  the  complex  visibilities  comes  from  the 
fact  that  the  product  of  complex  visibilities  is  not  biased  in  the 
photon  starved  region.  The  closure  phase  is  simply  the  imaginary 
part  of  the  log  of  the  triple  product  (see  Appendix).  The  triple 
product  was  first  used  in  high  resolution  optical  imaging  for 
reconstruction  of  speckle  images  (Weigelt.  1977).  a  technique 
called  speckle  masking.  This  technique  extends  the  sensitivity  of 
closure  phase  measurements  more  than  an  order  of  magnitude 
beyond  conventional  photon  rich  measurements. 

For  long  ba.seline  interferometers,  the  triple  product  must  be 
calculated  every  1 0ms.  for  each  spectral  channel,  and  each  r„ 
aperture.  The  signal-to-noise  ratio  (SNR)  of  the  triple  product  is 
qualitatively  different  from  the  SNR  for  amplitude  measure¬ 
ments.  At  high  photon  fluxes  .V>1.  the  SNR  of  all  these 
quantities  varies  as  N .  At  low  photon  fluxes.  .V  <  1 .  the  SNR  for 
amplitude  measurements  varies  as  N.  and  that  of  the  triple 
product  varies  as  (see  Appendix).  Again,  assuming  an  SNR 
of  30  in  1 0'*  s  (/f  (/)  measured  to  /  1 80).  the  magnitude  limit  with  10- 
cm  apertures  is  approximately  12.5  mag.  or  a  photon  flux  of  0.08 
photons  per  r^.  per  1 0  ms.  per  50  A  bandwidth.  The  limiting 
magnitude  for  1-m  apertures  is  14.2. 

Although  three  siderostats  will  be  in  operation  in  1987.  closure 
phase  measurements  are  not  planned  for  the  immediate  future.  In 
order  to  operate  three  interferometers  simultaneously  the  second 
delay  line  must  be  converted  to  high  precision  operation.  Oper¬ 
ation  as  an  absolute  interferometer  will  then  enable  closure  phase 
measurements  to  be  made  with  a  one  to  two  order-of-magnitude 
increa.se  in  sensitivity  over  active  phase  coherent  operation. 


4.  Conclusioa 

In  the  past  seventy  years  a  number  of  stellar  interferometers  have 
been  proposed  and  built.  Starting  with  Michelson.  the  instruments 
have  gotten  progressively  more  complex  and  sophisticated 
(Koechlinetal..  1985;  Labeyrieet  al..  1986;  Davis  et  al..  1986).  The 
Mark  III  incorporates  many  of  the  ideas  proven  by  these  previous 
instruments.  The  major  differences  between  this  interferometer 
and  others  reside  primarily  in  its  operational  characteristics  it  is 
fully  automated,  it  is  easy  to  operate,  and  it  is  mechanically 
accurate  and  thermally  stable.  It  is  also  capable  of  mak  ing  literally 
hundreds  of  observations  in  one  night;  both  astrometric  and 
stellar  diameter  measurements  are  possible.  Last  of  all.  this 
instrument  has  the  potential  to  operate  as  an  absolute  inter¬ 
ferometer  which,  from  an  astronomical  point  of  \iew.  holds 
tremendous  promise  for  observation  of  faint  objects 


I 
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To  fully  exploit  the  possibilities  of  absolute  interferometry, 
and  to  construct  images  of  faint  extragalactic  objects  at  milliarcsec 
resolution,  even  more  complex  interferometers  with  large,  mul¬ 
tiple  apertures  will  be  needed.  During  the  next  several  years  the 
Mark  III  should  demonstrate  much  of  the  necessary  technology, 
as  well  as  produce  some  of  the  scientific  results  that  would 
motivate  the  construction  of  a  large  optical  array. 
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Appendix:  photon-starved  interferometry 

The  photon-starved  signal-to-noise  ratio  for  the  three  inter¬ 
ferometric  quantities;  amplitude,  group  delay,  and  closure  phase 
are  discussed  in  this  Appendix.  We  present  limiting  magnitudes 
for  these  quantities  extrapolated  from  our  experience  with  the 
Mark  III  interferometer.  Recently,  there  has  been  considerable 
interest  in  the  use  of  light  scattered  from  the  ionosphere  from  a 
pulsed  laser  in  order  to  provide  a  probe  of  atmospheric  turbulence 
to  operate  an  active-optics  system  (Foy  and  Labeyrie.  1986).  The 
application  of  this  technique  is  not  considered  in  our  calculations. 

Estimates  of  amplitude,  group  delay,  and  closure  phase  each 
have  slightly  different  magnitude  limits.  Qualitatively,  different 
effects  degrade  the  sensitivity  for  the  different  quantities.  In  all 
cases,  we  assume  that  the  fringe  detector  is  a  dispersive  spec¬ 
trometer  followed  by  an  imaging  photon-counting  detector.  We 
also  assume  that  the  resolution  of  the  spectrometer  is  a  few  times 
greater  than  what  is  needed  to  maintain  high  fringe  contrast  in  the 
presence  of  the  optical  path  fluctuations  caused  by  atmospheric 
turbulence.  Magnitude  limits  are  calculated  assuming  an  unre¬ 
solved  star,  a  12-m  baseline.  lO-cm  collecting  apertures,  a  10ms 
coherent  integration  lime,  a  60-channeI.  50  A  resolution  spec¬ 
trometer.  unit  fringe  visibility,  a  total  integration  time  of  1 0*  s.  and 
a  final  signal-to-noise  ratio  of  .30.  The  detector  dark  count  is 
neglected,  and  a  total  system  efficiency  of  0.1  is  assumed.  We 
ignore  the  reduction  in  fringe  visibility  and  hence  signal-to-noise 
ratio  due  to  the  use  of  coherent  apertures  of  10  cm  and  integration 
times  of  10  ms.  In  practice,  the  seeing  would  have  to  be  somewhat 
better  than  I"  for  these  magnitude  limits  to  be  achieved.  Scaling 
formulae  are  given  for  the  case  of  different  seeing  conditions. 

At  high  photon  fluxes,  when  the  number  of  photons  per  /„ 
aperture,  per  atmospheric  coherence  time  Tq  .  per  spectral  channel 
i/..  is  much  greater  than  one.  the  signal-to-noise  ratio  is  propor¬ 
tional  to  the  square  root  of  the  photon  flux  for  estimates  of  all 
three  interferometric  quantities.  We  define  the  quantities  N  and  M 
as 

iV=  No.  of  delected  photons  per  frame,  i.e..  per  (•„.  To-  d7. . 

W  =  No.  of  frames 
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where  D  =  diameter  of  telescope,  T  =  total  integration  time,  and 
BW=  total  optical  bandwidth  of  the  spectrometer  (assumed  0.4  to 
0.7  pm).  The  terms  photon-rich  and  photon-starved  apply  to  the 
cases  N  >  1  and  \  <  1 ,  respectively. 

The  signal-to-noise  ratio  for  fringe  amplitude  measurements 
has  been  extensively  studied  for  many  types  of  interferometers.  It 
is  usually  more  convenient  to  consider  estimates  of  rather 
than  of  V.  The  signal-to-noise  ratio  for  estimates  of  can  be 
written  (Tango  and  Twiss,  1980) 

S)V^(K^)=  1^1 (Al) 

For  the  case  of  a  fringe-scanning  interferometer  with  4  bins  per 
scan,  N  is  replaced  by  (8/a^)Af.  The  qualitative  feature  of  this 
expression  is  that  at  low  photon  fluxes.  A/  <  1 ,  the  signal-to-noise 
ratio  is  proportional  to  A,  rather  than  to  Thus,  if  a  star  is 
dimmer  by  a  factor  of  10,  the  number  of  frames  M  must  increase 
by  a  factor  of  100  to  give  the  same  signal-to-noise  ratio. 

The  noise  behavior  of  the  group-delay  estimate  is  similar  to  the 
estimate.  The  “fringed”  spectrum  as  detected  by  the  spec¬ 
trometer  is 

/(‘')  =  ^o('')(l +sin2nvT'),  (A2) 

where  /o(v)  is  the  spectrum  of  the  star,  r  is  optical  frequency,  and  r' 
is  the  group  delay.  The  group  delay  changes  due  to  atmospheric 
turbuience.  Over  long  time  intervals,  r'  is  a  gaussian-distributcd 
random  variable  with  mean  r, ,  determined  by  the  position  of  the 
star  and  optical  delay  line,  and  variance  (t}  .  which  depends  on  the 
seeing  and  baseline  length.  The  variance  is  given  by 

<r»  =  6.88(2/2»)2(Z,/ro)’'\  (A3) 


where  L  is  the  baseline  .ength  and  is  the  atmospheric  coherence 
diameter. 

Each  coherence  time  Tq.  the  power  spectrum.  P{t).  of  /(f)  is 
calculated.  The  power  spectrum  is  then  averaged  for  the  length  of 
the  observation.  There  are  three  important  features:  a  peak  at  the 
origin  attributable  to  the  Fourier  transform  of  /o(f),  a  uniform 
bias  due  to  photon  noise,  and  a  gaussian  distribution  at  the  mean 
group  delay  r,  due  to  atmospheric  turbulence.  The  averaged 
power  spectrum  is  given  by 

</*(t)>  =  <1  J<--""/o(v)  (1  +  sin2s  vr  )  </fl^>.  (A4) 

where  <  >  represents  an  ensemble  average  over  both  group  delay 
x'  and  photon  statistics. 

A  better  estimator  is  the  unbiased,  normalized  power  spectrum 
Q(t),  given  by 


G(t)  = 


</»(T)-f/(v)</v> 

<//(v)dv>" 


(A  5) 


where  N,=  \l(v)dv  is  the  total  photon  flux  in  a  tOms  interval. 
Comparison  of  this  expression  with  the  estimator  (11)  shows  that 
(A  5)  is  just  the  estimator  evaluated  at  a  particular  group  delay. 
Thus,  the  performance  of  the  group  delay  estimator  can  be 
calculated  using  (A  1 ). 

Given  a  spectrometer  with  a  bandwidth  of  0.4  pm  to  0.7  pm. 
the  delay  resolution  is  dx  =  (l/0.4-  l/0.7)‘‘  ~  1  pm.  Since  the 
group  delay  fluctuates  with  a  standard  deviation  of  a,,  the 
gaussian  distribution  of  group  delay  in  Q(x)  will  be  spread  over 
approximately  aJAx  bins,  while  the  amplitude  of  the  gaussian  will 
be  reduced  by  approximately  aJAx.  Photon-starved  statistics 
apply  for  (1/2)  N,  F’  <  <r,/dx.  and  using  (A  1 ).  the  signal-to-noise 


ratio  for  an  estimate  of  the  amplitude  of  the  gaussian  using  a 
matched  filter  is  given  by 

SIV/t((?(T))=  (A6) 

I  .j  V 

where  N,  is  the  total  photon  flux  in  the  spectrum.  In  (A6)  sampling 
losses  are  neglected,  and  it  is  assumed  that  a,  p  Ax.  Using  (9)  and 
(A  3)  to  write  (A6)  in  terms  of  N,  rather  than  N, ,  yields  a  signal-to- 
noise  ratio  for  group  delay  estimates  which  is  nearly  identical  to 
the  signal-to-noise  ratio  of  a  conventional  amplitude  estimate. 

For  the  case  of  group-delay,  an  estimate  of  the  centroid  of  the 
group-delay  distribution  is  desired.  The  standard  deviation  of 
such  an  estimate  can  be  written 


(A7) 


where  the  S.V/?  is  given  by  (A6).  A  Monte  Carlo  simulation  was 
performed  which  verified  the  above  analytically  derived 
expressions. 

Under  photon-starved  conditions,  the  measurement  of  closure 
phase  is  qualitatively  dilTerent  from  the  measurement  of  ampli¬ 
tude.  The  definition  of  the  triple  product.  Eq.  (14).  is  repeated 
below: 


triple  product  =  (4,  c'*’'  /Ijc''’  4jc"''>.  (A8) 

In  photon-counting  measurements,  if  zero  photons  are  detected  in 
any  one  of  the  three  interferometers,  the  triple  product  is  zero.  For 
a  frame  to  contain  some  information  on  the  closure  phase,  at  least 
one  photon  must  be  detected  in  all  three  interferometers.  The 
Poisson  probability  function  is  given  by 

A 

/>(«!)=  —  .  (A9) 

nr. 

where  P{m)  is  the  probability  of  detecting  m  photons  given  an 
average  arrival  rate  of  N.  For  A  1  the  probability  of  detecting 
one  photon  is  A.  Thus,  the  probability  of  detecting  one  photon  in 
all  three  interferometers  is  A^.  so  that  the  variance  for  estimates  of 
closure  phase  should  be  proportional  to  A"'*  under  photon 
starved  conditions. 

More  precisely,  the  variance  of  the  closure  phase  calculated  via 
the  triple  product  can  be  written 

"i- -’(a.;-) 

for  a  signal-to-noise  ratio  SNR  =Mri^.  In  this  expression  A  is  the 
average  number  of  detected  photons  in  each  interferometer.  Equal 
visibilities  T  have  been  assumed  for  each  of  the  three  inter¬ 
ferometers.  Thus,  at  high  photon  fluxes,  the  variance  is  three  times 
the  variance  in  a  photon-rich  phase  measurement.  At  low  photon 
fluxes,  the  variance  is  proportional  to  A"  ’,  as  given  by  the  last 
term  of  (A  10).  For  this  term.  I'  is  the  geometric  mean  of  the 
visibilities  in  the  three  interferometers.  A  Monte  Carlo  simulation 
for  closure  phase  measurements  was  also  conducted,  which 
verified  the  above  expression. 

The  integration  time  needed  to  achieve  a  signal-to-noise  ratio 
of  30  versus  photon  flux  A  is  given  in  Fig.  II  for  measurements  of 
visibility  squared  and  group  delay,  and  closure  phase.  For 
astronomical  measurements,  a  photon  flux  corresponding  to  an 
SNR  of  30  on  an  unresolved  star  allows  a  diameter  measurement 
of  4%  on  a  partially-resolved  source.  An  SNR  of  30  for  group 
delay  corresponds  to  an  astrometric  measurement  of  0.6  pm  in 


FIf.  II.  Integration  time  vs  magnitude  for  photon-starved  measurements 


delay,  or  Oroi  in  angle.  An  SNF  of  30  for  the  triple  product 
corresponds  to  an  error  of  /./1 80  in  closure  phase.  The  limiting 
magnitudes  for  a  single-rg  system  are  approximately  14.3  and  12.5 
for  visibility/group  delay  and  closure  phase.  With  l-m  aterpures. 
the  limiting  magnitudes  are  approximately  16.8  and  14.2. 
respectively. 

If  rg  and  Xq  are  different  from  the  nominal  values  assumed,  the 
magnitude  limits  will  change.  In  general,  an  instrument  will  have  a 
fixed  collecting  aperture.  Thus,  if  rg  increases,  the  number  of 
photons.  N.  per  Xg  increases,  but  the  number  of  frames.  M. 
decreases.  Similarly,  given  a  fixed  10'*  s  for  a  measurement,  an 
increase  in  Xg  will  also  increase  y  and  decrease  M.  If  the  signal-to- 
noise  ratio  were  dependent  on  /V  and  M  evenly,  the  final  SJVR 
would  be  independent  of  seeing.  This  is  the  case  for  photon-rich 
measurements  where  the  SJV/l  is  proportional  to  y^NM . 

For  photon-starved  measurements,  this  is  not  the  case.  The 
signal-to-noise  ratio  for  amplitude  and  delay  measurements 
follow  the  scaling  relation  NM'  while  the  SNR  (ox  closure  phase 
follows  the  relation  N^'^  M'  Two  additional  facts  should  be 
considered.  As  rg  improves,  the  rms  path  length  fluctuation 
between  the  apertures  decreases  proportionally.  Hence,  the  num¬ 
ber  of  required  spectral  channels  that  must  be  separately  analyzed 
decreases  according  to  (9).  In  addition,  since  the  temporal 
fluctuations  of  phase  are  usually  modeled  by  assuming  frozen 
turbulence  translating  past  the  instrument  at  a  constant  wind 
speed.  Xg  should  increase  with  rg .  Thus,  the  scaling  of  N  and  M 
with  Xg  is  approximately  given  by 

Nxxt.  (AM) 

Mxxg*.  (A  12) 


The  scaling  of  the  magnitude  limits  for  photon-starved  measure¬ 
ments  with  seeing  is  then  given  by 

Amplitude/delay  sensittivity  oc  rg .  (A  1 3) 

Phase  closure  sensitivity  oc  rj'-’ .  (A  14) 

These  results  should  be  compared  to  the  case  of  a  wideband  phase 
tracking  system.  In  this  case,  only  N.  not  M  is  important,  and  since 
the  fringe  is  tracked,  there  is  no  need  for  the  spectrometer.  Thus, 
the  sensitivity  scales  with  seeing  as  Xg . 
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ABSTRACT 

We  report  initial  stellar  diameter  measurements  made  with  the  Mark  III  stellar  interferometer.  The  visibility 
measurements  showed  fluctuations  of  1  %-j%  for  a  10  s  integration  consistent  with  photon  sutistics  as  the 
predominant  source  of  noise.  Measurements  of  the  unresolved  star  ^  Tau  were  used  to  calibrate  the  instru¬ 
mental  and  atmospheric  biases  which  reduce  the  measured  visibility.  We  report  angular  diameters,  at  706  nm, 
of  6.15  ±  0.25  mas,  8.70  ±0.15  mas,  5.0  ±  0.5  mas,  and  7.2  ±  0.4  mas  rms  for  the  stars  «  CMi,  Gem,  a  Per, 
and  a  UMa,  respectively.  These  angular  diameters  were  calculated,  assuming  a  uniform  disk,  from  the  visibil¬ 
ity  measurements,  the  wavelength  of  observation,  and  the  projected  interferometer  baseline,  llte  angular  diam¬ 
eters  corrected  for  limb  darkening  are  also  presented.  In  the  next  observing  season,  the  instrument  should  be 
able  to  make  very  precise  diameter  measurements  on  a  large  number  of  stars  brighter  than  7  mag. 

Subject  headings:  instruments  —  interferometry  —  stars;  diameters 


1.  INTROOUCnON 

The  Mark  III  interferometer  is  a  joint  project  among  the 
Smithsonian  Astrophyskal  Observatory,  the  Massachusetts 
Institute  of  Technology,  the  Naval  Research  Laboratory,  and 
the  U.S.  Naval  Observatory.  The  instrument  is  located  on 
Mount  Wilson,  California.  First  fringes  were  observed  in  1986 
September  with  a  12  m  baseline. 

Long-baseline  optical  stellar  interferometers  have  tradi¬ 
tionally  been  used  to  measure  the  diameter  of  stars.  The  Mark 
III  interferometer  (Shao  et  at.  1988),  unlike  other  operational 
stellar  interferometers  (Labeyrie  1975;  Davis  and  Tango  1986), 
was  designed  principally  for  fundamental  astrometry 
(Mozurkewich  et  at.  1987).  Consequently,  the  instrument  was 
designed  with  such  characteristics  as  extreme  thermal  stability 
and  submicron  mechanical  accuracy  and  makes  extensive  use 
of  laser  metrology.  These  features  are  not  usually  necessary  for 
stellar  diameter  measurements.  In  addition,  the  instrument  was 
optimiad  for  high  astrometric  sensitivity  through  the  use  of 
very  wide  optical  bandpatnes,  at  the  expense  of  some  accuracy 
in  the  measurement  of  visibilities.  It  was  subsequently  modified 
to  enable  accurate  visibility  measurements  to  be  made  simulta¬ 
neously  with  astrometric  measurements,  and  a  set  of  prelimi¬ 
nary  stellar  diameter  measurements  were  made  in  1986 
November.  The  next  section  describes  the  hardware  used  for 
making  these  initial  visibility  measurements,  and  the  third 
section  describes  the  observation  and  analysis  procedure.  The 
last  section  summarizes  these  initial  measurements  of  a  CMi, 
p  Gem,  a  Per,  and  a  UMa. 

n.  the  mark  m  interferometer 

The  instrument  is  an  active-fringe-tracking  long-baseline 
optical  interferometer.  Figure  I  shows  a  schematic  of  the 


instrument.  A  star  tracker,  or  autoguider,  keeps  the  interfering 
wavefronts  parallel  to  a  fraction  of  an  arcsec.  The  vacuum 
delay  tines  are  monitored  with  laser  interferometers  with  5  nm 
resolution  and  can  be  set  with  10  nm  accuracy.  An  active 
fringe-phase-tracking  servo-system  keeps  the  optical  paths  in 
the  two  arms  equal  to  a  fraction  of  a  wavelength  of  light  in  the 
presence  of  mechanical  error,  thermal  drift,  and  atmospheric 
turbulence.  Presently,  light  from  one  side  of  the  beamsplitter  is 
used  for  fringe  tracking  and  light  from  the  other  side  is  used  for 
visibility  measurements  through  a  narrow  bandpass  filter.  In 
1986,  the  magnitude  limit  of  the  instrument  was  set  by  an 
int-rim  star  tracker  rather  than  by  the  fringe  tracker.  In  1987, 
an  improved  'tar-tracking  system  (Shao  et  at.  1988)  will  enable 
visibility  measurements  of  F  =  7  objects. 

The  fringe-detection  scheme  is  based  on  a  path-length 
modulation  technique  that  was  also  used  on  two  previous 
long-baseline  stellar  interferometers  on  Mount  Wilson.  A 
mirror  on  a  piezoelectric  stack  (PZT)  is  moved  in  a  500  Hz 
triangle-wave  pattern  with  a  one-wavelength  p-p  motion. 
During  each  I  ms  sweep  of  the  PZT,  the  number  of  detected 
photons  in  four  equally  spaced  time  bins  is  recorded,  as  shown 
in  Figure  2.  Successive  I  ms  scans  are  co-added  to  yield  the 
desired  coherent  integration  time.  These  quantities — A,  B,  C. 
and  D — are  used  to  estimate  the  visibility.  The  amplitude  of  the 
path-length  modulation  is  monitored  by  a  laser  interferometer 
and  is  controlled  by  a  low-bandwidth  servo  with  a  precision  of 
^0.1%  to  reduce  the  systematic  errors  that  arise  when  the 
modulation  amplitude  is  different  from  the  wavelength  of 
observation.  In  order  to  minimize  the  systematic  erron  in  the 
visibility  measurements  attributable  to  atmosfflwric  turbu¬ 
lence,  the  amplitude  measurements  were  made  with  an  effective 
aperture  of  2.5  cm,  an  optical  bandwidth  of  10  nm,  and  a 
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Fig.  1. — Schematic  diagram  of  the  Mark  III  stellar  interferometer 


coherent  integration  time  of  4  ms.  For  a  complete  description 
of  the  iiutniment,  see  Shao  et  ai.  (1988). 

The  detectors  used  for  the  fringe  measurements  are  GaAs 
PMTs  cooled  to  —  SO^C.  The  dark  count  at  this  temperature  is 
approximately  130  hz.  The  PMTs  are  fed  with  200  pm  optical 
fibers  and  are  sensitive  to  sky  background  in  a  l6Cr  diameter 
field.  Thus,  to  make  accurate  measurements  of  visibility,  the 


detector  dark  count  and  the  sky  background  must  be  accu¬ 
rately  measured.  The  measurement  of  these  quantities  was 
made  part  of  the  automated  observing  sequence. 

in.  observations  and  analysis 

The  measurements  presented  here  were  made  on  1986 
November  21.  The  stars  to  be  observed,  all  FK4  stars,  were 
programmed  for  observation  at  the  beginning  of  the  night.  The 
control  computers  then  directed  the  interferometer  to  make 
observations  in  the  preprogrammed  sequence.  The  program 
specified  which  stars  were  to  be  observed,  in  what  order,  and 
how  much  time  was  to  be  spent  on  each  star.  Data  were 
recorded  whenever  the  fringe  tracker  was  locked  on  a  fringe. 
These  data,  recorded  at  a  4  ms  rate,  included  the  time,  the 
delay-line  position,  and  the  quantities  A,  B,  C,  and  D  for  the 
amplitude  and  astrometric  channels.  After  each  stellar  obser¬ 
vation.  while  the  siderostats  were  slewing  to  the  next  star,  a  S  s 
measurement  of  oark  count  and  sky  background  was  also 
recorded.  Amplitude  measurements  were  made  using  a 
narrow-bandpass  filter  at  706  nm;  some  preliminary  measure¬ 
ments  at  604  nm  were  corrupted  by  scattered  laser  light  and 
are  not  included  here. 

Using  the  bin  counts  A,  B,  C,  and  D,  the  square  of  the  fringe 
visibility  is  calculated  as 

<(/t  -  -I-  (B_-  D)^  -  \> 

2  <^  -  Dry 

where  AI  =  A  -i-  B  -i-  C  -i-  D,  Dr  is  the  average  dark  count  per 
4  ms,  and  angle  brackets  denote  time  average.  This  estimator 
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is  unbiased  with  respect  to  photon  noise  and  dark  count  In 
addition,  it  can  be  shown  that  this  estimator  is  not  biased  by 
scintillation  on  time  scales  slower  than  the  coherent  integra¬ 
tion  time  (Tango  and  Twiss  1980). 

Atmospheric  turbulence  produces  systematic  reductions  in 
the  measured  visibility.  The  use  of  Bnite-diaweter  collecting 
apertures  causes  a  reduction  in  visibility  attributable  to  the 
rdative  tilts  and  curvatures  of  the  interfering  wavefronts.  The 
use  of  a  finite-length  coherent  integration  time  causes  a 
reduction  in  visibility  attributable  to  motion  of  the  fringe 
during  the  interval.  The  reduction  in  visibility  due  to  these 
effects  can  be  approximated  as  (cf.  Tango  and  Twiss  1980) 

*  Kn..  exp  1^ -  ‘  03(0’  ']  <2) 

and 

*=  Kru.  exp  [ -0-5(0’"]  <3) 

for  D/fo  and  t/to  <  1.  In  these  expressions,  D  and  r  are  the 
aperture  diameter  and  coherent  integration  time;  ro  is  the 
coherence  diameter,  i.e.,  that  diameter  aperture  over  which  the 
rms  wavefront  fluctuations  are  approximately  equal  to  1 
radian;  and  to  is  the  coherence  time,  i.e.,  that  time  interval 
during  which  the  rms  fluctuations  of  the  interference  fringe  are 
equal  to  1  radian.  Equation  (2)  assumes  the  use  of  a  slow  angle 
tracker — one  which  removes  only  systematic  pointing  errors. 
With  a  fast  angle  tracker,  which  removes  some  of  the  wave- 
front  tilt,  the  reduction  in  visibility  can  be  made  smaller.  For  I" 
seeing  with  a  10  m  s~‘  wind,  «  14  cm  and  Tq  %  12  ms  at 
700  nm;  with  O  ^  2.5  cm  and  r  ==  4  ms,  the  visibility  reduction 
is  approximately  13.4%. 

An  additional  source  of  error  is  finite  temporal  coherence. 
For  example,  with  a  bandpass  filter  which  is  40  nm  wide  and  a 
one-fringe  servo  tracking  error,  the  reduction  in  fringe  visibility 
is  0.3%.  With  the  10  nm  wide  filter  used  for  these  measure¬ 
ments,  the  visibility  reduction  due  to  partial  temporal  coher¬ 
ence  is  very  small.  This  large  usable  optical  bandpass  is  a 
significant  advantage  of  stellar  interferometers  that  track  the 
central  fringe. 

A  final  source  of  bias  is  the  instrument  itself  An  asym¬ 
metrical  beamsplitter  and  the  wavefront  distortion  caused  by 
imperfect  optics  will  both  lower  the  visibility.  Our  calibration 
procedure  for  systematic  errors  due  the  atmosphere  and  the 
instrument  consists  of  observing  a  star  whose  intrinsic  visibility 
is  over  95%.  Thus,  as  part  of  the  observing  sequence,  unre¬ 
solved  stars  are  periodically  measured. 

The  data  analysis  procedure  begins  with  a  data  qualification 
step.  Depending  on  the  seeing,  the  fringe  tracker  will  remain 
locked  for  a  period  of  a  few  tenths  of  a  second  to  several 
minutes.  Data  are  recorded  only  when  the  fringe  tracker  is 
locked.  In  the  data  analysis,  a  short  segment  of  data  at  the 
beginning  and  end  of  each  lock  representing  the  acquisition 
and  loss  of  the  fringe  is  discarded.  In  addition,  if  the  length  of 
the  contiguous  track  of  data  is  less  than  0.5  s,  the  data  is  also 
discarded.  From  this  qualified  data,  is  calculated  and  aver¬ 
aged  to  10  s. 

After  the  quantity  V*  is  averaged  to  10  s,  the  next  step  is  to 
correct  for  the  biases  described  earlier.  This  is  done  by  normal¬ 
izing  the  measured  visibilities  by  the  visibility  of  an  unresolved 
star.  P  Tau  was  used  as  the  reference  star  for  these  observa¬ 
tions.  Its  diameter  can  be  estimated  from  its  spectral  type. 


luminosity  class,  and  magnitude  as  1.3  mas,  which  yields  an 
assumed  visibility  of  98.6%  for  a  projected  baseline  of  12  m. 
(Comparable  diameter  estimates  of  1.1  mas  and  1.6  mas  were 
found  by  Wesselink,  Paranya,  and  De  Vorkin  1972;  and 
Hertzsprung  1922,  respectively.)  From  measurements  of  p  Tau, 
the  normalization  constant  for  these  measurements  was  deter¬ 
mined  as  0.751,  i.e.,  a  visibility  reduction  of  24.9%.  x  Leo, 
which  has  a  measured  diameter  (Hanbury  Brown,  Davis,  and 
Allen  1974)  was  not  used  as  a  calibrator  because  it  was  at  a 
much  larger  zenith  distance  and  was  observed  only  for  a  short 
period  of  time. 

Figure  3  shows  the  normalized  visibilities  for  the  seven  stars 
observed.  The  gaps  in  the  plot  are  caused  by  temporarily  poor 
seeing  during  which  the  fringe  tracker  would  not  stay  locked 
on  the  white-light  fringe  for  more  than  0.5  s;  these  data  are 
rejected  by  the  data  selection  algorithm.  Estimates  of  the 
coherence  parameters  for  this  data  set  ate  tq  «  20  cm  and 
To%11.4ms  (cf.  Cola  vita,  Shao,  and  Staelin  1987),  corre¬ 
sponding  to  slightly  better  than  1”  seeing.  With  these  values, 
the  visibility  reduction  due  to  the  atmosphere  is  approximately 
11.3%.  Thus,  given  the  normalization  constant  determined 
above,  the  reduction  in  visibility  attributable  to  the  instrument 
is  approximately  15.4%. 

In  order  to  evaluate  the  stability  and  accuracy  of  our  mea¬ 
surements,  we  measured  the  variance  of  the  raw  measure¬ 
ments,  integrated  for  10  s,  over  a  period  of  3  hr.  This  was 
compared  to  the  theoretical  variance  due  to  photon  noise.  The 
photon  noise  variance  for  the  F^  estimator  of  equation  H )  is 
given  by 


a 


2 


5!  J-_L 

4  M 


(4) 


where  F  is  the  raw  visibility,  N  is  the  number  of  photons  per 
frame  (4  ms),  and  M  is  the  number  of  frames.  Table  1  shows  the 
measured  and  theoretical  variance  of  F*  for  the  three  stars 
P  Tau,  P  Gem,  and  x  CMi.  There  were  insufficient  data  on  the 
other  stars  for  a  meaningful  variance  measurement. 

Table  I  shows  that  the  fluctuations  of  the  raw  F^  measure¬ 
ments  are  consistent  with  photon  noise,  which  is  at  the 
l%-3%  level  for  a  10  s  integration.  As  discussed  earlier,  the 
measured  visibility  is  dependent  upon  the  atmospheric  seeing 
Thus,  if  the  seeing  changes  during  a  series  of  observations,  i.e., 
if  rg  and  Tq  do  not  remain  constant,  the  variance  of  the  mea¬ 
sured  F^  wil!  be  larger  than  the  theoretical  photon  noise 
variance  given  by  equation  (4).  The  good  agreement  of  the 
measured  and  theoretical  values,  which  correspond  to  stan¬ 
dard  deviations  of  only  1.2%  for  P  Gem  and  x  CMi.  indicates 
that  atmospheric  turbulence  should  not  be  a  significant  source 
of  error  in  this  set  of  measurements. 

In  general,  the  sensitivity  of  the  measured  visibility  to 
changes  in  tq  and  Tq  is  easily  determined  from  equations  (2) 
and  (3): 


Kneisured  \^oJ  \  ^0  / 


(5i 


TABLE  I 

Measlxed  asu  THEORE-nCAL  Variance  of  t'^ 


Star 

s 

I  -  Irawl 

„2 

^  mtMured 

Tau  . 

4.22 

2500 

054* 

10.6  X  10- * 

8.4  X  10“ 

P  Gem  . 

13.45 

2500 

0.137 

0.94  X  10’* 

1.33  X  10-' 

jCMi  . 

16.26 

2500 

0,337 

1.31  X  10-* 

140  X  10“ 
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and 


^mtmr*4  \^oJ  \  ^0  , 


(6) 


Using  the  estimated  parameters  for  this  data  set,  10%  changes 
in  ro  and  Tq  lead  to  changes  in  the  measured  visibility  of  O.S% 
and  2.0%,  respectively.  In  addition,  the  precise  values  of  the 
coherence  parameters  depend  on  zenith  angle.  The  coherence 
diameter  rg  is  inversely  proportional  to  the  ^  power  of  the 
turbulence  strength,  and  thus  has  a  zenith  angle  dependence 
Tq  oc  oos*^*  r.  The  coherence  time  Tq  is  approximately  pro¬ 
portional  to  r^/W,  where  W  is  the  wind  spewl  normal  to  the 
propagation  path.  Averaging  W  with  respect  to  wind  direction 
yields  for  to  an  approximate  zenith  angle  dependence  to  x 
2  CM*'*  z/(l  +  cos  z).  Using  the  estimated  parameters  for  this 
data  set,  the  systematic  reduction  in  measured  visibility  at  a 
30°  zenith  angle  is  approximately  0.7%.  These  systematic  visi¬ 
bility  reductions  have  been  factored  into  the  normalized  visi¬ 
bilities. 

Beginning  with  our  next  observing  season,  when  we  begin 
observations  of  fainter  stars,  apertures  larger  than  2.S  cm  and 
coherent  integration  times  longer  than  4  ms  will  be  used.  Con¬ 
sequently,  turbulence-induced  errors  will  increase,  and  careful 
calibration  of  atmospheric  effects  will  be  crucial.  One  analysis 
procedure  that  could  be  used  to  reduce  these  errors  is  to 
analyze  the  data  with  two  different  coherent  integration  times 
and  two  different  coherent  aperture  sizes.  With  these  measure¬ 
ments  it  will  then  be  possible,  in  theory,  to  solve  for  both  the 
true  visilHlity  as  well  as  ro  and  Tq. 

Stellar  diameters  are  estimated  from  the  normalized  visibil¬ 
ity,  the  effective  wavelength  of  observation,  and  the  projected 
baseline.  The  filter  used  for  these  observations  was  centered  at 
706  nin,  and  the  projected  baseline,  routinely  determined  to 
one  part  in  10*  as  part  of  the  astrometric  data  reduction,  was 
nominally  12  m  for  observations  at  zenith.  Assuming  a 
aaform  stellar  disk,  the  diameters  of  several  of  the  stars  whose 


normalized  visibilities  are  shown  in  Figure  3  are  presented  in 
column  (3)  of  Table  2.  The  uniform  disk  angular  diameters 
were  corrected  for  limb  darkening  based  upon  the  appropriate 
model  atmospheres  from  Carbon  and  Gingerich  (1969)  for 
*  CMi.jJ  Gem,andx  UMa,and  Parsons (1971) for s  Per. 

A  cosine  law  was  adopted  and  appropriate  values  for  the 
limb-darkening  coefficient  were  deriv^  from  the  model  atmo¬ 
spheres.  The  angular  diameters  were  then  corrected  using  the 
relationship  between  the  limb-darkened  diameter,  the  uniform 
disk  diameter,  and  the  limb-darkening  coefficient  given  by 
Hanbury  Brown  and  Twiss  (1958).  The  estimated  “true "diam¬ 
eters  of  the  stars  thus  derived  (Table  2,  col.  [4])  are  5%-7% 
larger  than  the  equivalent  uniform  disks. 

The  errors  given  for  the  uniform  disk  diameters  (col.  [3])  are 
the  formal  rms  errors  in  the  measurement  and  in  the  normal¬ 
ization  process.  These  errors  include  the  effects  of  photon 
noise,  as  well  as  the  uncertainties  in  the  assumed  angular  diam¬ 
eter  of  the  calibration  star  P  Tau.  As  noted  before,  a  true  visi¬ 
bility  of  98.6%  was  assumed  for  P  Tau.  We  believe  our 
estimate  for  the  diameter  of  P  Tau  is  good  to  ±30%.  A  30% 
error  in  the  estimate  of  its  diameter  would  produce  a  1% 


TABLE  2 

Sttllar  DiAMrTWS  at  706  NANOiHFrens 


Total  Asoula*  Diameter 

Observed  - 

Time  Uniform  Disk  Limb  Dtrkened  Disk 


Star  (s|  (mas)  ±  <r  (rms)  (mas)  ±  a  (rms) 

(I)  (2)  (3)  (4) 


iCMt  .  80  6.15  ±  0.25  6.44  ±  0.25 

^Gem  .  100  8.70  +  0.15  9.26  ±0.15 

» Per  .  20  5.0(  ±  O  S)'  5J(  ±  O.Sf 

I  UMa  .  10  7.2(  ±  0.4)‘  7.6(  ±  0.4r 

%  Gem .  10  (multiple  star) 

X  Leo .  20  (unresolved)  >95%  visibility 


*  Parentheses  indicate  theoretical  photon  noise  error  only. 
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change  in  the  visibility  used  to  nonnalize  the  other  measure¬ 
ments.  In  addition,  the  errors  include  the  effects  of  lOVo  rms 
fluctuations  in  the  values  of  and  to  between  observations  of 
the  calibrator  of  fi  Tau  and  the  other  stars.  The  errors  in  the 
limb-darkening  corrections  are  estimated  to  be  of  the  order  of 
±0.5%  due  to  deviations  in  the  (a|q>roximate)  cosine  relation 
used  from  the  model  atmospheres,  litis  latter  effect  is  included 
in  the  errors  listed  in  column  (4)  of  Table  2.  While  we  believe  all 
significant  sources  of  error  have  been  induded,  additional  data 
with  repeated  observations  at  different  baselines  and  at  differ¬ 
ent  wavdengths  will  be  needed  to  confirm  the  accura^  of  these 
diameter  measurements. 

IV.  DISCUSSION 

A  limb-darkened  angular  diameter  of  S.SO  ±  0.17  mas  was 
determined  for  a  CMi  by  Hanbury  Brown,  Davis,  and  Allen 
(1974)  from  measurements  at  440  nm  made  with  the  Intensity 
Interferometer.  The  limb-darkened  angular  diameters  of 
fi  Gem  and  a  Per  found  in  the  literature  are  7.8  ±  0.6  mas  (at 
SSO  nm)  by  Faucherre  et  al.  (1983)  for  fi  Gem,  and 
3.29  ±  0.2S  mas  (at  600  nm),  and  2.9  ±  0.4  mas  (at  SSO  nm)  for 
a  Per  by  Koechlin  and  Rabbia  (1985)  and  Bonneau  et  al. 
(1981),  respectively.  Our  measurement  for  a  UMa  shows  rea¬ 
sonable  agreement  with  an  indirect  estimate  of  the  limb- 
darkened  diameter  based  on  absolute  infrared  photometry  of 
6.64  ±  0.21  mas  by  Blackwell  and  Shallis  (1977). 

Our  diameter  measurements  are  systematically  larger  than 
previous  interferometric  measurements.  Consequently,  we 
have  looked  for  possible  systematic  biases  to  our  measure¬ 
ments.  a  Leo  was  measured  by  Hanbury  Brown  to  be  1.2  mas 
and  should  have  given  a  visibility  of  98.5%.  Our  measurement 
of  a  Leo  gave  a  result  of  96%  with  2%  rms  photon  noise.  The 
discrepancy  is  1.2  <r.  a  Leo  was  not  used  as  a  calibrator 
because  its  photon  statistics  were  poor.  However,  it  does 
provide  a  check  on  our  diameter  measurements  which  are 
otherwise  systematically  different  from  previous  measure¬ 
ments. 

There  are  many  possible  sources  of  systematic  error  for  our 
current  and  previous  measurements.  The  seeing,  r,,,  xg  can 
change  with  the  seasons,  the  time  of  night,  the  average  wind 
speed,  and  there  may  be  effects  that  correlate  with  the  magni¬ 
tude  and  color  of  the  star  since  the  star  tracker  and  fringe 
tracker  work  better  on  brighter  stars  and  use  detectors  most 
sensitive  to  red  light.  The  most  important  of  these  errors  are 
included  in  our  estimated  errors  in  Table  2.  We  cannot  explain 
the  systematic  difference  between  our  recent  measurements 
and  previous  measurements.  However,  these  measurements 
were  made  at  a  longer  wavelength  than  previous  diameter 
measurements.  In  the  past,  observation  of  stellar  fringes  was  a 
major  accomplishment  in  itself.  In  the  future,  we  expect  the 


emphasis  will  be  on  systematic  errors  such  as  the  ones  dis¬ 
cussed  in  this  paper  and  differences  in  measured  diameters  at 
different  times  would  have  astrophysical  significance. 

Future  operation  of  the  instrument  will  be  different  in  three 
ways.  Most  important,  with  the  installation  of  a  new  photon¬ 
counting  camera  as  a  part  of  the  star-tracking  system  (1987 
July),  active  fringe  tracking  and  hence  accurate  visibility  mea¬ 
surements  will  be  possible  for  stars  to  at  least  V  =■  5.0 
(probably  to  F  =  7).  These  improvements  will  allow  a  dra¬ 
matic  increase  in  the  number  of  stars  available  for  study.  We 
plan  to  conduct  a  long-term  program  to  repeatedly  observe  the 
angular  diameters  of  approximately  200  stars  at  many  wave¬ 
lengths. 

In  addition  to  providing  direct  measurements  of  photo- 
spheric  angular  diameters  and  limb  darkening,  we  intend  to 
observe  a  number  of  stars  which  are  likely  to  possess  extended 
atmospheres  or  shells  through  narrow-band  filters  or  polari¬ 
zers  for  evidence  of  resolved  structures.  Extended  emission  in 
Ha  from  an  envelope  surrounding  the  Be  star  y  Cas  has  been 
reported  by  Thom,  Granes,  and  Vakili  (1986)  from  observa¬ 
tions  made  with  the  CERGA  I2T  interferometer.  We  plan  to 
observe  this  source  as  well  as  other  Be  stars  with  known  IR 
excesses  (Gehrz,  Hackwell,  and  Jones  1974);  supergiants  such 
as  p  Cas,  HR  8752  (Lambert,  Hinkle,  and  Hall  1981),  and 
P  Cyg:  and  the  Mira  variables  R  Leo  and  y  Cyg.  With  the 
extension  of  the  maximum  baseline  length  to  20  m  in  1987- 
1988,  diameter  measurements  of  the  nearer  Cepheid  variables 
will  be  attempted. 

Second,  we  plan  to  install  a  low-resolution  spectrometer  and 
photon-counting  camera  detector  so  that  fringe  visibility  at  all 
wavelengths  can  be  measured  simultaneously.  Third,  we  plan 
to  add  software  for  real-time  baseline  calibration.  This  capabil¬ 
ity  will  allow  operation  as  an  absolute  interferometer  in  order 
to  make  diameter  measurements  of  objects  up  to  30  times 
fainter  than  the  7  mag  limit  of  the  fringe  tracker  using  5  cm 
apertures. 

V.  CONCLUSIONS 

Stellar  diameters  of  a  CMi,  p  Gem,  a  Per,  and  a  UMa  mea¬ 
sured  with  the  Mark  Ill  stellar  interferometer  are  reported. 
Values  are  given  for  the  cases  of  uniform  and  limb-darkened 
disks.  The  measurements  were  performed  using  a  totally 
automated  instrument  which  will  be  capable  of  making  mea¬ 
surements  of  stars  brighter  than  V  =  7  with  diameters  from 
07003  to  07015  and  located  between  5°  and  65°  declination. 


This  work  was  made  possible  in  part  by  NRL  contract 
N000l4-86-C-2t05  and  the  Smithsonian  Scholarly  Studies 
Program. 
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Atmospheric  phase  measurements  with  the  Mark  III 
stellar  Interferometer 

M.  Mark  Cotavita,  Michael  Shao,  and  David  H.  Staeiin 


The  Mark  III  interferometer  is  a  phase-coherent  stellar  interferometer  designed  for  astrometry.  Operating 
through  the  turbulent  atmosphere,  the  instrument  is  also  a  sensitive  detector  of  atmospheric  phase  fluctua¬ 
tions.  The  effect  of  phase  fluctuations  on  astrometric  accuracy  is  reviewed,  and  phase  measurements 
obtained  with  the  instrument  at  Mt.  Wilson  using  a  12-m  base  line  are  presented.  These  measurements  agree 
well  with  the  predictions  of  a  simple  Kolmogorov  spatial  spectrum  over  the  frequency  range  of  0.001-100  Hz. 
From  these  measurements,  the  outer  scale  of  turbulence  for  propagation  through  the  entire  atmosphere  is 
estimated  to  be  >2  km.  The  standard  deviation  for  an  absolute  astrometric  measurement  estimated  from 
these  measurements  is  ~0.147^'^  sec  of  arc  for  long  integration  times  for  conditions  of  0.5-sec  of  arc  seeing. 
For  star-switched  relative  measurements,  this  error  should  decrease  as  the  square  root  of  the  number  of 
switching  cycles. 


i«  vivoQucni 

The  Mark  III  aatrometric  interferometer  is  a  phaae- 
ooherent  white-light  stellar  interferometer  designed 
for  astrometry.  It  became  operational  on  Mt.  Wilson 
in  Sept  1986  using,  at  present,  a  12-m  N-S  base  line. 
It  succeeds  the  3.1-m  Mark  II  interferometer*  which 
operated  on  Mt  Wilson  from  Aug.  1982  until  Sept. 
1984  and  incorporates  a  number  of  improvements  over 
that  instrument,  especially  in  the  area  of  thermal  sta¬ 
bility.  The  instrument  measures  the  phase  of  the 
white-light  fringe  in  two  wide  spectral  channels  cover¬ 
ing  ~0.6-0.9  and  0.4-0.6  tan.  The  fringe  phase  in  the 
red  channel  serves  as  the  error  signal  for  a  white-light 
fringe  servo  which  controls  the  position  of  a  IsMr- 
monitorsd  delay  line,  while  the  fringe  phase  in  the  blue 
channel  is  used  with  the  two-color  method.  Continu¬ 
ous  fringe  tracks  vary  from  a  fraction  of  a  second  to 
tens  of  seconds,  depending  on  the  seeing.  Figures  1 
and  2  present  a  photograph  of  the  Mark  III  interferom¬ 
eter  site  and  a  schematic  diagram  of  the  instrument. 
The  instrument  is  discussed  in  detail  in  a  separate 
paper.* 
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Atmospheric  turbulence  is  the  fimdamental  factor 
limiting  ^e  performance  of  all  ground-based  astrome¬ 
tric  instruments,  affecting  both  sensitivity  and  accura¬ 
cy.  While  in  practice  instrumental  limitations  can 
dominate  performance,  the  intent  in  the  design  of  the 
Mark  III  instrument  was  to  keep  systematic  errors 
sufficiently  small  that  atmospheric  turbulence  would 
be  the  only  significant  error  mechanism  when  rapid 
star  switching  was  employed.  Section  11  describes  the 
effects  of  turbulence  on  measurements  of  stellar  posi¬ 
tion,  both  for  the  case  of  a  simple  Kolmogorov  spatial 
spectrum,  and  for  the  case  of  a  model  incorporating  a 
finite  outer  scale.  Section  III  presents  atmospheric 
phase  measurements  obtained  with  the  Mark  EH  inter¬ 
ferometer  at  Mt  Wilson  in  the  fall  of  1986.  Data  from 
the  Mark  II  instrument  are  also  presented.  Only  one- 
color  observations  are  consider^  in  this  paper.  The 
use  of  two-color  methods  to  reduce  atmospheric  errors 
is  an  important  part  of  the  operation  of  the  Mark  III 
instrument  and  is  discussed  in  detail  in  a  separate 
paper.* 

AUiiOSpffVfsC  BViVClS 

After  propagating  through  the  turbulent  atmo¬ 
sphere,  an  initally  coherent  wave  will  exhibit  a  finite 
coherence  area  and  a  flnite  coherence  time.  The  co¬ 
herence  diameter  ro  is  that  aperture  diameter  over 
which  the  rms  phase  fluctuations  ot  the  wavefront  are 
rad.  With  a  slow  angle  tracker,  this  parameter  is 
given  by* 

r„-1.68i**jrdliC»(«l-»«,  (1) 
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Pig.  1.  Site  of  the  Mark  III  stellar  interferometer  at  Mt.  Wilson. 


where  h  is  a  wavenumber,  is  the  structure  parame¬ 
ter,  and  the  integration  is  over  the  propagation  path. 
The  spatial  variance  of  the  phase  fluctuations  across 
the  interference  fringe  is  approximately  twice  that  of 
the  fluctuations  across  a  single  wavefront.  The  coher¬ 
ence  time  To  is  that  time  interval  during  which  the  rms 
phase  fluctuations  of  the  interference  fringe  are  equal 
to  1  rad.  With  a  slow  fringe  tracker,  and  for  the  case  of 


a  long-base  line  interferometer,  this  parameter  is  given 
by* 

1.361*^  Jf^"  (f/iC;(h)V^(h)I-®" 

=  0.81(r,/D.  V{h)  -  V,  (2) 

where  V  is  the  wind  speed.  Finite  values  for  ro  and  To 
restrict  the  maximum  usable  subaperture  size  and  the 
maximum  usable  coherent  integration  time  for  a  stel¬ 
lar  interferometer. 

While  atmospheric  phase  fluctuations  limit  the  sen¬ 
sitivity  of  a  ground-based  interferometer,  the  more 
serious  problem  for  astrometry  is  that  the  turbulent 
atmosphere  corrupts  the  angle  of  arrival  of  the  stellar 
wavefront,  i.e.,  the  apparent  position  of  the  star.  The 
statistics  of  the  angle-of-arrival  process  a(t)  severely 
limit  the  accuracy  possible  from  a  single-color  absolute 
astrometric  measurement.  The  angle  of  arrival  a(t)  is 
related  to  fringe  position  x(t)  via  x(t>  *  a(t)B,  where  fi 
is  the  interferometer  base  line  and  to  phase  difference 
4(f)  via  4(t)  kx(t),  where  k  >  2t/X  and  X  is  the 
wavelength  of  observation. 

The  most  useful  description  of  the  angle-of-arrival 
process  is  its  power  spectrum  Waif).  In  the  geometric- 
optics  approximation,  the  contribution  to  the  power 
spectrum  attributable  to  each  atmospheric  layer  can 
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be  computed  independently^;  the  total  spectrum  is  the 
sum  of  all  such  terms,  taking  into  account  the  height 
dependence  of  and  V.  Using  a  spectral  representa¬ 
tion,^  the  power  spectrum  W„Q)  can  be  written 

Wjn  -  r  dhV-'\  du 

/O  JO 

X  *  ((u*  +  (u>/V0^1''^)«(u,<i)),  (3) 

where  B  and  ^  are  the  base  line  and  wind  vectors,  $(k) 
is  the  refractivity  spatial  spectrum,  and  <o  =  2irf.  The 
filter  function  g(u,a))  is  given  by  2[1  -  cos(S<<j/ V)]  when 
V  and  B  are  parallel  ( V  ||  B)  and  by  2[1  -  co8(u)]  when 
they  are  perpendicular  (V 1  B).  Equation  (3)  is  valid 
in  the  geometric-optics  regime,  i.e.,  for  frequencies  /  « 
V/ylUL,  where  is  the  Fresnel  zone  size  and  L  is  the 
height  of  the  turbulent  layer.  The  representation  (3) 
can  also  be  used  to  calculate  the  power  spectra  of 
relative  and  time-lagged  relative  astrometric  measure¬ 
ments  through  the  appropriate  choice  of  the  function 

For  the  case  ^  ||  B  and  a  Kolmogorov  spatial  spec¬ 
trum,  ♦(«)  ■  0.033C^((~*’''^,  W„{fi  is  given  by® 

W,{/)  -  0.0654B-2  f*  'll  -  coa^2rfB/V)\^'*'^  (4) 

io 

where  the  nulls  in  Eq.  (4)  are  artifacts  of  the  frozen 
turbulence  assumption.  For  the  case  V  1.B,  the  spec¬ 
trum  Wa(/)  is  smooth  with  a  high-frequency  asymp¬ 
tote  equal  to  that  of  Eq.  (4)  and  a  low-frequency  as¬ 
ymptote  larger  than  that  of  Eq.  (4)  by  the  factor  of  %. 
Berauae  of  the  uncertainty  in  the  wind  direction  for 
moat  astronomical  applications,  Lindegren'®  suggests 
averaging  with  respect  to  wind  direction,  which  yields 
the  average  puwer  spectrum 

W„lf)  -  -  1.29  f*  dhCl\r'"r‘^'\  f  «  0.20V/B,  (5a) 

4  Jo 

JV.(/)  -  0.0664B-^  f*  /»0.20V/B.  (5b) 

Jo 

The  error  variance  of  a  T-second  absolute  astrome¬ 
tric  measurement  can  be  written 

r  dfWjn«nfywfT),  (6) 

Jo 

where  sincix)  ■  sin(x)/x.  For  T  »  B/V,  Eq.  (6)  can  be 
written*® 

2.92  dhClV-''^^T''\  T»B/V.  (7) 

For  the  case  of  l-sec  of  arc  seeing  and  a  constant  10-m/ 
sec  wind  speed,  the  value  of  / Cjd/i  is  ~8.4  X 10“*®  m*'®, 
and  the  error  standard  deviation  can  be  written 

«  0.257^''*  *ec  of  arc,  T»B/V.  (8) 

These  results  are  independent  of  the  interferometer 
base  line  B  and,  implicit  in  the  geometric-optics  ap¬ 
proximation,  independent  of  wavelength.  In  addi¬ 
tion,  these  results  depend  only  weakly  on  integration 
time  and  wind  speed,  a  consequence  of  the  nearly  Uf 
nature  of  the  atmospheric  statistics.  Very  long  inte¬ 


gration  times  are  thus  required  to  obtain  high  absolute 
accuracies. 

If  a  finite  value  for  the  outer  scale  Lo  is  incorporated 
into  the  refractivity  spatial  spectrum,  the  astrometric 
error  standard  deviation  in  an  absolute  measurement 
will  decrease  more  quickly  than  7^*/s  fg,  iQQg 
integration  times.  For  the  case  ^  ||  B,  and  a  von 
Karman  spatial  spectrum,**  ♦(«)  *  0.033C®(<t®  + 
the  angle-of-arrival  power  spectrum  (3)  is 

given  by'® 

W„(f)  =  1.29  r  dhCy-''Y^-^"' 

Jo 

X  |I  +  (V’/2ir/Lu)-)-*^’,  /«0.23V'/B.  (9) 

The  low-frequency  asymptote  predicted  by  Eq.  (9) 
makes  strict  demands  on  the  validity  of  the  frozen 
turbulence  assumption,*®  and  this  behavior  would  not 
be  expected  to  be  observed  in  practice. 

A  more  conservative  result  of  a  flat  low-frequency 
asymptote  is  obtained  for  the  case  V  ±  B: 

W„(/)=6.6  I'dhCiy-' 

Jo 

X  |L,;®  +  (2ir//V)=l-''\  /«0.18V/B:  B«Lo  (10a) 
iWJ/}  =  6.6  I"  dJiC^y-'  B-"- 

X  +  (2)r//V0'^)--'''':  B»4.  (10b) 

Averaging  with  respect  to  wind  direction  as  in  Eq.  (5) 
yields  low-frequency  asymptotes  equal  to  approxi¬ 
mately  one-half  of  the  vidues  given  in  Eq.  (9),  i.e., 

W',.(/)  =  3.3 

Jo 

f«  ynwl^;  B«Lq,  (lla) 

W,Af)  »  3.3 1^'  dhC^V-'  (Lo/B?L^, 

f«V/2irL^  B»L„.  (11b) 

Division  of  these  results  by  2T  yields  the  error  vari¬ 
ance  of  an  absolute  astrometric  measurement  for  inte¬ 
gration  times  long  compared  with  Lo/V.  Thus,  with  a 
finite  outer  scale,  the  error  standard  deviation  de¬ 
creases  as  white  noise  lOr  very  long  integration  times. 
In  addition,  the  standard  deviation  is  smaller  with  a 
base  line  much  longer  than  the  outer  scale.  Unfortu¬ 
nately,  it  is  likely  tlmt  such  base  lines  would  be  imprac- 
tically  large  and  would  also  resolve  most  of  the  stars  of 
interest.  Also,  the  base  line  necessary  for  a  significant 
improvement  in  accuracy  is  dependent  on  the  precise 
modeling  of  the  low-spatial-frequency  region  of  the 
refractivity  spectrum.  Using  a  Greenwood-Tarazano 
spatial  sp^rum*®'*^  rather  than  a  von  Karman  spec¬ 
trum  to  model  this  region  yields  a  somewhat  more 
pessimistic  result  than  that  presented  in  Eq.  (11b).® 

It  seems  highly  likely  that  the  appropriate  value  for 
the  outer  scale  for  vertical  propagation  through  the 
atmosphere  is  significantly  larger  than  the  several- 
hundred-meter  value  often  assumed.  Dravskikh  and 
Pinkelstein*®  estimate  that  the  outer  scale  should  be  of 


4108  APPUEO  OPTICS  !  Vol.  26.  No.  19  /  1  Octobw  1987 


the  order  of  the  scale  height  of  the  troposphere,  and 
that  below  this  spatial  frequency  the  3-D  spectram 
should  behave  as  k~^,  a  considerably  more  pessimistic 
result  than  either  a  von  Karman  or  Greenwood-Tara- 
zano  model  From  a  number  of  observations  of  the 
phase  structure  function  over  base  lines  from  1  to  10 
km  obtained  with  the  VLA  radio-interferometric  ar¬ 
ray,  Armstrong  and  Sramek'^  estimated  a  mean  spatial 
spectrum  obeying  k~^  *  for  those  observations  during 
which  turbulent  fluctuations  were  detectable.  While 
the  estimated  *  power  law  represents  a  slight  depar¬ 
ture  from  a  Kolmogorov  spectrum  at  these  base  lines, 
their  measurements  do  not  appear  to  indicate  an  outer 
scale  signiflcantly  smaller  than  that  estimated  above. 
While  radio  measurements  are  primarily  sensitive  to 
water  vapor  fluctuations,  as  opposed  to  temperature 
fluctuations  for  most  optical  measurements,  it  is  ex¬ 
pected  that  the  outer  scale  for  temperature  turbulence 
should  be  at  least  as  great  as  the  outer  scale  for  water- 
vapor  turbulence. 

The  discussion  above  concerns  only  a  single  one- 
color  absolute  astrometric  measurement.  A  signifr- 
cant  improvement  in  error  performance  is  obtained  for 
relative  measurements.  The  power  spectra  for  such 
measurements  have  flat  low-frequency  asymptotes,'** 
even  assuming  an  inHnite  outer  scale,  due  to  the  rejec¬ 
tion  of  low-spatial-frequency  atmospheric  energy. 
Relative  astrometry  with  a  single-beam  interferometer 
is  accomplished  by  star  switching,  in  which  the  posi¬ 
tions  of  a  number  of  stars  are  repeatedly  measured 
over  a  number  of  switching  cycles.  For  such  measure¬ 
ments,  the  error  in  the  difference  of  two  star  positions, 
or  in  the  difference  between  one  star  position  and  the 
centroid  of  the  other  positions,  should  decrease  as  the 
square  root  of  the  number  of  switching  cycles.  Since 
the  absolute  error  has  only  a  weak  dependence  on 
integration  time,  short  dwell  times  per  star  per  switch¬ 
ing  cycle  are  acceptable  to  maximize  the  number  of 
cycles.  In  addition,  short  switching  cycles  are  desir¬ 
able  to  reject  instrumental  instabilities  with  time 
scales  longer  than  the  cycle  time.  Finally,  it  should  be 
noted  that  a  significant  improvement  in  the  error  be¬ 
havior  of  either  an  absolute  or  relative  measurement 
can  be  obtained  with  the  two-color  method.-' 

■.  Data  from  the  Mlark  ■  and  Maik  ■  Mailoromatars 

The  data  from  the  Mark  III  interferometer  to  be 
presented  below  were  obtained  in  the  fall  of  1986  at  Mt. 
Wilson.  Common  instrumental  parameters  for  these 
observations  were  a  12-m  N-S  b^  line,  2.5-cm  diam 
collecting  apertures,  and  a  4-m8ec  coherent  integra¬ 
tion  time  for  the  red  channel.  The  effective  wave¬ 
length  in  the  red  channel  for  these  observations  was 
~0.7  fim.  The  observations  below  are  of  bright  stars 
which  provided  at  least  500  photons/4  msec  at  a  visibil¬ 
ity  of  M).6.  The  high  SNR  simplified  the  phase  un¬ 
wrapping  necessary  to  estimate  the  high-frequency 
atmospheric  phase  fluctuations.  The  Hrst  data  set  is 
from  a  short  scan  of  a  Aql  starting  at  2.8-h  UT  on  6 
Nov.  1986.  The  second  data  set  is  from  a  3-h  observa- 
tionofdCasstartingat2.8-hUTon25Nov.  1986.  The 
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Fig.  3.  Fringe  position  for  a  Aql,  6  Nov.  1986. 


third  data  set  is  from  a  1.4-h  observation  of  0  Gem 
starting  at  12.2-h  UT  on  30  Oct.  1986.  These  data  are 
excerpted  from  a  much  larger  data  set  designed  to 
evaluate  instrumental  stability  and  demonstrate  rapid 
star  switching;  these  astrometric  observations  are  re¬ 
ported  separately.''* 

The  Mark  111  interferometer  employs  active  fringe 
tracking.  The  white-light  fringe  position  x(t)  is  given 
by 

x(t)  =  ut)  ■¥  ifUUk,  (12) 

where  L{t)  is  the  position  of  the  laser-monitored  delay 
line,  Mt)  is  the  fringe  phase,  and  k  =  2r/X.  The  fringe 
phase  <t>  is  the  tracking  error  of  the  white-light  fringe 
servo.  It  is  measured  via  a  path-length  modulation 
scheme  and  used  to  control  the  delay-line  position. 
The  following  three  figures  illustrate  the  operation  of 
the  fringe  tracker. 

Figure  3  is  a  time  domain  plot  of  fringe  position  for  a 
Aql  at  a  resolution  of  100  msec/point  after  subtraction 
of  the  systematic  sidereal  motion.  Peak-to-peak  ex¬ 
cursion-  over  ~72  sec  are  ~16  jim.  This  level  of  fluctu¬ 
ation,  corrected  for  the  high-pass  Hltering  of  the  finite 
segment  length,  corresponds  to  ~0.3-sec  of  arc  seeing 
for  this  segment  of  data.  In  general,  the  variance  of 
the  fringe  position  fluctuations  would  be  expected  to 
be  given  b^ 

o'i  =  e.SSk--[B/r„f'\  (13) 

Figure  4  is  a  plot  of  fringe  phase,  i.e.,  the  error  of  the 
fringe  servo.  These  data  behave  as  a  highpass-filtered 
version  of  the  fringe  position  and  exhibit  rms  fluctua¬ 
tions  of  ~0. 1  fim.  A  clearer  indication  of  the  operation 
of  the  fringe  servo  is  provided  by  the  power  spectra  of 
these  data  presented  in  Fig.  5.  These  sp^ra  are 
consistent  with  a  first-order  servo  with  a  closed-loop 
bandwidth  of  20  Hz,  as  implemented  by  the  real-time 
system. 

The  next  figure  is  for  the  3-h  track  of  0  Cas.  These 
data  consisted  of  a  number  of  separate  fringe  locks 
with  a  somewhat  lower  duty  cycle  toward  the  end  of  the 
data  due  to  changing  seeing.  Again,  a  power  spectrum 


1  October  1987  /  Vo(.  26,  No.  19  /  APPLIED  OPTICS  4109 


Fig.  4.  Fringe  servo  error  for  a  Aql,  6  Nov.  1986. 


Fig.  5.  Power  spectra  of  fringe  position  and  servo  error  for  a  Aqi,  6 
Nov.  1986. 

provides  the  moet  information  about  the  data,  and  one 
is  presented  in  Fig.  6.  Asymptotes  corresponding  to 
and  power  laws  are  indicated,  and  good  corre¬ 

spondence  to  a  Kolmogorov  spatial  spectrum  over  the 
0.001-100-Hz  frequency  range  is  observed.  Good 
agreement  with  the  Kolmogorov  model  has  also  been 
obtained  by  others  in  laser  propagation  experiments 
over  short  horizontal  paths  ( <1M  m)  with  base  lines  of 
<30  cm.**-'® 

From  the  frequency  of  intersection  of  the  high-  and 
low-frequency  asymptotes  of  Fig.  6,  an  effective  wind 
speed  of  14  m/sec  is  estimated  using  Eq.  (5).  From  the 
high-frequency  asymptote  W„(f)  «  1.1  X  a 

value  of  ^e  coherence  time  To  15  msec  at  the  effec¬ 
tive  wavelength  of  0.7  urn,  or  11  msec  at  0.55  /nn,  is 
estimated  using  Eq.  (2).  By  comparison,  for  tlw  data 
of  Fig.  5,  the  coherence  time  is  estimated  as  18  msec  at 
0.55  iim. 

The  low-frequency  asymptote  of  Fig.  6  is  given  by 
W„{f)  ■  2.0  X 10“'®/“®^,  valid  for  frequencies  between 
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Fig.  6.  Power  spectrum  of  fringe  position  for  0  Cas,  25  Nov.  1986. 

~0.001  and  0.1  Hz  Using  this  value  and  the  effective 
wind  speed  of  14  m/sec  in  Eqs.  (5)  and  (1)  yields  a  value 
of  the  coherence  diameter  ro  25  cm  at  0.7  nm,  or  19 
cm  at  0.55  pm,  for  seeing  (arA/ro)  of  0.5  sec  of  arc. 
Using  Eq.  (7),  this  value  of  the  low-frequency  asymp¬ 
tote  corresponds  to  an  atmospherically  limited  abso¬ 
lute  astrometric  accuracy  of  0.147^*^  sec  of  arc  for 
long  integration  times  or  M).065  sec  of  arc  for  a  100-sec 
integration  time.  This  value  is  somewhat  smaller  than 
the  value  estimated  in  Eq.  (8),  which  assumed  worse 
seeing,  and  should  decrease  as  the  square  root  of  the 
number  of  switching  cycles  for  a  star-switched  relative 
measurement. 

At  frequencies  below  10~®  Hz,  fluctuations  with  a 
power  spectrum  much  steeper  than  are  observed 
in  the  data  of  Fig.  6.  These  can  be  identified  as  ran¬ 
dom  instrumental  errors,  especially  thermal  drifts. 
Prior  experiments  investigating  the  thermal  stability 
of  the  instrument  characterized  these  instabilities  as 
obeying  power  laws  between  and  f~*  at  very  low 
frequencie  s.  Rapid  star  switching,  with  cycles  times 
much  shorter  than  1000  sec,  should  be  adequate  to 
eliminate  astrometric  errors  due  to  these  instabilities. 
It  is  possible  that  some  fraction  of  these  errors  are  due 
to  imperfections  in  the  siderostat  bearings,  an  error 
which  cannot  be  reduced  by  star  switching.  Addition¬ 
al  laser  metrology  of  the  siderostats  is  required  to 
eliminate  this  error  source.  Such  a  system  should  be 
installed  on  the  Mark  III  instrument  in  1987.® 

At  frequencies  as  low  as  10~®  Hz,  good  correspon¬ 
dence  is  observed  to  a  simple  Kolmogorov  model,  i.e., 
one  with  an  infinite  outer  scale.  Using  this  frequency 
with  a  14-m/sec  wind,  Eq.  (10)  places  a  lower  bound  on 
the  value  of  the  outer  scale  at  '>^2  km.  As  mentioned 
earlier,  it  is  likely  that  the  true  value  is  still  somewhat 
greater  than  this.  Hopefully,  future  measurements 
with  the  Mark  III  interferometer  will  provide  a  better 
estimate  of  this  parameter. 

Figure  7  presents  a  power  spectrum  of  fringe  posi¬ 
tion  for  0  Gem.  Behavior  similar  to  that  of  Fig.  6  is 
observed,  with  an  identical  high-frequency  asymptote, 
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Fig.  7.  Power  spectrum  of  fringe  position  for  0  Gem,  30  Nov.  1986. 


Fig.  8.  Composite  power  spectrum  of  fringe  position  obtained  with 
the  Mark  II  interferometer,  Aug.  and  Sept.  1983. 


and  a  slightly  larger  low-frequency  asymptote  of  '^3.5 
X  corresponding  to  ro  a  14  cm  at  0.55  Mm. 

Finally,  for  comparison,  Fig.  8  presents  a  composite 
power  spectrum  of  several  nights  of  data  obtained  in 
Aug.  and  Sept  1983  with  the  earlier  Mark  II  interfer¬ 
ometer  at  Mt  Wilson.  Instrumental  parameters  were 
similar  to  Uiose  for  the  Mark  III  instrument  with  the 
ezc^ition  of  a  3.1-m  N-S  base  line.  Low-  and  high- 
firequency  asymptotes  corresponding  to  a  Kolmogorov 
spectnun  are  in^cated.  Th^  data  are  characterized 
values  of  rp  and  To  of  ~25  cm  and  24  msec  at  0.56 
Mm,  with  an  effective  wind  speed  oi~8.6  m/sec.  Ther¬ 
mal  instabilities  in  the  instiument  are  worse  than  for 
the  Mark  m  and  are  evident  at  frequencies  below 
~10~‘  Hz.  These  arose  because  the  Nfark  n  instru- 
mmit  was  constructed  on  a  steel  optical  table  exposed 
to  the  nighttime  air. 


IV.  Conclusions 

A  phase-coherent  stellar  interferometer,  such  as  the 
Mark  III  interferometer,  in  addition  to  being  a  power¬ 
ful  astronomical  instrument,  is  also  capable  of  provid¬ 
ing  high-resolution  measurements  of  the  atmospheric 
phiue  process.  Measurements  of  phase  difference  at  a 
12-m  base  line  obtained  with  the  Mark  III  interferome¬ 
ter  agree  well  with  the  predictions  of  a  Kolmogorov 
spatial  spectrum  over  the  0.001-100-Hz  frequency 
range.  On  the  basis  of  these  measurements,  tiie  outer 
scale  for  vertical  propagation  through  the  entire  atmo¬ 
sphere  does  not  appear  to  be  smaller  than  a  few  kilome¬ 
ters.  Measurements  at  Mt.  Wilson  indicate  a  one- 
color  atmospherically  limited  absolute  error  standard 
deviation  of  ~0.065  sec  of  arc  for  a  100-sec  integration 
in  conditions  of  0.5-sec  of  arc  seeing.  This  error 
should  decrease  as  the  square  root  of  the  number  of 
switching  cycles  for  a  star-switched  relative  measure¬ 
ment.  Greater  improvement  should  be  possible  with 
two-color  techniques. 
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Two-color  method  for  optical  astrometry:  theory  and 
preliminary  measurements  with  the  Mark  IH 
stoHar  Interferometer 

M.  Mark  Colavita,  Michael  Shao,  and  David  H.  Staelin 


The  two-color  method  for  interferometric  aatrometry  provides  a  means  of  reducing  the  error  in  a  stellar 
position  measurement  attributable  to  atmospheric  turbulence.  Tbe  primary  limitation  of  the  method  is 
shown  to  be  turbulent  water  vapor  fluctuations,  which  limit  the  amount  of  improvement  over  a  one-color 
measurement  obtainable  arith  a  two-color  estimate.  Secondary  atmospheric  effects  caused  by  diffraction 
from  small  refractive-index  inhomogeneities  and  differential  refraction  for  the  observation  of  stars  away  from 
zenith  are  shosm  to  introduce  errors  that  behave  as  white  noise  and  which  should  usually  not  be  significant. 
Other  potential  error  sources  due  to  photon  noise,  systematic  instrumental  effects,  and  imperfect  data 
reduction  are  alto  considered.  Tbe  improvement  in  accuracy  possible  with  the  two-color  method  is  estimated 
as  a  factor  of  5-10  over  the  corresponding  one-color  measurement.  Some  preliminary  two-color  measure¬ 
ments  with  the  Mark  III  stellar  interferometer  at  Mt.  Wilson  are  presented  which  demonstrate  a  factor  of  ~5 
reduction  in  the  amplitude  of  the  atmospheric  fluctuations  in  a  stellar  position  measurement 


mVOQUCQOIl 

The  accuracy  of  conventional  ground-baaed  aatro- 
metric  meaaurementa  suffers  from  the  deleterious  ef¬ 
fects  of  atmospheric  turbulence.  Because  of  the  na¬ 
ture  of  the  low-frequency  differential  atmospheric 
phase  fluctuations,  which  behave  almost  as  1//  noise,^ 
very  long  integration  times  are  usually  required  to 
obt^  highly  accurate  measurements  of  stellar  posi- 
tions.^-3  However,  the  accuracy  of  such  long-duration 
measurements  is  typically  limited  by  systematic  in¬ 
strumental  errors.  For  the  case  of  relative  astro¬ 
metry,^  the  situation  is  somewhat  improved.  Howev¬ 
er,  the  ultimate  limitation  to  high  accuracies  remains 
the  turbulent  atmosphere.  Hie  application  of  the 
two-color  method  to  interferometric  astrometry^'^  pro¬ 
vides  a  means  for  reducing  tiie  atmospheric  contribu¬ 
tion  to  the  error  in  a  stellfu  position  measurement 

Two-color  techniques  for  the  precise  measurement 
of  long  paths  through  the  c^n  atmosphere  have  been 
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used  successfully  for  geodetic  ranging^ ’’  and  have  ap¬ 
plication  to  sateUite  ranging^  ^  and  altimetry.’*’  In  this 
case,  the  problem  which  the  two-color  techniques  ad¬ 
dress  is  the  determination  of  the  path-averaged  refrac¬ 
tive  index  needed  to  convert  an  apparent  measure¬ 
ment  to  the  proper  vacuum  distance.  In  astrometric 
applications  of  two-color  techniques,  the  interest  is  in 
the  measurement  of  angles.  Large  atmospheric  scales 
are  thus  rejected,  and  the  nonunity  atmospheric  re¬ 
fractive  index  does  not  introduce  errors.  In  this  case, 
the  two-color  method  is  employed  to  correct  the  errors 
caused  by  turbulent  fluctuations  in  refractive  index, 
which  usually  contribute  only  a  small  amount  to  the 
error  in  the  range-flnder  case. 

The  use  of  two-color  techniques  is  similar  in  both  the 
range-flnder  and  astrometric  cases.  For  astrometry, 
let and  az  be  the  stellar  positions,  i.e.,  the  angles  of 
arrival,  observed  at  wavelengths  Xi  and  Xz.  The  two- 
color  estimate  of  the  stellar  position  azc  can  be  written 

"zr  “  “i  ~  D(o2  -  O)),  (1) 

where  D  is  the  atmospheric  dispersion, 

D  •  ln(X,)  -  ll/|n(Xj)  -  n(X,)l,  (2) 

and  n(X))  and  n(Xz)  are  the  atmospheric  refractive 
indices  at  wavelengths  Xi  and  Xg.  For  Xi  >  0.7  jun  and 
Xz  ■  0.6  fim,  the  dispersion  is  ~87.1.  In  Eq.  (1),  the 
quantity  ai  is  a  measure  ofthe  true  stellar  position  plus 
^e  instantaneous  atmospheric  error,  whUe,  in  princi¬ 
ple,  Z>(az  -  ai)  is  a  measure  of  the  atmospheric  error 
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only.  Thus  the  quantity  a^c  should  be  a  measure  of 
the  true  stellar  position. 

In  practice,  the  two-color  estimate  (1)  will  be  affect¬ 
ed  by  photon  noise  in  the  measurement  of  the  angle  of 
arrival  a  «  x/fi,  where  x  is  the  fringe  position  and  B  is 
the  base  line  length.  Let  be  ^e  photon-noise 
variances  in  the  measurement  of  fringe  phase  at  wave¬ 
lengths  Xi,X2  during  a  coherent  integration  time  of 
length  Tc.  Then  the  astrometric  error  variance  of  a  T- 
second  two-color  measurement  due  only  to  photon 
noise  can  be  written  for  large  D  as 

4^  “  (§y  (3) 

where  A  is  a  wave  number  2t/X.  The  photon-noise 
error  is  minimized  for  long  base  lines  and  small  values 
of  D,  the  latter  requiring  both  a  large  separation  be¬ 
tween  the  two  colors,  and  also  that  one  of  the  colors  be 
located  in  the  blue  region  of  the  spectrum,  where  the 
atmosphere  is  most  dispersive.  To  estimate  the  mag¬ 
nitude  of  the  photon-noise  error,  let  Xi  «  0.7  /im  u>d  X2 
»  0.5  pm,  for  which  D  “  87.1.  I^t  aj*  0.3*  rad*, 
which  corresponds  to  '^100-photons/coherent  integra¬ 
tion  time  for  a  fringe  of  0.5  visibility.  With  fi  «  12  m 
and  Tr «  4  msec,  the  rms  error  in  a  100-sec  integration 
due  to  photon  noise  is  0.4  msec  of  arc;  for  longer  obser¬ 
vations,  this  error  decreases  as  the  square  root  of  the 
integration  time  according  to  Eq.  (3).  Thus  despite 
the  luge  multiplier  D  in  Eq.  (1),  the  photon-noise  error 
is  quite  small  in  even  a  short  two-color  observation. 

However,  there  are  other  errors  in  the  two-color 
estimate  attributable  to  atmospheric  turbulence 
which  are  usually  more  significant  than  the  photon- 
noise  error.  The  primary  atmospheric  error  is  due  to 
turbulent  water  vapor,  which  obeys  a  different  disper¬ 
sion  relation  than  dry  air.  Essentially,  turbulent  wa¬ 
ter  vapor  introduces  an  error  which  has  the  appearance 
of  low-level  leakage  of  the  one-color  atmospheric  error 
into  the  two-color  measurement.  Secondary  atmo¬ 
spheric  effects  include  diffraction  from  small  refrac¬ 
tive-index  inhomogeneities  and  differential  refraction 
for  the  observation  of  stars  away  from  zenith.  These 
effects  will  usually  be  smaller  than  the  effect  due  to 
water  vapor  and  in  addition  decrease  as  the  square  root 
of  the  integration  time.  Section  11  considers  the  ef¬ 
fects  of  turbulent  water  vapor,  while  Secs.  Ill  and  IV 
consider  the  effects  of  diffraction  and  refraction.  Sec¬ 
tion  V  presents  some  requirements  on  the  quality  of 
the  instrument  and  on  the  data  reduction  necessary  to 
yield  accurate  two-color  position  estimates.  Finally, 
Sec.  VI  presents  some  preliminary  two-color  observa¬ 
tions  from  the  Mark  III  astrometric  interferometer  at 
Mt.  Wilson.  These  observations  show  a  factor  of  5 
reduction  in  the  amplitude  of  the  atmospheric  fluctua¬ 
tions  in  an  angle-of-arrival  measurement. 

vJipOv  CWOV 

Atmospheric  turbulence  is  driven  by  microfluctua¬ 
tions  in  temperature,  pressure,  and  water  vapor  with 
temperature  fluctuations  making  the  dominant  contri¬ 
bution  to  refractivity  fluctuations  in  the  optical  region. 


The  two-color  method  can  eliminate  the  astrometric 
errors  caused  by  microfluctuations  in  temperature  and 
pressure;  however,  the  method  cannot  simultaneously 
eliminate  the  errors  caused  by  microfluctuations  in 
water  vapor. 

At  optical  frequencies,  the  atmospheric  refractivity, 
N(\)  •»  n(X)  -  1,  can  be  written” 

N(Xt  =  AMP/T  -  B(X)e/T,  (4) 

where  P  is  total  pressure  in  millibars,  e  is  the  partial 
pressure  of  water  vapor  in  millibars,  and  T  is  tempera¬ 
ture  in  kelvins.  The  wavelength  dependence  of  the 
refractive  index  is  given  by  the  dispersion  relations 
A(X)andB(X)'': 

,4(X)  =  10-"|2372.4.34  +  684255.24/(130  -  a*) 

+  4549.40/(.38.9  -  »*)),  (5a) 

W(X)  =  4(X)  -  10'*|64S7.31  +  58.058ir* 

-  0.71  ISOa*  +  0.08851ff*(,  (5b) 

where  a  =  1/X  in  i/m"'.  At  X  *  0.55  ^m,  Eq.  (4)  can  be 
written 

(V  =  I0'“(79.01f>/7’  -  12.26e/T).  (6) 

Linearizing  Eq.  (6)  for  small  parameter  fluctuations 

and  using  the  relation  Q  2l7e/T  between  water 
vapor  partial  pressure  and  water  vapor  density  Q,  in  g/ 
m*.  yield’* 

AN  <=  -MlXlB/T-lAr  +  M(X)/7lAP  -  |B(X)/217)AQ 

- -a(X)Ar+6(X)AB-c(X)Ae-  (7) 

At  optical  wavelengths,  this  expression  is  usually  dom¬ 
inate  by  temperature  fluctuations.  At  X  *  0.55  urn,  P 
»  1013.25  mbar,  and  T  »  288.16  K,  the  values  of  the 
parameters  a,  b,  and  c  in  Eq.  (7)  are  0.96  X  10~®,  0.27  X 
10-8,  and  0.056  X  lO”®. 

Consider  now  the  two-color  measurement  given  by 
Eq.  (1).  In  the  geometric-optics  approximation,  angle 
of  arrival  is  proportional  to  the  difference  in  integrate 
refractivity  between  the  two  interferometer  beams. 
With  this  interpretation  of  AN,  the  instantaneous  an¬ 
gle  of  arrival  a  at  wavelength  X,  can  be  written 

-a(X,)Ar  +  b(X,)aP  -  c(X,)AQ.  (8) 

Thus  the  two-color  measurement  (Eq.  (Dj  can  be  writ¬ 
ten 

r  -fi(.\|>(i  -  /)/o,)ar  +  b(x,)(i  -  n/n,.\AP 

-  i  -  n/n^\AQ.  (9) 

where 

1>I  =  n,.  •  /1(X,)/|/4(X.,)  -  .4(X,)I,  (10a) 

11, ,  =  «(X,)/|«(X.,)  -  ft(X,)l,  (10b) 

with  A(X),B(X)  as  given  in  Eqs.  (5).  Evaluation  of  Eqs. 
(10)  at  Xi  =  0.7  Mm  and  X^  ■  0.5  am  yields  Dr  ^  Dp 
87.1  and  D(^  ^  -61.1.  Thus,  while  the  dispersion 
constants  for  temperature  and  pressure  fluctuations 
are  equal,  the  dispersion  constant  for  water  vapor  fluc¬ 
tuations  is  different  from  both  of  these  and  of  the 
opposite  sign.  Thus  the  effects  of  temperature,  pres- 
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sure,  and  water  vapor  on  the  refractive  index  cannot  be 
corrected  simultaneously  when  only  two  wavelengths 
are  observed. 

At  visible  wavelengths,  Eq.  (9)  will  usually  be  mini¬ 
mized  by  choosing  D  »  Dr,  leaving  the  residual 

-c(X|)(l +|D^O,,|)AQ.  (10 

Thus  the  ratio  of  the  variance  of  the  error  in  a  two-color 
measurement  to  the  variance  of  the  error  in  a  one-color 
measurement  can  be  written 


(4)  ^  1  + 1  ?r  I Y  r  W'*  1 


(J2) 


Using  the  previous  values  of  a,  c,  Dt,  and  Dq,  the 
coefficient  of  the  structure  constant  ratio  in  Eq.  ( 12)  is 
equal  to  0.021. 


The  ratio  of  variances  in  Eq.  (12)  could  just  as  well  be 
a  ratio  of  power  spectral  densities,  since  the  spectral 
behavior  of  temperature  and  water  vapor  fluctuations 
is  identical.*^-'^  Thus  the  residual  error  in  a  two-color 


measurement  will  behave  like  the  error  in  a  one-color 


measurement,  reduced  in  energy  according  to  Eq.  (12). 
It  should  be  noted  that  a  three-color  technique  could 
simultaneously  correct  errors  due  to  temperature, 
pressure,  and  water  vapor  fluctuations.  However, 
with  all  three  colors  in  the  optical  region,  the  appropri¬ 
ate  dispersion  constants  are  of  the  order  of  10-'’-10^, 
and  implementation  would  be  somewhat  difficult. 

To  estimate  the  improvement  in  accura<^  possible 
with  the  two-color  method,  the  value  of  the  simulta¬ 


neous  ratio  of  the  integrated  water  vapor  turbulence  to 
the  integrated  temperature  turbulence  is  required. 
However,  from  available  simultaneous  measurements, 
only  the  point  ratio  at  low  altitudes  can  be 

determine.  Since  the  water  vapor  scale  height  is  only 
'»2  km  vs  km  for  the  dry  atmosphere,  these  mea¬ 
surements  should  provide  a  conservative  estimate  of 
the  integrated  ratio.  While  none  of  these  measure¬ 
ments  was  performed  at  night  at  an  observatory-quali¬ 
ty  site,  the  measurements  of  (>ossard‘‘  seem  most  ap¬ 
plicable.  These  measurements  were  made  slightly 
inland  from  the  Pacific  coast  at  an  elevation  of  several 


hundred  meters.  From  comparison  of  the  simulta¬ 
neously  measured  temperature  and  humidity  spectra, 
a  value  of  0.4  is  determined  for  the  ratio  The 

values  of  this  ratio  estimated  from  other  inland  mea¬ 


surements**-**  range  from  0.06  to  0.67. 

Some  additional  measurements  of  interest  are  those 


of  Armstrong  and  Sramek**  of  the  atmospheric  phase 
structure  function  as  measured  with  the  VLA  for  base 


lines  between  1  and  10  km.  Approximating  their  em¬ 
pirical  structure  functions  with  a  Kolmogorov  model 
matched  at  a  1-km  separation  yields  a  mean  value  for 
the  integrated  radio  refractivity  of  / C^(radio)dA  «  3.4 
X  10~*2  m*^.  This  mean  value  applies  to  approxi¬ 
mately  two-thirds  of  their  observations;  there  was  no 
clear  detection  of  the  atmosphere  in  their  remaining 
observations,  which  appeand  to  have  a  detection 
threshold  of  approziiutely  /(]^(radio)dh  -  10~*^ 
m*^.  A  typical  ^ue  of  the  intei^ted  optical  refrac- 
tivity  can  Iw  estimated  using,  for  example,  a  Hufhagel 


turbulence  profile,'*'  -'  which  corresponds  to  ~1  sec  of 
arc  seeing,  yielding  /C^(optical)dh  *  6.2  X  10"' *  m'^*. 
At  radio  frequencies,  assuming  dominant  humidity 
fluctuations,  U^iradio)  ==  3.6  X  10"  "C^,  while  at  opti¬ 
cal  frequencies,  assuming  dominant  temperature  fluc¬ 
tuations,  C^(optical)  “  9.3  X  10"' The  ratio  of 
these  quite  different  nonsimultaneous  measurements 
yields  a  value  for  the  ratio  of  the  integrated  water 
vapor  turbulence  to  the  integrated  temperature  turbu¬ 
lence  of  0.15,  somewhat  smaller  than  the  point  mea¬ 
surement  cited  above. 

Thus,  using  the  above  parameters  in  Eq.  (12),  it  is 
estimated  that  a  factor-of-5-10  improvement  in  accu¬ 
racy  (standard  deviation)  over  a  one-color  measure¬ 
ment  should  be  possible  using  the  two-color  method. 
For  example,  the  error  in  a  one-color  absolute  astrome¬ 
tric  measurement  under  the  assumption  of  Kolmo¬ 
gorov  turbulence  and  1-sec  of  arc  seeing  is  of  the  order 
of  0.1  sec  of  arc  for  a  lOO-sec  integration.'’-^^  Thus, 
assuming  a  stable  instrument,  a  two-color  measure¬ 
ment  of  the  same  quantity  should  provide  an  accuracy 
of  0.01-0.02  sec  of  arc  for  the  same  integration  time. 
Similar  improvements  should  also  be  obtained  for  rela¬ 
tive  measurements  and  star-switched  measurements, 
which  have  smaller  one-color  atmospheric  errors.*'^^ 

While  water  vapor  fluctuations  are  the  primary  at¬ 
mospheric  limitation  of  the  two-color  method,  ^ere 
exist  several  smaller  sources  of  error  caused  by  atmo¬ 
spheric  turbulence;  these  are  discussed  in  the  next  two 
sections.  In  these  sections,  a  dry  atmosphere  is  as¬ 
sumed  so  that  there  is  no  error  due  to  turbulent  water 
vapor. 

M.  DMractive  Effects 

When  diffraction  is  considered  in  the  solution  of  the 
wave  equation  for  propagation  through  turbulence,  it 
is  found  that  phase  fluctuations,  observed  at  different 
optical  wavelengths,  become  decorrelated  when 
caused  by  small  refractive-index  inhomogeneities.^*-^ 
With  widely  spaced  wavelengths,  the  size  of  the  inho¬ 
mogeneities  below  which  decorrelation  occurs  is  of  the 
order  of  the  Fresnel  zone  size  where  X  is  the  mean 

wavelength  and  L  is  the  height  of  the  particular  turbu¬ 
lent  layer,  llius,  in  the  absence  of  water-vapor  fluctu¬ 
ations,  this  decorrelation  limits  the  effectiveness  of  the 
two-color  method  at  high  temporal  frequencies.  Re¬ 
lated  problems  arise  in  adaptive  optics  systems  when 
the  pointer  and  tracker  be^s  are  at  different  wave¬ 
lengths.'*'’** 

The  po’  ver  spectrtun  of  the  error  introduced  into  a 
two-color  astrometric  measurement  by  diffraction 
K’diiK/)  can  be  determined  as  follows.  LetS|,52bethe 
optical  phases  at  wavelengths  Xi,X2  observ^  at  one 
aperture  of  the  interferometer.  U  will  be  convenient 
to  neglect  atmospheric  dispersion  in  calculating  the 
statistics  of  the  error;  thus  let  >  Si  and  S^  >  [1  + 
(1/Z>)]S2,  where  D  is  the  dry-air  dispersion  and  81,82 
are  the  phases  which  would  1^  obeerWl  in  a  nondisper- 
sive  atmosphere.  Thus  the  path-length  error  d(t)  in  a 
two-color  distance  measurement  at  a  single  aperture 
can  be  written 
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d(t)  -  kj%u)  -  Dike's:, it)  -  *r‘s;(‘)i 

-  (D  +  l)(*r'S,(0  -  (13) 

where  k-, «  2t/X„  ki<k2.  The  error  cl(t)  has  the  power 
spectrum 

Wjifi  -(D  +  +  *J’'»Vs2S2(/) 

-2kr'ki‘Ws,siini.  (lO 

where  the  terms  Wsisjif)  are  the  optical  phase  power 
spectra  and  cross  spectra.  These  spectra  can  be  writ- 
ten» 

M'ssyi/)  -  8TV/J'2*.//v*d«l(«V0“  -  (2r/)T''^ 

X  f;dhgi»M*i'‘.h).  (15) 

where  gUJi)  is  the  spectral  Hlter  function  for  phase 
fluctuations 

-  C08[(l  -  k,/k,)h,y2/,,] 

+  CO«|(l  +  fc,/fc^)/llrV2*,J.  (16) 

In  these  equations,  V  is  the  transverse  wind  velocity 
and  4(K,h)  is  the  refractivity  spatial  spectrum. 

The  power  spectrum  Wjif),  Eq.  (14),  can  be  rewrit¬ 
ten  as  a  single  integral  of  the  form  of  Eq.  (15),  with 
filter  function 

gj(ljl)  "  |1  +  C08(/l<:V*,)|  +  II  +  C08((*,/*2))lJrV*,|l 

-  2|co*((l  +  *,/fc,)/i»2/2A|| 

+  Md  -  *,/*2)A-(V2*,||,  (17) 

which  can  be  rewritten  as 

gi(itji)  “  4»m^I(l  -  *,/fc2)li»'’/4A,| 

+  4  sin^Id  +  *|/*2)5«V4*,) 

-  2  8in'*I(*,/*2)(i»V2*,|  -  2  sin^(/i«-/2)i,).  (18) 

Each  term  of  the  filter  function  gd(K,/i)  has  the  same 
form  as  the  Hlter  function  for  log  intensity  fluctua- 
tions^ 

«i„,(«,)i)  =  8  8in^(/i«-/2*).  (19) 

Thus  it  is  expected  tha  the  power  spectrum  Wd(/)  will 
have  the  same  form  as  the  log  intensity  power  spec¬ 
trum;  in  particular,  the  spectrum  will  be  white  for 
frequencies  /  <(  V/^JXL,  which  implies  an  astrometric 
error  which  decreases  as  the  square  root  of  the  integra¬ 
tion  time.  Substituting  a  smoothly  varying  Kolmo¬ 
gorov  spectrum  into  Eq.  (15),  combining  Eqs.  (14), 
(15),  and  (18),  and  solving  for  the  low-frequency  as¬ 
ymptote  yiel(i 

Wj{f)  -  +  l)-/,T«d»(«V)-' 

X  f;dhgJ,J,)0.(KaCl(hU-''"  l«VI^\L.  (20) 

Evaluating  the  x  integral  yields 

W^(/)  -  2.0(0  +  xmktM;) 

X  |/JdWi<«C^()i)V-'(/i)l,  /«V/VXL.  (21) 

where  the  auxiliary  function  F{k\M)  is  given  by 

m,  w  2id  - /i,/)!2)/2i<'^' 

+  21(1  +  *  -  (*  i/*-.-)^'  -  1  •  (22) 


This  function  can  be  approximated  by  its  first  term  to 
an  accuracy  of  12%  for  wave  number  ratios  hi/hz  ^  %. 

The  power  spectrum  for  log  intensity  fluctuations  is 
computed  similarly,  using  the  filter  function  (19), 
yielding^® 

=  7.94*^'' 

X  l/JdWi^^C^OilV^'Oi)),  /«0.67V/V)1.  (23) 

Combining  Eqs.  (21)  and  (23)  allows  the  power  spec¬ 
trum  of  the  path-length  error  at  one  aperture  to  be 
written  in  terms  of  the  log  intensity: 

Wjlf)  -  0.25(0  +  l)‘^*r"^{A,yltj)»4'^,(0,*,),  K/vOT.  (24) 

Diffractive  effects,  acting  like  scintillations,  will  be 
uncorrelated  for  transverse  separations  in  excess  of  a 
few  Fresnel  zones.'  Thus  the  path-length  errors  at 
each  aperture  of  the  interferometer  can  be  considered 
essentially  independent,  so  that  the  error  in  angle  of 
arrival  for  an  interferometer  of  base  line  B  can  be 
written 

W'diffI/)  =  0..50(O  +  U^H-^krF(k„k2)W^^(0M,  /«  V/ypZ,  (25) 

where  W)„y(0,fci)  and  F{kj,k2)  are  given  by  Eqs.  (23) 
and  (22).  For  realistic  integration  times,  only  the  low- 
frequency  asymptote  of  Wdin(/)  is  need^.  However, 
the  high-frequency  asymptote  is  easily  computed  by 
assuming  uncorrelated  fluctuations  at  each  wave¬ 
length.  This  yields  a  high-frequency  power  spectnun 
proportional  to 

The  astrometric  error  variance  in  a  T-second  inte¬ 
gration  due  to  diffractive  effects  can  be  written 

4i.  H'j^ff(0)/27’.  T»J^/V  (26) 

for  an  error  standard  deviation  which  decreases  as  the 
square  root  of  the  integration  time.  Usi^  a  Hufnagel 
turbulence  profile,  the  integral  Sdhh*^^d~{h)  is  equal 
to  6.2  X  10~®  m'’/ ’.  Assuming  B  =  12  m,  V{h)  =  V  *  10 
m/sec,  and  X|  =  0.7  ^m,  X2  =  0.5  ^m,  for  which  D  =  87.1 
and  F  =  0.14,  the  low-frequency  asymptotes  (23)  and 
(25)  have  values  IFin/lO)  =  0.002  and  IFdifKO)  =  10"'®. 
These  yield  a  two-color  error  of  ~0.15  msec  of  arc  rms 
for  a  100-sec  integration.  This  should  be  a  conserva¬ 
tive  result  due  to  the  assumptions  of  point  receivers 
and  a  zero  inner  scale.  Thus  the  astrometric  error  in  a 
two-color  position  measurement  due  to  diffraction 
should  be  negligible. 

IV.  Relractiva  Effects 

A  first-order  analysis  of  the  two-color  method  as¬ 
sumes  that  the  beam  paths  in  each  color  are  identical. 
However,  for  observation  of  stars  away  from  zenith, 
dispersive  refractivity  gradients  transverse  to  the  di¬ 
rection  of  propagation  cause  the  beams  to  bend  and  the 
different  color  beams  to  separate.  Since  the  beams 
then  pass  through  different  turbulence,  the  effective¬ 
ness  of  the  two-color  correction  is  diminished.  The 
effect  just  described  is  purely  geometric.  It  occurs 
even  in  the  geometric-optics  limit  of  wave  propagation 
where  the  previously  derived  diffractive  effects  are 
nonexistent.  Related  problems  arise  in  two-color  geo- 
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detic  and  satellite  ranging;  previous  studies,^^  howev¬ 
er,  have  concerned  path-length  rather  than  angle-of- 
arrival  eCfects. 

A  full  wave-optics  solution  to  the  propagation  prob¬ 
lem  for  separate  two-color  beam  paths  shows  that  the 
astrometric  error  caused  by  the  combination  of  diffrac¬ 
tive  and  refimctive  effects  is  not  separable.  However, 
it  can  be  shown  that  the  effect  of  the  two  mechanisms  is 
bounded  above  by  the  rms  sum  of  the  contribution  due 
to  diffraction,  calculated  assuming  coincident  beam 
paths,  as  in  the  last  section,  plus  the  contribution  due 
to  geometric  beam  separation,  calculated  in  the  geo¬ 
metric-optics  propagation  Umit.^^  Therefore,  the  fol¬ 
lowing  a^ysis  will  be  carried  out  assuming  geometric 
optics. 

Light  from  a  star  at  zenith  angle  8  undergoes  a  sys¬ 
tematic  astronomical  refraction  which  displaces  the 
refracted  beam  from  an  imrefracted  beam  path.  As¬ 
suming  an  exponential  plane-parallel  atmosphere,  the 
transverse  displacement  of  these  beams  as  a  function 
of  height  h  above  ffie  observatory  can  be  written^'-^ 

4(/»)  "  /ioN(X,)  tan(A)  aedt)  e»p(-/ii/ho)  eip(-/i//io),  (27) 

where  *  n(Xi)  - 1  is  the  atmospheric  refractivity 
at  wavelength  Xi,  hi  is  the  height  of  the  observatory 
above  sea  level,  and  ho  8.3  km  is  the  scale  height  of 
atmospheric  density.  Since  the  atmosphere  is  disper¬ 
sive,  different  colors  undergo  different  amounts  of 
bending.  Assuming  that  the  different-color  beams  are 
coincident  at  the  earth’s  surface,  the  transverse  beam 
separation  as  a  function  of  height  can  be  written 

Mh)  -  (/ioA((X|)/Di  tan(9)  8ec(#) 

X  eip(-/i,/(i(|)(l  -  exp(-X//io)|.  <28) 

where  D  is  the  atmospheric  dispersion.  For  example, 
for  a  star  at  a  zenith  angle  of  15°,  a  sea-level  observa¬ 
tory,  and  a  dispersion  of  100,  the  beam  separation  at 
the  top  of  the  atmosphere  is  6.4  mm. 

Let  Wzif)  denote  the  power  spectrum  of  the  angle- 
of-arrival  error  caused  by  refractive  effects  as  calculat¬ 
ed  in  the  geometric-optics  limit.  Except  for  a  factor  of 
(D  -f  1)2,  the  desired  power  spectrum  has  the  same 
form  as  the  power  spe^rum  of  a  relative  astrometric 
measurement  Wr(J)  (Refs.  3,33): 

W^f)  =  10.5B"'''^/Jd/iC’(h)V-‘(h)ii'^(/i), 

f«0.nV/B;  0«B,  (29) 

where  A(h)  is  the  separation  of  the  beams  from  one 
aperture  as  a  function  of  height.  Therefore,  the  power 
spectrum  of  the  astrometric  error  will  be  white  at  low 
fluencies.  Substituting  Eq.  (28)  for  A(h)  and  let¬ 
ting  (D  +■  1)/D  at  1  yield 

W2(/)  -  10.5B-<«|/lo^(X,)  e»p(-/i,//io)|2tan2(«)8ec^«) 

X  J'J<«»c;(h)V-'(h)ll  -  m|3(-(i/M^  /«V/B,  (30) 

a  result  which  is  nearly  independent  of  wavelength. 
This  result  includes  the  factor  sec2(9)  to  account  for 
the  lengthening  of  the  propagation  path  and  a  worst- 
caae  reduction  of  transverse  wind  velocity  with  zenith 
angle. 


The  astrometric  error  variance  of  a  T-second  inte¬ 
gration  due  to  refractive  effects  can  be  written 

4  a  W;5(0)/2r.  T»B/V  (31) 

for  an  error  standard  deviation  which  decreases  as  the 
square  root  of  the  integration  time.  Using  a  Hufnagel 
turbulence  profile,  the  value  of  the  integral 
SdhC^{h)[l  —  exp(-h//»o)12  is  1.3  X  lO”*®  m*''®.  As  in 
the  last  section,  assume  R  -  12m,  V{h)  *  lOm/sec, 
and  Xi  »  0.7  ttm.  For  hi  »  1.74  km  (the  elevation  of 
Mt  Wilson),  Jie  astrometric  error  in  100-sec  integra¬ 
tion  is  0.5, 1.5,  and  3.8  msec  of  arc  rms  for  sources  at  15, 
30,  and  45°  from  zenith.  Again,  these  results  should 
be  conservative  due  to  the  assumptions  of  point  aper¬ 
tures,  a  zero  inner  scale,  and  geometric  optics.  Thus 
refractive  effects  should  be  negligible  for  observation 
of  stars  near  zenith,  and  even  for  sour<»s  at  45°,  the 
refractive  error  in  a  100-sec  integration  should  still  be 
at  least  a  factor  of  25  below  the  one-color  absolute 
astrometric  error.  Thus,  for  most  two-color  observa¬ 
tions,  refractive  effects  should  not  be  a  limiting  factor. 

V.  kwinimenlal/Data  Procaaaing  ConaMoraliona 

The  previous  sections  have  considered  fundamental 
sources  of  error  in  applications  of  the  two-color  meth¬ 
od  attributable  to  photon  noise  and  atmospheric  tur¬ 
bulence.  This  section  discusses  some  practical  consid¬ 
erations  relating  to  the  instrument  and  the  data 
processing  necessary  to  yield  accurate  two-color  posi¬ 
tion  estimates. 

A  potentially  serious  source  of  systematic  error  in  a 
two-color  measurement  is  due  to  mismatches  in  the 
thicknesses  of  dispersive  materials  in  each  arm  of  the 
interferometer.  L«t  t  denote  the  thickness  mismatch 
of  a  material  with  refractivity  N„  =  «„-!.  Then  the 
residual  error  in  a  two-color  fringe-position  measure¬ 
ment  at  wavelengths  Xj  and  Xz  can  1^  written 

=  (32) 

where  D  is  the  value  of  the  atmospheric  dispersion 
assumed  in  the  data  reduction,  and  D„  »  Nm(Xi)/ 
[iVmlXz)  -  iV,„(Xi)]  is  the  material  dispersion.  In  prin¬ 
ciple,  with  relative  measurements,  this  residual  error  is 
not  significant  ae  long  as  it  is  stable  with  time.  Howev¬ 
er,  with  a  wideband  instrument,  the  wavelengths  Xi 
and  Xz,  and  hence  the  dispersion  Dm,  will  be  dependent 
on  the  color  of  the  star.  To  estimate  the  magnitude  of 
fluctuations  in  ezc,  model  the  interferometer  channels 
as  rectangular  with  passbands  of  0.44-0.56  fun  for  blue 
and  0.56-0.86  for  red.  Furthermore,  assume  that 
the  starlight  can  be  modeled  with  a  blackbody  spec¬ 
trum  over  that  range  of  wavelengths.  Finally,  assume 
that  the  dispersive  material  is  fu^  silica,  llien,  for  a 
range  of  temperatures  from  2500  to  25,(X)0  K,  which 
models  most  of  the  stars  to  be  observed,  the  range  of 
D„  is  from  63.1  to  65.3.  Setting  6  equal  to  its  mean 
value  of  87.1,  the  fluctuations  in  ezc  over  this  range  of 
blackbody  temperatures  are  approximately  0.02t  P-P. 
To  keep  this  error  smaller  than  0.1  ^m  requires  t  <  5 
urn.  Fortunately,  this  degree  of  matching  is  readily 
obtained  by  adjusting  the  angle  of  a  compensating 
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plate.  With  the  instrument  configured  as  a  Fourier 
transform  spectrometer,  this  optic  is  adjusted  untU  the 
central  Mnges  in  both  channels  are  aligned  to  within 
approximate^  one-fourth  of  a  fringe. 

A  random  error  results  from  the  use  of  the  incorrect 
wavelength  or  dispersion  in  the  two-color  data  reduc¬ 
tion,  which  introduces  a  leakage  of  the  one-color  spec¬ 
trum  into  the  two-color  result.  Let  xi(t),X2(t)  be  the 
fringe  positions  at  wavelengths  Xi,X2.  Each  of  these 
positions  is  given  by  the  sum  of  the  position  of  the 
laser-monitored  delay  line  L  and  the  phase  of  stellar 
fringe  0,  i.e.,  x  >  L  +  Using  a  tilde  to  denote  the 
value  of  a  parameter  assmned  in  the  data  reduction,  let 
if,  ■  k,(l  -  ei),  where  «,  is  the  fractional  wavelength 
error  at  X,,  and  let  D  ■  D(1  +  to)-  Then,  to  first  order 
the  residual  error  in  a  two-color  measurement  can  be 
written 

«}(•«)  a  -(<2  +  “  '2^  “  <o^^tUy  (33) 

Tbe  processes  Ht)  and  behave  as  low-pass  and 
high-pass  Altered  versions  of  the  angle-of-arrival  pro¬ 
cess  a(t).  Thus  for  long  integration  times,  the  phase 
leakage  term  should  become  negligible,  and  tbe  astro- 
metric  error  of  the  two-color  measurement  should  ap¬ 
proach  that  of  a  one-color  measurement,  reduced  in 
amplitude  by  <2  fth  Therefore,  even  with  several 
percent  errors  in  D  and  k,  it  is  possible  to  realize  a 
greater  than  order-of-magnitude  improvement  in  as¬ 
trometric  accuracy  over  a  one-color  measurement. 

VI.  Two-Color  Obeervsflom  wWi  the  Marie  W 
Maiferonialar  at  Ml.  WSsoit 

The  Mark  III  astrometric  interferometer  is  a  phase- 
coherent  white-light  stellar  interferometer  designed 
for  astrometry.  Employing  at  present  a  12-m  N-S 
base  line,  the  instrument  became  operational  at  Mt. 
Wilson  in  Sept.  1986.  It  succeeds  the  Mark  II  instru¬ 
ment  which  operated  on  Mt.  Wilson  from  Aug.  1982 
until  Sept.  1984.  The  interferometer  measures  the 
phase  of  the  white-light  fringe  in  two  wide  spectral 
channels  covering  ~0.6-0.9  /im  and  0.4-0.6  The 
fringe  phase  in  the  red  channel  serves  as  the  error 
signal  for  a  white-light  fringe  servo  which  controls  the 
position  of  a  laser-monitored  delay  line,  while  the 
fringe  phase  in  the  blue  channel  is  used  off-line  with 
the  two-color  method.  A  coherent  integration  time  for 
the  red  channel  of  4  msec  is  employed,  during  which 
time  the  quadrature  components  of  the  fringe  phase  in 
each  channel  are  recorded,  along  with  the  delay-line 
position  and  other  information.  Continuous  fringe 
tracks  range  in  length  from  a  fraction  of  a  second  to 
tens  of  seconds,  depending  on  the  seeing.  The  instru¬ 
ment  is  discussed  in  detail  in  a  separate  paper  as  are 
some  one-color  atmospheric  phase  measurements” 
and  some  preliminary  astrometric  measurements.'*^ 

The  data  below  are  from  a  3-h  track  of  0  Cas  starting 
at  2.8-h  UT  on  25  Nov.  1986.  With  2.5-cm  diam  aper¬ 
tures,  photon  fluxes  of  500  and  650  per  4  msec  were 
observed  in  the  red  and  blue  channels,  respectively, 
with  average  fringe  visibilities  of  0.5  and  0.3.  The 
mean  wavelengths  in  the  red  and  blue  channels,  ob¬ 


tained  from  a  least-squares  At  to  several  seconds  of 
data,  were  Xi  =  0.704  ^m  and  X2  »  0.529  ^m  for  a  dry-air 
dispersion  D  »  110.  Good  seeing  characterized  most 
of  this  observation,  with  approximate  values  for  the 
coherence  diameter  and  coherence  time  of  ro  at  25  cm 
and  To  15  ms  at  0.7  /on,  correspondiog  to  ''^.5-Bec  of 
arc  seeing.  The  next  several  paragraphs  describe  the 
data  processing  applied  to  the  raw  data  collected  by 
the  r^-time  system. 

Every  4  msec,  during  a  fringe-lock,  the  phase  of  the 
white-light  fringe  in  the  red  chaimel  is  estimated  as 

=  tan-’S|(V', -X,)/(y,  +  X,)l,  (34) 

where  Xj  and  Vi  are  the  quadrature  components  for 
the  red  channel.  These  are  given  by  X  ■  C  -  A  and  Y 

D  —  B,  where  A,  B,  C,  and  D  are  the  four  bins  into 
which  photons  are  accumulated  during  a  fringe  scan.^ 
The  near-unity  constant  S  corrects  for  the  difference 
between  the  effective  red  wavelength  and  the  leng^  of 
the  path-length  modulation  stroke  s;  for  8  »  1,  ® 

tan'HVi/Xi)  -  t/4.  The  quantity  6  is  given  by 

=  (35) 

where  t  =  2t(s  -  \e)/K  is  the  number  of  radians  by 
which  the  path-length  modulation  stroke  exceeds  the 
effective  wavelength.  A  value  of  5  ^  1.17  was  used  in 
these  reductions  corresponding  to  s  0.633  Mm-  The 
fringe  phase  <t>  is  the  error  of  the  fringe  servo;  fringe 
position  is  actually  given  by 

jt(0  =  /-(()  + <^(^)/(i,  (36) 

where  L  is  the  position  of  the  laser-monitored  delay 
line  and  k  =  2r/\. 

The  two-color  phase  difference  0.^  is  02  — 
where  02  is  the  fringe  phase  in  the  blue  channel.  As 
discussed  in  the  Appendix,  the  coherence  time  of  this 
signal  can  be  several  seconds.  A  two-color  coherent 
integration  time  of  100  msec,  comprising  25  4-msec 
one-color  coherent  integration  times,  was  used  in  these 
reductions.  Such  long  integration  times  allow  the  ap¬ 
plication  of  two-color  techniques  to  cool  stars  which 
provide  relatively  small  photon  fluxes  in  the  blue  chan¬ 
nel.  The  quantities  accumulated  each  4  msec  are  the 
quadrature  components  for  the  two-color  phase  differ¬ 
ence 

.\  =  >' .  sinl(\|/X.)<;>il  +  Xj  cos((X,A.,)<i||,  (37al 

>  =  -  X,sin|(X,/X,,)((.,|,  (37b) 

where  X.>  and  are  the  quadrature  components  for 
the  blue  channel.  The  error  due  to  the  mismatch 
between  the  effective  blue  wavelength  and  the  length 
of  the  path -length  modulation  stroke  tends  to  average 
to  zero  when  these  quantities  are  accumulated  due  to 
the  randomizing  effect  of  the  fluctuations  in  0t(t). 
The  estimated  phase  in  the  red  chimnel  is  incoherently 
averaged  during  this  time  to  yield  0i.  Thus,  every  100 
msec,  the  average  two-color  phase  difference  is  esti¬ 
mated  as 
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(38) 

and  the  two-color  fringe  position  is  given  as 

ijc  «  L  +  ij/fci  -  D^2/*2-  *39) 

From  a  preliminary  first  pass  through  the  data  using 
only  a  coarse  sampling  of  the  laser  positions,  an  ap- 
proaimate  base  line  is  determined  and  used  to  subtract 
the  diurnal  component  s(t)  from  X2c(t)  to  yield  the 
discriminant  ivAt)-  In  the  absence  of  atmospheric 
and  instrumental  instabilities,  Xacit)  should  be  exactly 
zero.  However,  the  two-color  technique  is  sensitive  to 
errors  in  the  measurement  of  the  fringe  phases,  and 
some  degree  of  data  qualification  is  necessary  to  reject 
anomalous  data  to  prevent  the  introduction  of  biases 
into  the  two-color  position.  Additional  qualification 
is  necessary  to  deal  with  the  problems  of  fringe  hop¬ 
ping  and  phase  unwrapping. 

Implicit  in  the  definition  of  fringe  phase  is  proper 
reference  to  the  central  white  light  fringe.  While  the 
true  fringe  phase  is  an  unambigous  quantity,  equal  to 
zero  at  the  midpoint  of  the  central  fringe,  this  phase 
must  be  reconstnicted  from  the  estimates  provided  by 
Eq.  (34),  which  are  necessarily  modulo  2t.  Due  to  the 
finite  tracker  bandwidth,  excursions  of  can  exceed 
hx,  and  thus  some  sort  of  phase  unwrapping  is  neces¬ 
sary  to  maintain  an  analytic  estimate.  In  addition, 
there  is  an  ambiguity  when  initializing  the  unwrapper 
memory.  A  simple  first-difierence  phase  unwrapper 
for  01  was  employed  in  the  data  reduction,  and  refer¬ 
ences  to  01  above  actually  refer  to  the  unwrapped 
quantity. 

To  simplify  the  data  processing  and  minimize  the 
propagation  of  unwrapping  errors,  the  unwrapper 
memory  is  reset  at  the  beginning  of  each  lOO-msec 
integration.  This  is  equivalent  to  assuming  that  the 
initial  value  of  0i  is  within  of  true  zero  phase.  If 
this  is  not  the  case,  due  to  an  undetected  unwrapping 
error  by  the  real-time  system  which  caused  a  second¬ 
ary  fringe  to  be  temporarily  tracked,  or  just  due  to  a 
la^e  excursion  from  the  proper  track  point,  the  aver¬ 
age  phase  01  will  be  in  error  by  2xn,  where  n  is  the 
number  of  fringes  the  initial  value  of  0i  was  removed 
from  the  central  fringe.  The  two-color  phase  differ¬ 
ence  02  will  then  be  in  error  by  mod(— 2irnAi /X2),  where 
the  modulo  function  is  over  the  interval  — x  to  x,  and 
the  two-color  position  iic  will  be  in  error  by 
mod(— 2xraAiA2).  Most  fringe-hop  errors  correspond 
to  n  ^  ±1  which  lead,  for  the  parameters  determined 
above,  to  two-color  errors  of  approximately  ±20  ftm. 
Most  other  errors,  such  as  phaM-unwrapping  errors 
during  the  two-color  coherent  integration  time,  do  not 
have  such  a  characteristic  signature  and  are  more  diffi¬ 
cult  to  reject. 

Data  selection  was  performed  every  1(X)  msec.  For 
the  quantity  if  v( ,  windows  of  approximately  A 1  - 

A2)  «  ±10  Min  were  used  to  detect  most  of  t  he  fringe 
hop  errors  described  above.  While  these  points  could, 
in  principle,  be  edited  for  inclusion  into  the  data  set, 
their  frequency  was  such  that  it  was  simpler  to  just 
reject  them.  Data  were  also  selected  on  the  basis  of 


the  fringe-amplitude  SNR  in  the  red  and  blue  chan¬ 
nels.  This  quantity  is  estimated  as 

/■  =  +  y*).  (40) 

where  ^  is  the  estimated  photon  fl(u  for  that  scan. 
Rejecting  those  data  points  with  very  low  SNRs  elimi¬ 
nates  a  number  of  anomalous  data  points.  After  per¬ 
forming  a  second  pass  through  the  data  using  some¬ 
what  tighter  windows  for  X2C>  the  data  were  averaged  in 
groups  of  10  to  provide  nominal  l-sec  estimates  of  X2c- 
A  modified  median  filter  was  employed  which  aver¬ 
aged  the  median  eight  out  of  ten  lOO-msec  points, 
selected  again  on  the  basis  of  i2c-  Due  to  data  selec¬ 
tion,  the  lOO-msec  points  were  not  necessarily  continu¬ 
ous,  and  times  were  also  averaged  with  the  fringe  posi¬ 
tions  to  provide  a  mean  time  for  each  nominal  l-sec 
data  point.  The  maximum  time  interval  over  which 
this  incoherent  averaging  was  performed  was  limited 
to  10  sec  to  eliminate  sm(x>thing  errors  prior  to  sub¬ 
traction  of  the  sidereal  rate.  For  comparison,  the  one- 
color  fringe  position  was  estimated  simultaneously 
with  the  two-color  position,  using  the  same  qualified 
data  and  the  same  averaging  scheme. 

Figures  1  and  2  show  the  one-  and  two-color  fringe 
positions,  nominally  1 -sec/data  point,  after  subtrac¬ 
tion  of  the  sidereal  rate.  Figures  3  and  4  show  this 
same  data  after  simple  averaging  in  ten-point  groups 
to  provide  a  resolution  of  nominally  10  sec/point. 
From  these  figures,  it  is  clear  that  the  use  of  the  two- 
color  method  has  significantly  reduced  the  intermedi¬ 
ate-frequency  noise,  presumably  mostly  atmospheric 
in  origin.  However,  a  low-frequency  residual  is  clearly 
evident.  The  presence  of  this  residual  in  both  the  one- 
and  two-color  plots  suggests  that  it  is  a  systematic 
instrumental  drift.  (The  slight  difference  between 
the  residuals  at  the  ends  of  the  one-  and  two-color  data 
is  due  to  systematic  residual  astronomical  refraction  in 
the  one-color  data,  which  is  not  present  in  the  two- 
color  position.) 

The  relative  performance  of  the  one-  and  two-color 
estima'  es  becomes  much  clearer  in  the  frequency  do¬ 
main.  Power  spectra  were  obtained  using  the  l-sec 
averaged  data,  interpolated  where  necessary  to  ap¬ 
proximate  uniform  sampling.  These  spectra  are 
shown  in  Fig.  5.  All  the  data  were  used  to  estimate  the 
spectra  at  the  very  lowest  frequencies,  while  only  the 
first  half  of  the  data,  where  the  fraction  of  interpolated 
points  was  smallest,  was  used  at  intermediate  frequen¬ 
cies.  The  power  spectra  are  terminated  at  a  frequency 
below  that  at  which  interpolation  artifacts  become 
significant.  Systematic  instrumental  effects  are  visi¬ 
ble  in  both  the  one-  and  two-color  spectra  at  frequen¬ 
cies  below  ~10^  ’  Hz.  These  errors  typically  have  a 
power  spectrum  obeying  power  laws  between  f~'  and 
A '  and  are  likely  due  to  thermal  instabilities. 

At  frequencies  above  ~10~ '  Hz,  the  power  spectra 
are  mostly  atmospheric  in  nature.  The  one-color 
spectrum  obeys  the  expected  /'  power  law,  whose 
asymptote  is  indicated  in  the  figure;  the  one  color 
spectral  height  corresponds  to  a  turbulence-liiuiied 
absolute  astrometric  accuracy  ol  ~Go  msec  of  at.  )<  r  a 
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Pig.  1.  One-color  fringe  position  for  ,i  Cas.  25  Nov.  1986,  at  3  One-color  fringe  position  for  (1  Cas,  25  Nov.  1986.  at 

10  sec/point. 
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Fig.  2.  Two-color  fringe  position  for  8  Cas,  25  Nov.  1986.  at  Two-color  fringe  position  for  i?  Cas,  25  Nov.  1986,  at 

1  sec/point.  10  sec/point. 


lOO-sec  integration.^^  Over  most  of  the  range  above 
10~^  Hz,  the  two-color  spectrum  shows  about  a  15-dB 
improvement  over  the  one-color  spectrum  or  a  factor 
of  reduction  in  error  amplitude.  The  bump  in  the 
two-color  spectrum  near  0.0015  Hz  is  believed  to  be  a 
siderostat  drive  error  which  is  masked  by  atmospheric 
noise  in  the  one-color  spectrum.  This  error  should  be 
eliminated  with  additional  laser  metrology  on  the  in¬ 
terferometer  to  be  installed  in  the  near  future. 

Thus,  despite  the  limitations  of  the  two-color  meth¬ 
od,  which  were  discussed  in  the  previous  sections, 
these  data  indicate  that  a  two-color  absolute  astrome¬ 
tric  accuracy  of  ~13  msec  of  arc  should  be  obtainable 
from  a  single  lOO-sec  measurement.  For  star- 
switched  relative  measurements,  this  error  should  de¬ 
crease  as  the  square  root  of  the  number  of  switching 
cycles,  so  that  wide-field  relative  accuracies  of  sever^ 
msec  of  arc  should  be  obtainable  in  a  single  night. 
While  the  effect  of  systematic  errors  and  the  appiica- 
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Fig.  5.  Power  apectnun  of  one-  and  two-color  fringe  positions  for 
rf  fas.  25  Nov.  1986. 
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tion  of  the  two-color  method  to  a  large  data  set  need  to 
be  explored  experimentally,  these  preliminary  results 
sugg^  that  a  significant  improvement  in  the  accuracy 
of  ground-based  astrometric  measurements  is  possible 
wiU)  two-color  techniques. 

Appendbr:  Coheranoe  of  the  Two-Color  Phase 

Consider  the  two-color  phase  difference 

02(t)  -  02  -  (X,A,)0„  (Al) 

where  and  ^2  are  the  fringe  phases  at  wavelengths  Xi 
and  X2.  Because  of  rorrelations  between  and  4>2,  the 
coherence  time  of  02  be  much  greater  than  that  of 
either  0i  or  02.  In  the  geometric-optics  limit  for  obser¬ 
vation  of  stars  at  zenith,  02  is  equid  to  (l/D)(Xi/X2)0i, 
where  D  is  the  dispersion.  As  02  is  primarily  a  measure 
of  the  one-color  atmosphere,  water  vapor  can  be  ig¬ 
nored  in  this  calculation.  The  atmospheric  coherence 
time  of  the  one-color  phase  0i  can  be  written^^'^^ 

Hence  the  coherence  time  Tg  for  the  two-color  phase 
difference  02  can  be  written 

To  «  (DXj/X, ,  (A3) 

Using  a  Hufnagel  turbulence  profile  with  V(h)  =  V  * 
10  m/sec  yields  To  *  13  msec  at  Xi  *  0.7  ^m.  Thus, 
with  X2  «  0.5  >im,  which  gives  D  *  87.1,  the  two-color 
phase  difference  coherence  time  Tq  is  equal  to  1.8  sec. 
With  a  coherent  intention  time  for  the  two-color 
phase  difference  of  T^,  the  variance  of  the  phase  fluc¬ 
tuations  about  the  interval  mean  can  be  written  <rf  » 
these  fluctuations  cause  a  fractional  visibil¬ 
ity  reduction  of  approximately  <tV2  for  «  1.  Thus, 
even  with  Tj  *  100  msec,  the  reduction  in  fringe  visi¬ 
bility  is  <1%. 

As  discussed  in  Secs.  Ill  and  IV,  the  effects  of  dif¬ 
fraction  and  refraction  reduce  the  correlation  between 
01  and  02.  Consider  the  effects  of  diffraction  for  ob¬ 
servation  of  a  star  at  zenith.  Assuming  a  low  fringe- 
tracker  bandwidth,  the  process  02  as  given  above  will 
be  proportional  to  d(t)  of  Eq.  (13).  Therefore,  the 
power  spectrum  of  02  will  be  proportional  to  Waaif), 
Eq.  (25).  Using  the  notation  of  that  section,  the  vari¬ 
ance  of  02  due  to  diffraction  can  be  written 

aj  -  ().50(X|/X,)-f(*,,*2)ff, (A4) 

where  is  the  variance  of  the  log  intensity.'  The 
variance  (A4)  is  an  upper  bound  on  the  fluctuations 
about  the  interval  mean  which  would  occur  in  a  Hnite 
coherent  integration  time.  With  a  Hufnagel  turbu¬ 
lence  profile,  the  log  intensity  variance  at  0.7  ^m  is 
given  by  <r'f„i  »  0.2.  Thus,  with  F{ki,k2)  ■  0.14,  the 
variance  (A4)  is  given  by  »  0.03  rad-,  which  would 
lead  to  a  maximum  fractional  visibility  reduction  of 
only  1.5%.  Thus  diffractive  effects  can  be  safely  ig¬ 
nored. 

Consider  now  the  effect  of  refraction,  computed,  as 
in  Sec.  IV,  in  the  geometric-optics  limit.  Assuming 
independent  phase  errors  at  each  interferometer  aper¬ 


ture,  the  variance  of  02  attributable  to  beam  separation 
for  a  finite  zenith  angle  is  given  by  twice  the  phase 
structure  function'  evaluated  at  A(h),  i.e., 

4  =  2  X  2.92klf;dhCl(h)A^'Hh),  (A5) 

where  A(h)  is  given  by  Eq.  (28)  of  Sec.  IV.  Substitut¬ 
ing  Eq.  (28)  into  Eq.  (A5)  and  incorporating  an  air 
mass  correction  yield 

4  =  5.84*.1|A(/»((X,)  eip(-A,/A„)/Z))*'' 

X  (1  -  exp(-A/Ao)l*".  (A6) 

Using  a  Hufnagel  turbulence  proflle,  the  final  integral 
has  a  value  of  1.5  X  10~'^  m'^^.  Using  the  parameters 
in  Sec.  IV,  the  variance  v|'  is  equal  to  0.03,  0.13,  and 
0.57  rad-  for  ff  =  15,  30,  and  45®,  which  yields  upper 
bounds  to  the  fractional  visibility  reductions  of  ~2,  7, 
and  28%,  respectively.  These  results  were  calculated 
assuming  geometric  optics,  point  apertures,  and  a  zero 
inner  scale  for  the  turbulence.  However,  even  at  45® 
zenith  angle,  the  maximum  beam  separation  is  only  ~3 
cm;  the  Fresnel  zone  size,  aperture  diameter,  and  inner 
scale  are  all  of  this  order.  Therefore,  these  results 
should  be  fairly  conservative,  and  refraction  should 
not  be  a  major  limitation  to  long  coherent  integration 
times  of  the  two-color  phase  difference. 
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ABSTRACT 

Preliminary  astrometric  measurements  with  the  Mark  III  optical  stellar  interferometer  using  a  12,005 
m  north-south  baseline  were  made  during  seven  nights  in  November  1986.  A  total  of  22  stars  were 
observed  over  large  angles  in  a  3000-A  wide  band  centered  at  7000  A.  As  expected,  the  estimated 
precision  of  the  relative  positions  increases  with  hour-angle  coverage.  For  the  ten  surs  observed  over  at 
least  90*  of  hour  angle,  the  right  ascensions  are  consistent  at  the  S  ms  of  time  (]a)  level,  while  the 
uncertainties  in  declination  are  at  the  22  mas  (lo)  level.  With  this  small  data  set,  the  question  of 
systematic  errors  cannot  be  fully  examined.  A  discussion  of  observational  techniques,  errors,  limiutions 
of  this  system,  and  future  plans  is  presented. 


1.  INTRODUCTION 

Determining  accurate  positions  of  surs  over  wide  angles 
allows  imporunt  measurements  not  otherwise  possible.  For 
example,  fundamental  constants,  such  as  those  of  preces¬ 
sion,  nuution,  and  galactic  rotation  require  wide-angle  as¬ 
trometry.  Measurements  at  optical  wavelengths  are  being 
supplanted  by  those  made  in  the  radio  bands.  Positions  of 
extragalactic  radio  sources  accurate  to  a  few  milliarcseconds 
are  now  available  over  the  northern  hemisphere  from  radio 
observations.  A  similar  precision  is  needed  in  the  visible,  if 
we  hope  to  independently  align  radio  maps  with  optical  im¬ 
ages  made  with  such  systems  as  the  Hubble  Space  Telescope, 
which  will  have  a  resolution  of  80  mas. 

Gassical  astrometry  uses  the  meridian  transit  circle  to 
determine  the  position  of  stars  over  large  angles  to  accura¬ 
cies  of  order  0.1  arcsec.  The  new  FK5  optical  reference 
frame  is  believed  to  be  internally  consistent  at  the  0.05  arcsec 
level  at  epoch  1987.  The  alignment  of  this  reference  frame 
with  respect  to  the  radio  reference  frame  has  yet  to  be  deter¬ 
mined.  We  believe  that  the  radio  reference  frame  is  probably 
internally  consistent  at  the  level  of  a  few  milliarcseconds 
(Johnston  and  DeVegt  1986;  Ma  etal.  1986). 

Modem  techniques  are  evolving  that  promise  improved 
accuracy  in  determining  sur  positions.  Space  techniques 
such  as  the  HlPPARCOStkXtV&le  promise  a  reference  frame 
of  stars  down  to  magnitude  10  with  accuracies  of  a  few  mil- 
liarcseconds.  As  space  techniques  are  very  costly  and  de¬ 
ployment  in  orbit  is  uncertain  at  this  time,  the  development 
of  ground-based  techniques  is  very  important.  Initial  mea¬ 
surements  with  an  interferometer  at  Mount  Wilson  Obser¬ 
vatory  gave  stellar  positions  accurate  to  3  arcsec  (Shao  et  al. 
1987).  A  new-generation  instrument,  the  Mark  III  stellar 
interferometer,  it  in  the  process  of  fabrication.  Here  we  re- 
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port  preliminary  measurements  with  that  instrument  using 
the  initial  12  m  north-south  baseline. 

Further,  interferometry  appears  to  be  a  technology  that 
offers  improved  accuracy  for  both  ground  and  space  obser¬ 
vations.  Before  space-based  observations  are  planned  in  de¬ 
tail  using  this  technique,  experience  must  be  gained  through 
the  development  of  ground-based  systems. 

II  THE  INSTRUMENT 

The  Mark  III  interferometer  is  the  third  in  a  series  of  opti¬ 
cal  interferometers  built  on  Mount  Wilson,  California.  Stel¬ 
lar  fringes  were  first  tracked  with  this  instrument  in  Septem¬ 
ber  1986,  By  November,  the  system  was  operating  well 
enough  to  allow  stellar  observations  to  continue  all  night. 
This  instrument  was  specifically  built  to  do  fundamental  as¬ 
trometry.  Several  features  were  incorporated  into  its  design 
to  allow  it  to  accomplish  this  goal.  The  fringes  were  tracked 
in  a  closed-loop  mode,  to  allow  accurate  measurement  of  the 
delay.  Sur  tracking  and  fringe  tracking  were  completely 
automated,  allowing  rapid  switching  between  stars.  One 
hundred  sixty  observations  were  accomplished  during  one 
night.  The  light  from  the  siderostats  was  mixed  in  a  central, 
temperature-controlled  building,  and  the  siderostats  were 
mounted  on  massive  concrete  piers.  As  a  result,  the  baseline 
is  fairly  suble  at  the  micron  level.  Other  features  designed 
and  built,  but  not  yet  operational,  are  expected  to  further 
improve  the  performance  of  the  system.  We  will  limit  our 
description  of  the  instrument  to  the  points  relevant  to  the 
astrometric  performance.  A  full  description  of  the  instru¬ 
ment  is  given  in  Shao  ei  al.  ( 1988 ) . 

The  interferometer  consists  of  two  2S-cm-diameter  sider- 
osuts,  mounted  on  massive,  insulated  concrete  piers  and 
currently  operating  with  a  north-south  baseline  of  12.005  m 
Four  additional  piers  have  been  constructed,  allowing  future 
observations  with  a  20  m  north-south  baseline,  as  well  as 
with  two  oblique  baselines.  From  the  siderosuts.  the  light 
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puses  through  evKuated  tubu  into  a  central,  temperature- 
controlled  building.  The  light  from  each  siderosut  is  direct¬ 
ed  into  delay  lines  containing  nsovable  retroreflectors.  Next, 
the  light  is  combined  at  a  beam  combiner  and  the  multing 
fringes  are  detected  with  a  photomultiplier  tube  using  the 
technique  described  below.  Fw  the  observations  discussed 
here,  tte  clear  aperture  wu  len  than  2  in.  and  light  from 
only  one  side  of  the  beam  combiner  wu  detected.  To  main¬ 
tain  parallel  wavefronu  at  the  beam  combiner,  li^t  from  an 
annulu  surrounding  the  central  aperture  wu  deflected  into 
a  star  tracker  which  adjusts  the  anglu  of  two  flat  mirrors. 
All  of  the  optical  hardware  in  the  central  building  is  on  a 
supeiinvv  table  except  the  delay  linu.  However,  the  posi- 
tiou  of  the  retroreflectors  in  the  delay  lino  are  meuured 
relative  to  the  superinvar  table  with  a  laser  metrology  sys¬ 
tem. 

To  detect  fringo,  the  length  of  the  delay  line  is  dithered  in 
a  SOO  Hz  triangle  wave  with  an  amplitude  of  one  optical 
wavelength.  The  ‘stroke*  is  divided  into  four  equal  s^ments. 
The  number  of  photons  observed  in  each  of  the  segments  is 
used  to  determine  the  position  of  the  fringe  for  real-time 
tracking.  A  vride-bandpus  filter  wu  used  to  increase  the 
sensitivity  and  to  decrease  the  contrut  of  all  but  the  central 
fringe.  Four  milliseconds  of  data  are  coherently  combined. 
This  time  constant  wu  chosen  to  be  short  relative  to  the  time 
constant  of  the  atmosphere,  so  that  atmospheric  turbulence 
is  frozen  during  a  meuurement. 

III.  DATA 

The  nightly  observing  procedure  consisted  of  cycling 
through  a  list  of  stars.  Each  observation  consisted  of  track¬ 
ing  the  central  fringe  on  one  star  for  100  s  aAer  the  fringe 
position  wu  found.  During  periods  of  poor  seeing  or  at  large 
zmth  angles,  the  fringe  tracker  would  occuionally  lose  the 
fringe  and  take  several  seconds  to  reacquire  it.  During  good 
seeing,  the  fringe  tracker  would  successfully  track  a  fringe 


for  several  minutes.  As  many  u  160  observations  were  possi¬ 
ble  during  a  single  night. 

In  reducing  the  data,  an  observation  wu  first  averaged  to 
1  s.  These  1  s  delays  were  averaged  by  differencing  the  ob¬ 
served  and  predicted  delays  and  adding  the  mean  difference 
to  the  predicted  delay  at  the  mean  time.  If  an  observation 
consisted  of  fewer  than  ten  1  s  averages,  it  wu  discarded. 
The  data  were  not  edited  in  any  other  way.  A  plot  of  the 
observed  delays  u  a  funetkm  of  time  is  shown  for  a  typical 
night  in  Fig.  1.  The  large  dots  are  the  100  s  averages.  Multi¬ 
ple  observations  of  the  same  sUr  are  connected  with  small 
dots. 

All  of  the  utrometric  information  is  contained  in  the  ob¬ 
served  delay.  It  is  related  to  the  star’s  position  through 

d  =  Co-)-B-S,  (1) 

where  B  is  the  baseline  vector,  S  is  a  unit  vector  pointing  in 
the  unrefracted  direction  of  the  star,  and  Q  is  a  constant 
which  depends  on  the  length  of  the  optical  path  inside  the 
central  building  and  on  the  differential  atmospheric  path 
length  between  the  receiving  apertures.  Note  that  for  a  hori¬ 
zontal  baseline  with  a  vacuum  dday  fine;  the  air  path  is  near¬ 
ly  identical  for  both  arms  of  the  interferometer,  causing  the 
effects  of  atmospheric  refraction  to  cancel  out  (see  the  Ap¬ 
pendix).  We  adopt  the  convention  of  meuuring  the  delay  u 
a  length.  Expanding  this  in  a  right-handed  coordinate  sys¬ 
tem.  with  Z  axis  parallel  to  the  Earth's  axis,  and  F  axis  to¬ 
ward  the  eut  point  on  the  horizon,  gives 

Cf)  +  sin(S)  —  cos(£)sin(/t) 

+  B,  co&{6)coi(h),  (2) 

where  A  =  /  -  a  is  the  instantaneous  hour  angle,  f  is  the  local 
sidereal  time,  and  S  is  the  declination.  As  the  coordinate 
frame  depends  upon  the  roution  of  the  Earth,  the  declina¬ 
tions  are  with  respect  to  the  instantaneous  celestial  pole 
while  the  right  ascension  zero  point  is  arbitrary.  Observa- 
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Fic.  1.  Observed  delay  (in  mi¬ 
crons)  venus  uaiversal  time  (se¬ 
conds)  on  5  November  1986  Note 
that  each  targe  dot  is  the  meanofa 
too  s  integration.  Small  dou  con¬ 
nect  observations  of  the  same  sur. 
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lions  of  the  delay  of  a  star  as  a  function  of  time  provide  three 
independent  measurements,  the  amplitude  and  phase  of  a 
sinttsoid  and  a  constant  dfset.  A  zero  point  for  right  ascen¬ 
sion  cannot  be  determined  fhmi  these  observations  since 
then  is  no  tie  to  either  the  ecliptic  or  to  a  well-dellned  local 
longHude  reference.  Observation  of  N  stars  at  various  decli¬ 
nations  can  be  described  by  specifying  2Ar  -t-  3  parameters, 
the  three  components  of  the  baseline,  the  dday  constant,  N 
declinations,  and  S—  1  right  ascensions.  Wi^  three  inde¬ 
pendent  measurements  each,  at  least  three  start  must  be  ob¬ 
served  to  determine  all  of  the  parameters.  In  practice,  this 
solution  is  not  well  constrained.  Expanding  Eq.  (2)in  terms 
of  assumed  baseline  components  and  star  position  plus  small 
corrections, 

md  —  Co  —  B,  sin(d)  +  B,  cos(d)sin(A) 

—  B„  cos(5)cos(A) 

»=  ACo  -f  A^,  sin(5)  +  ^B,  cos(d) 

-i-  sin(A)(  —  AHycosCd)  -t-  &SBysm(S) 

-1-  Aa.0x  cos(d)  ] 

-f  cos(A)(AArcos(^)  +  ^Bycm{S) 
-ASBxm(S)].  (3) 

where  S  and  A  are  now  the  assumed  declination  and  hour 


angle.  Since  ^x~6.8  m  and  Bz-~lOm  and  By^i  cm,  the 
residuals  are  dominated  by  the  B^  and  B,  terms.  If  we  want 
to  determine  equally  good  values  for  Ad  and  Aa,  the  coeffi¬ 
cients  of  sin(  A)  andcos(A)  must  be  determined  to  the  same 
precision.  With  6  hr  of  coverage,  centered  on  transit,  cos(A ) 
changes  by  only  ±0.15,  whereas  sin(A)  varies  by  ±  0.7. 
The  cos(A)  term  can  be  separated  from  the  constant  term 
with  only  about  one-fifth  of  the  precision  of  the  other  terms. 
If  the  baseline  or  the  delay  varies  slowly  with  time,  more 
parameters  must  be  added  to  the  lit,  further  hindering  their 
determination. 

We  used  a  two-step  procedure  to  determine  star  positions 
from  the  observed  delay.  First,  we  solved  for  the  baseline 
using  a  priori  estimates  of  the  star  positions.  Then  we  used 
the  residual  delays  from  this  fit  to  determine  improved  sur 
positions.  Since  we  used  a  priori  star  positions  to  determine 
the  baseline,  the  improved  star  positions  are  not  absolute, 
but  are  relative  to  a  mean  position  of  the  stars  in  the  fit. 

For  the  baseline  fit,  we  used  star  positions  taken  from  the 
FKS  catalog,  transformed  to  apparent  topocentric  place  of 
date  in  accordance  with  recent  lAU  resolutions  (see  Kaplan 
1981 ).  For  the  two  binaries  in  our  sample,  the  positions  in 
the  catalog  correspond  to  the  center  of  mass  of  the  systems. 
We  used  the  orbit^  positions  determined  by  Worley  ( private 
communication )  to  determine  the  positions  of  the  primaries. 
The  parameters  used  are  shown  in  Table  I.  In  addition,  the 
differential  optical  path  length  through  the  atmosphere  was 


Tasle  I.  Astrometric  peremeten  for  the  observed  stars 


Name 

RA 

DEC 

■g 

■■ 

nmn 

RV 

IBBi 

a 

And 

00 .08  .-23. 265 

29:05:25.58 

1.039 

-16.33 

0.0240 

-12 

fi 

Cas 

00:09:10.695 

59:08:59.18 

6.827 

-18.09 

0.0720 

11 

a 

Cas 

00:40:30.450 

56:32:14.46 

0.636 

-3.19 

0.0160 

-4 

Ts 

And 

02:03:53.963 

42:19:46.99 

0.404 

-5.17 

0.0130 

-12 

X 

O 

Arl 

02:07:10.403 

23:27:44.66 

1.383 

-14.83 

0.0430 

-14 

Q 

Per 

03:24:19.365 

49:51:40.34 

0.246 

-2.46 

0.0290 

-2 

7 

Ori 

05:25:07.857 

6:20:58.74 

-0.059 

-1.39 

0.0260 

18 

B 

Tau 

05:26:17.511 

28:36:26.67 

0.169 

-17.51 

0.0180 

9 

P 

Aur 

05:59:31.720 

44:56:50.78 

-0.541 

0.03 

0.0370 

-18 

y 

Gea 

06:37:42.725 

16:23:57.43 

0.293 

-4.16 

0.0310 

-13 

a 

Gen  a 

07:34:35.997 

31:53:18.53 

-1.347 

-9.87 

0.0720 

-1 

a 

CMi« 

07:39:18.113 

5:13:30.06 

-4.755 

-102.29 

0.2880 

-3 

P 

Gea 

07:45:18.946 

28:01:34.26 

-4.740 

-4.59 

0.0930 

3 

a 

Leo 

10:08:22.315 

11:58:01.89 

-1.693 

0.64 

0.0390 

6 

P 

UMa 

11:01:50.482 

56:22:56.65 

0.988 

3.4 

0.0420 

-12 

a 

UMa 

11:03:43.666 

61:45:03.22 

-1.675 

-6.65 

0.0310 

-9 

y 

Dra 

17:56:36.367 

51:29:20.21 

-0.081 

-1.94 

0.0170 

-28 

a 

Lyr 

18:36:56.332 

38:47:01.17 

1.726 

28.61 

0.1230 

-14 

a 

Aql 

19:50:47.002 

8:52:06.03 

3.629 

38.63 

0.1981 

-26.3 

y 

Cyg 

20:22:13.704 

40:15:24.14 

0.03 

0.0030 

-8 

a 

Cyg 

20:41:25.917 

45:16:49.31 

Bin 

0.23 

0. 

-5 

c 

Per 

21:44:11.164 

9:52:29.92 

0.207 

-0.06 

0.0060 

5 

Position  is  for  ths  cantor  of  mss  of  the  binary  systuB.  To  find  the 
position  of  the  observed  star  at  the  tlae  of  the  observations  add  to  a  Gem 
Aa— 0.  0793  and  At— 0. ”158  and  to  a  CMl  Aa— 0.*008S  and  A8— 1.*378 
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ettimated  and  removed  frun  the  data.  A  pivot-point  correc¬ 
tion  to  account  fw  the  siderosut  geometry  was  tested  but  not 
included,  since  it  had  no  significant  effect  on  the  results.  The 
form  and  magnitude  of  th^  corrections  are  discussed  in  the 
Appendix. 

^ectively,  then,  the  baseline  was  determined  relative  to 
the  FKS  system.  The  baseline  was  not  constant,  but  drifted 
slowly  during  the  observations.  On  four  nights  during  the 
fall  of  1986,  a  laser  interterometer  was  operated  in  a  vacuum 
tube  between  one  of  the  unused  piers  a^  the  central  build¬ 
ing,  although  not  connected  to  the  superinvar  table.  These 
dau  were  taken  before  some  of  the  final  improvements  were 
made  to  the  system  and  may  not  be  relevut  to  the  stellar 
observations.  Nevertheless,  there  were  large  variations  in  the 
path  length  for  the  first  few  hours  after  sunset.  Later  in  the 
night,  the  length  drifted  by  about  1  pm  per  hour.  A  constant 
drift  rate  provides  an  acceptable  fit  (less  than  1  pm  residu¬ 
als)  to  the  data  later  in  the  night.  The  large  variations  in  the 
early  evening  are  almost  certainly  due  to  large  temperature 
changes  that  occur  at  sunset  and  should  take  place  prior  to 
the  start  of  our  observations  in  November.  Therefore,  we 
assumed  that  the  delay  constant  and  the  componenu  of  the 
baseline  varied  linearly  during  the  stellar  observations.  In 
Eq.  (2),  Co  was  replaced  with  Q  +  C,  t,  with 
+  ^1.  A  etc. 

There  are  two  major  sources  of  noise:  variations  in  the 
baseline  and  atmospheric  turbulence.  ( 1 )  Noise  due  to  base¬ 
line  variations  is  caused  by  effects  such  as  thermal  drift, 
bearing  errors,  and  flexure,  and  will  occur  over  periods  of 
hours  or  when  comparing  observations  taken  in  different 
parts  of  the  sky.  (2)  The  atmospheric  effects  include  both 
fluctuations  ateut  the  mean  delay  and  a  systematic  offset. 
The  small  offset  occurs  because  the  path  length  through  the 
atmosphere  is  slightly  different  for  the  two  siderostats.  With 
a  north-south  baseline,  this  offxt  is  mainly  a  function  of 
declination.  It  has  been  estimated  and  removed  using  the 
procedure  outlined  in  the  Appendix.  The  fluctuations  have 
an  amplitude  of  as  much  as  15  pm  at  a  frequency  of  a  few 
tenths  of  a  Hertz.  Using  data  taken  with  this  interferometer, 
Colavita  et  al.  ( 1987)  have  shown  that  the  power  spectrum 
of  these  fluctuations  goes  Since  this  is  close  to  1// 

noise,  we  assumed  that  the  uncertainty  of  a  100  s  average  due 
to  atmospheric  turbulence  was  the  rms  width  of  the  distribu¬ 
tion  of  1  s  averages  about  the  mean  value.  This  is  not  strictly 
correct,  since  the  noise  should  integrate  down  as  the  sixth 
root  of  the  integration  time.  Thus  tlw  noise  estimate  used  is 
roughly  a  factor  of  2  larger  than  the  noise  expected  from  the 
atmosphere  alone.  A  histogram  of  these  values  is  presented 
in  Fig.  2. 

The  baseline  parameters  were  determined  separately  for 
each  night’s  dau  using  standard  linear-regression  analysis. 
In  wei^ting  the  data,  we  assumed  that  the  atmosphere  con¬ 
tributed  all  of  the  noise.  A  plot  of  the  residual  delays  as  a 
function  of  time  is  shown  in  Fig.  3.  The  derived  parameters 
were  collected  in  Table  II.  The  observing  times,  number  of 
observations,  mean  uncertainties  of  the  observation,  and 
the  rms  and  of  the  fits  are  collected  in  Table  III.  The 
standard  deviation  of  the  residual  delays  varies  from  4.6pm 
on  14  November  to  7.7  pm  on  1 1  November.  When  data 
from  all  seven  nights  are  combined,  the  standard  deviation  of 
the  residuals  is  6.5  pm.  The  days  with  the  largest  systematic 
variations  in  the  residuals  ( 5  November,  10  November)  are 
also  days  with  large  values  of  the  reduced  On  four  of  the 
nights,  the  reduced  is  less  than  1.2,  and  the  rms  of  the 
residuals  after  the  fit  is  leu  than  the  mean  uncertainty  of  the 


HISTOGRAM  OF  ESTIMATED  UNCERTAINTIES 
140 


120 


100 


ao 


so 


40 


20 


0 

0  24aaioi2i4i«ia20 
MICRONS 

Fic.  2.  Hittotmn  of  the  tuumed  unoerttinty  (in  mi- 
croni)  of  the  obseivuions  of  all  data  (rom  November 
1986.  The  obtervaiiont  arc  I00>  in  IcttSth,  short  enough 
to  be  dominated  by  atmospheric  fluctuations.  The  dau 
were  averaged  to  I  s  before  determining  the  rms.  A 
change  in  declinationof  I  arcsec  corresponds  to  a  change 
in  delay  of  roughly  SS  ^m. 


lOOOO  20000  SOOOO  40000  60000 


UT  (seconds) 

Fic.  3.  Residual  delay  versus  time,  (a)  5  November  1986:  (b) 
6  November  1986.  Note  that  the  observationt  obtained  on  S 
November  have  systematic  eflccu  of  S-lOpm  peak-io-peak. 


-20782.2*M.4 
•207M.0*  S.6 
-207W.S*28.7 
•ZOTM.O*  6.8 
•20m.4>  9.6 
-199S8.8<11.2 


675270$. lalS.S 
6752697.2<  5.2. 
6752716.6*27.9 
6752707.9*  6.8 
6752699.5*  9.4 
6752692.5*10.7 


31362.0*  5.7 
31365.3*  1.7 
31364.3*  8.5 
31371.4*  1.9 
31386.6*  2.7 
31369.8*  3.5 


•9926024.4*12.9 
•9926050.4*  4.5 
•9926052.4*20.3 
•9926082.3*  5.8 
•9926004.3*  7.4 
•9926082.4*  8.5 


4.89*  2.^ 
-1.71*  1.8 
3.39*  4.0 
-3.28*  1.8 
2.18*  2.0 
2.U*  2.5 


■ferens/hour^ 

•3.0*  2.5 

1.3*  0.9 

•5.9*  2.1 

•1.6*  1.7 

1.6*  0.6 

•1.9*  1.5 

•3.0*  3.9 

2.6*  1.2 

•5.2*  2.8 

-1.8*  1.7 

3.1*  0.5 

-2.6*  1.5 

•1.1*  1.9 

0.0*  0.6 

-3.2*  1.6 

-0.7*  2.3 

1.0*  0.8 

-1.6*  2.0 

•5.9*  4.8 

3.0*  1.9 

•11.5*  3.8 

dau  points.  On  these  nights,  the  estimate  of  the  uncertainties 
seems  to  be  valid.  However,  recall  that  it  is  about  a  factor  of  2 
larger  than  that  expected  from  the  atmosphere.  On  the  three 
remaining  nights,  the;^  is  greater  than  1 .5  and  the  rms  of  the 
residuals  is  greater  than  the  mean  of  the  uncertainties.  Two 
of  these  nights  show  large  systematic  shifts  in  a  plot  of  the 
residual  versus  time  (Fig.  3).  These  are  three  of  the  longest 
nights,  but  if  dau  from  half  the  night  is  deleted  from  the  fit, 
the  is  improved,  but  not  to  the  level  of  the  four  good 
nighu.  Thus,  on  at  least  these  three  nights,  we  seem  to  be 
limited  by  a  source  of  noise  other  than  the  atmosphere,  pre¬ 
sumably  thermal  or  mechanical  drifts  of  the  baseline. 

If  dau  from  all  of  the  stars  are  combined,  there  is  no  sys¬ 
tematic  variation  of  the  residual  delay  with  hour  angle. 
However,  plots  of  individual  surs  ( Fig.  4 )  show  such  trends. 
This  implies  that  the  fit  could  be  improved  by  solving  for 
corrections  to  the  star  positions. 

Equation  (3)  was  used  to  determine  corrections  to  the 
sur  positions,  assuming  the  fitted  baseline  was  correct.  Er¬ 
rors  in  the  baseline  estimate,  caused  by  motions  of  the  base¬ 
line  more  complicated  than  a  linear  drift,  or  by  using  incor¬ 
rect  sur  positions  in  the  previous  step,  couple  directly  into 
the  estimates  of  star  positions.  This  two-step  fitting  does  not 
allow  absolute  astrometry,  but  sur  positions  are  determined 
relative  to  the  mean  of  the  stars  included  in  the  fit.  Again,  a 
standard  linear  regression  was  used  and  the  weighting  was 
based  on  the  estimate  of  the  atmospheric  noise.  Sur  posi¬ 


tions,  relative  to  the  assumed  positions,  determined  from 
each  night’s  dau,  as  well  as  from  a  global  fit  to  all  of  the  data, 
are  collected  in  Tables  IV  and  V.  The  stars  are  arranged  in 
order  of  increasing  right  ascension.  Stars  at  the  beginning 
and  end  of  the  list  were  only  observed  near  setting,  or  rising, 
respectively,  and  on  some  nights  have  very  few  observations. 
If  less  than  three  observations  were  made  of  a  star  during  a 
night,  a  fit  position  for  that  star  was  not  included  in  the 
Ubles,  although  those  observations  were  used  for  fitting  the 
baseline  and  for  the  fit  to  all  seven  nights'  dau. 

To  determine  at  what  level  the  data  are  self-consistent,  x' 
tests  were  performed  on  the  positions  in  Tables  IV  and  V. 
The  fit  to  al]  seven  nights'  data  was  uken  as  the  mean  posi¬ 
tion.  The  positions  and  errors  determined  from  each  night 
were  uken  as  the  data  points.  There  was  one  fit  parameter 
for  each  sur  in  the  uble.  The  reduced  came  out  to  be  1 .08 
for  the  right  ascensions  and  0.97  for  the  declinations,  with  9 1 
degrees  of  freedom.  Thus,  the  data  are  consistent  with  the 
( la)  errors  ubulated.  That  is,  the  formal  errors  from  the 
least-squares  fit  are  a  good  estimate  of  the  precision  (night- 
to-night  repeaubility)  of  the  system. 

Positions  determined  for  stars  observed  only  near  rising  or 
setting  have  very  large  errors.  However,  the  errors  are 
smaller  for  surs  with  total  hour-angle  coverage  of  at  least 
90*.  These  ten  stars  have  a  mean  error  in  declination  of  22 
mas  ( la)  and  in  right  ascensions  S  ms  of  time  using  seven 
nights  of  daU.  We  stress  that  these  numbers  show  the  preci- 


Table  III  SuRimvy  of  observations. 


Date 

no.  of 
oba. 

naa 

MB 

a 

mb 

oba  tlBe 
houra 

5  Nov 

98 

6.8 

6.5 

1.64 

9.5 

8 

6  Nov 

152 

5.6 

6.9 

10.8 

5 

9  Nov 

77 

5.1 

5.8 

1.20 

7.4 

2 

10  Nov 

142 

7.3 

7.0 

1.53 

0 

11  Nov 

144 

7.7 

7.3 

1.62 

10.5 

7 

12  Nov 

94 

7.4 

9.0 

0.98 

6.4 

0 

14  Nov 

68 

4.6 

6.2 

0.84 

5.4 

0 

All 

819 

6.5 

7.1 

1.28 

61.1 

2 

RESnUAL  DELAY  (i 
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HOUR  ANGLE 


Fic.  4.  Residual  dduys  for  individual  stars  versus 
hour  angle,  (a)^  Tan,  (b)aPer,  (c)  yOri,  (dl  y 
Gem.  These  four  stars  show  the  largest  systematic 
trends  in  the  data.  For  our  observinf  geometry,  er¬ 
rors  in  declination  correspond  to  an  offset  while  er¬ 
rors  in  right  ascension  correspond  to  a  nonzero 
slope.  Note  that  the  assumed  positions  for  Tau 
and  a  Per  are  off  in  declination,  while  the  positions 
for  y  Ori  and  y  Gem  are  off  in  right  ascension. 


Table  IV.  Derived  corrections  to  right  ascensions  for  observed  stars  m  seconds  of  time. 


star 

5  Nov 

6  Nov 

9  Nov 

10  Nev 

11  Nov 

12  Nov 

14  Ncv 

mesn 

a 

And 

•.OaOOe.1714 

.0081a. 0110 

o  •  • 

■.0045a. 01 74 

■.0221O.0204 

see 

a  a  a 

•.0002O.008S 

0 

Css 

•.1298*.06t4 

. 0134a. 02U 

•  ee 

■.02080.0341 

.09370.0593 

.33180.5224 

a  a 

.00S6O.0169 

« 

Css 

■.0053a. 0420 

•.0281a.0183 

■ .03340.0670 

■.0193a. 0290 

■.0163O.0327 

-.10870.1538 

a  a 

-.02310.0129 

Geis 

.0291a.02S0 

- .0036a.0091 

■.1247a.0892 

•.0015a. 0105 

.02040.0131 

.00280.0415 

-.31850.8692 

. 0031c. oosa 

a 

Ari 

•.0201a.0191 

.0ai1a.0084 

■.0072a.OS24 

.0027a.0065 

.0055O.0086 

.04600.0354 

a  a  a 

.0024O.0042 

a 

Per 

.0053S.01M 

■.0110a. 01 13 

.00920.0246 

■  .01340.0108 

■ .01550.0115 

■ .02290.0306 

-.03420.0458 

-.01050.0057 

7 

Ori 

■.00S3a.0092 

.0242a.0140 

.0069a.0091 

.0400a.0102 

.0090O.0081 

■.0040O.0143 

.02030.0194 

.0108O.0040 

0 

Tau 

.0106a.009e 

.0097a.0108 

.01960.0092 

.01310.0088 

.00550.0081 

.01930.0133 

.01950.0171 

.01250.0038 

0 

Aur 

.0009a.0106 

.0071a.010S 

.00670.0111 

.0044S.0102 

■ .00180.0115 

.0033O.0160 

.02870.0174 

.00520.0044 

7 

Gen 

.0128a. 0106 

.03S4a.01S6 

.0069O.0130 

.0002a.0120 

.0185O.0090 

.01470.0135 

.OOSOo.0130 

.01260.0044 

tt 

Gen 

■.0854a.0U2 

■.0834a.0204 

■ .09410.0145 

■.0624S.0128 

■.0833O.0092 

■ .07930.0161 

-.11800.0213 

-.08310.0052 

a 

CM! 

.0176a.01ia 

.OOSIs.0185 

■.0129S.0124 

.02850.0158 

.01430.0216 

.02450.0174 

■.0001O.0116 

.0084O.0054 

0 

Gael 

.0310a. 0128 

■.00S3a.0179 

.00370.0135 

.00330.0125 

.0060o  0101 

.0050O.0147 

.01690.0144 

.0093O.0049 

9 

Leo 

.0233a.0342 

.0608a.0624 

■ .00950.0237 

■ .06150.0382 

.007SO.0257 

.0904O.0840 

.00340.0321 

-.00030.0128 

0 

UMa 

•.0021a.0193 

... 

■ .01420.0241 

■.01640.0284 

.05610.0250 

.01210.0276 

.02660.0256 

.0095O.0099 

« 

UMa 

... 

... 

... 

.06740.0474 

a  •  a 

sea 

.06740.0474 

7 

Drs 

■.23280.3883 

woo 

... 

•  •  e 

a  a  a 

a  a  a 

-.23280.3883 

« 

Lyr 

.03120.2099 

ooo 

.02740.4014 

... 

aea 

a  a  a 

.03030.1860 

A 

Aqt 

.0366S.0774 

eoo 

■.13120.1512 

... 

e  e  a 

a  a  a 

-.00230.0688 

7 

Cy« 

■.0210O.0361 

•  •• 

■.1703a.0657 

.17740.1799 

•  a. 

a  a  a 

-.04860.0311 

« 

CVG 

■.04180.0397 

o.o 

■.0780a.056S 

.36980.1072 

a  a  a 

a  a 

-.01820.0310 

< 

Per 

.00270.0198 

*  0  o 

■ .08440.0978 

.07270.1214 

e  a  a 

sea 

.0010O.0191 
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TaM-E  V.  Derived  corrections  to  decliiutions  for  observed  stars  in  seconds  of  arc. 


•tar 

5  MOV 

6  Mv 

9  dov 

10  Hav 

11  Nev 

12  MOV 

14  Nov 

■  And 

-O.AdSaO.ASa 

o.oasao.osi 

a. a 

•0.069.0. 09S 

-0.096*0. 09S 

aaa 

a  . 

0.017*0.040 

0  Caa 

•0.U«a0.3S7 

e.U1a0.00A 

aaa 

0.026a0.M1 

0.167*0.171 

1.666*2.323 

a  a 

0.089*0.064 

a  Ca« 

-O.OOAaO.IBI 

•O.IMaO.OTI 

•0.0i7a0.291 

0.09S.0.106 

-0.139.0.133 

■0.394*0.647 

0.339*3.661 

•0.074*0.050 

1  tm 

O.ITSaO.IIS 

O.OaOaO.067 

•O.OStaO.A66 

O.niaO.OSI 

0.077*0.062 

0.106*0.193 

-1.693*3.301 

0.048*0.028 

.’an 

•o.aaAao.aw 

O.OSTaO.OIt 

O.OZSaO.260 

•O.OOIaO.OSA 

0.040*0.047 

0.091*0.162 

•  a. 

0.022x0.021 

m  P9r 

-e.noao.orr 

•O.iaiaO.OSS 

•O.OAOaO.IU 

-O.1SOaO.069 

-0.200*0.039 

•0.210*0.126 

-0.206*0.192 

•0.139*0.029 

t  Ori 

•O.KTaO.Ota 

•O.OIIaO.OOt 

•O.OTSaO.OAS 

0.061.0.044 

•0.022*0.037 

-0.064*0.063 

0.021*0.063 

•0.010*0.018 

$  Tau 

O.IMaO.OU 

O.ISIaO.OSr 

0.210a0.0t2 

0.106a0.04S 

0.129*0.042 

0.136*0.061 

0.200*0.073 

0.137*0.018 

0  hut 

O.OOSaO.nS 

•o.oasao.on 

•0.076a0.047 

-0.016a0.071 

0.040*0.032 

-0.013*0.071 

0.009*0.067 

•0.016*0.022 

1  Sau 

O.eOSaO.OAS 

-O.OAOaO.OOO 

-O.OSOaO.OM 

•0.046c0.047 

-0.031*0.046 

■0.003*0.061 

-0.069*0.031 

■0.031*0.018 

■  Can 

■O.OZtaO.Otr 

aO.OIOaO.IIO 

•O.eSSaO.046 

•0.060.0.069 

-0.063*0.042 

0.063*0.072 

-0.131*0.061 

-0.031*0.021 

a  CH< 

•O.OUaO.OAS 

■o.ooato.oaA 

•0.039a0.040 

•0.0S6a0.072 

-0.041*0.006 

0.010*0.067 

•0.026*0.040 

-0.043*0.020 

0  Saw 

•O.OISaO.mi 

o.ooato.iM 

•O.0Q6aO.O46 

O.006a0.0TS 

0.009*0.039 

0.031*0.064 

•0.069*0.040 

■0.011*0.022 

a  Lao 

•0.04680.216 

■O.SO»a0.616 

0.U1a0.1S6 

0.440a0.2S3 

0.067*0.126 

•0.306*0.366 

0.037*0.136 

0.081*0.075 

0  IMa 

•O.U1aO.U6 

•  a  a 

O.TSOaO.STO 

0.649.1.342 

-0.413*0.409 

•0.667*1.293 

•0.393*0.302 

•0.199*0.195 

a  UNa 

.  • 

a  a  a 

... 

•0.693*0.739 

-0.695*0.739 

1  Ora 

a  • 

•1.266a1.9M 

.a. 

.  a  a 

-1.266*1.996 

a  Lyr 

.  . 

0.332a1.160 

0.176.2.396 

a  a  a 

0.302*1.044 

a  Adi 

«  • 

•O.QS1aO.U2 

-1.329*1. 0S1 

a  >  a 

-0.226*0.408 

1  Cvo 

a. a 

•0.097a0.103 

•0.944*0.3U 

0.061*0.969 

-0.237x0.158 

•  cyg 

a  a  a 

•0.229a0,163 

•0.327*0.202 

1,969*0.460 

-0.077*0.135 

f  >ar 

... 

•0.012a0.090 

•0.370*0.S37 

0.122*0.634 

•0.019*0.087 

Sion  of  the  measurements  and  not  the  absolute  accuracy. 
Lsrger  systematic  errors  may  still  be  present  in  the  data. 

Figure  S  shows  the  difference  between  the  observed  and 
the  sur  positions  at  the  epoch  of  observation  for  the  ten 

sun  with  extensive  hour-angle  coverage.  The  error  bars  are 
1  standard  deviation.  Declinations  for  two  stan  deviate  sig¬ 
nificantly  from  the  assumed  positions  (more  than  5  standard 
deviations).  The  most  significant  deviation  in  right  ascen¬ 
sion  is  3  standard  deviations.  These  deviations  are  larger 
than  expected,  given  our  noise  estimates.  The  is  1.8  per 
degree  of  freedom  for  the  right  ascensions  and  2.6  for  the 
declinations.  If  we  assume  that  the  rms  errors  in  the  FKS 


Fic.  3.  Diflimiioe  between  the  fli  poeition  end  Uk  FKS  posi¬ 
tion  for  the  observed  sun  with  the  corretpondini  ( lo)  er¬ 
rors.  The  uneeitaittiiet  in  right  ascension  are  rooghly  flve 
timet  larger  than  the  uncertainties  in  declination  because  u>e 
uaad  a  north-eoMh  baseline.  A  large  easi-west  component 
win  decrease  the  errors  in  right  ascension. 


catalog  at  epoch  1987  average  around  0K)0S  in  right  ascen¬ 
sion  and  0’040  in  declination,  and  add  these  errors  quadrati- 
cally  to  our  formal  errors,  then  the  is  1  per  degree  of 
freedom,  as  expected  for  Gaussian  noise.  Alternatively,  the 
FKS  catalog  may  be  more  accurate  than  this  if  there  are 
systematic  errors  in  our  data. 

If  there  is  flexure  in  the  siderosut  systems,  or  if  the  differ¬ 
ential  atmospheric  delay  was  not  properly  estimated,  there 
may  be  a  systematic  error  in  the  star  positions  as  a  function 
of  declination.  The  dau  are  plotted  in  Fig.  6.  There  is  no 
apparent  systematic  trend  in  the  plot. 

IV  FUTURE  PLANS 

The  accuracy  of  our  measurements  of  sUr  positions  is  lim¬ 
ited  by  three  effects.  The  observed  delay  contains  an  addi¬ 
tional  term  proportional  to  the  index  of  refraction  of  the 


DECLINATION  (dagrats) 

Fic.  6.  Erron  in  decliiution  venus  tkclmuion.  There  are  no  tys- 
teinatic  efliecu  premnl  in  this  nnall  dau  let. 
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atmosphere.  Turbulence  produces  changes  in  the  index  of 
refractioo,  which  in  turn  causes  the  atmospheric  delay  to 
wander  a^t  the  expected  value.  The  expected  value  has 
been  estimated  and  removed,  but  the  fluctuations  limit  our 
ability  to  estimate  the  observed  delay.  This  noise  is 
which  implies  that  the  signal>to>noise  increases  as  the  sixth 
root  of  the  integration  time.  Second,  errors  in  the  estimate  of 
the  baseline  produce  a  residual  delay  which  is  then  interpret* 
ed  as  an  enw  in  the  calculated  star  positions.  The  baseline 
errors  are  dominated  by  rapid  drifts,  due  to  elTeca  such  as 
imperfect  ball  bearings  in  the  skferostats,  flexure,  and  the 
pivot*point  correction,  as  well  as  slower,  thermal  drifts.  Fi¬ 
nally,  correlations  in  the  fit  parameters  can  cause  the  wrong 
baseline  to  be  fit  to  the  data,  with  resulting  errors  in  the 
derived  star  positions. 

Four  major  improvements  will  be  incorporated  into  this 
system  during  1987.  Each  of  them  should  ^ve  a  significant 
rfect  on  the  sources  of  noise  discussed  above. 

First,  observations  will  be  performed  in  two  colors  simul¬ 
taneously.  After  the  light  from  the  two  siderostats  is  com¬ 
bined,  it  passes  through  a  dichroic  mirror.  This  directs  the 
blue  Ught  and  red  light  into  separate  photomultiplier  tubes. 
Because  the  atmosphere  is  a  dispersive  medium,  the  atmo¬ 
spheric  delay  is  different  in  red  and  blue  light,  with  the  differ¬ 
ence  proportional  to  the  atmospheric  delay  in  either  color 
alone.  Thus,  by  measuring  this  difference,  and  either  mea¬ 
suring  or  assuming  the  dispersion,  the  effects  of  the  atmo¬ 
sphere  can  be  greatly  reduced.  Some  preliminary  two-color 
dau  have  already  been  taken.  Colaviuer  nf.  (1987)  estimate 
that  two-color  observations  should  produce,  roughly,  a  fac¬ 
tor  of  S  decrease  in  the  atmospheric  noise. 

Second,  a  more  extensive  laMr  metrology  system  will  mea¬ 
sure  the  potion  of  each  siderostat  mirror  relative  to  the 
concrete  piers  on  which  they  are  mounted.  Changes  in  the 
baseline  due  to  motion  of  the  siderostats,  imperfect  ball  bear¬ 
ings,  and  thermal  motions  of  the  siderostats  will  be  measured 
and  corrected.  However,  drifts  due  to  changing  separation  of 
the  piers  will  still  have  to  be  estimated  from  the  data.  In 
addition,  an  internal  light  source  and  retroreflectors  mount¬ 
ed  on  the  siderostat  mirrors  will  provide  a  direct  measure¬ 
ment  of  the  delay  constant  up  to  20  times  each  night.  These 
measurements  should  reduce  the  systematic  variations  seen 
in  the  residual  delays.  However,  the  metrology  system  im¬ 
poses  limitt  on  siderostat  motion  which  will  restrict  sky  cov¬ 
erage. 

Third,  an  additional  baseline  with  a  substantial  E-W  com- 
iwnent  will  soon  be  in  operation.  Switching  between  base¬ 
lines  will  take  only  a  few  minutes  and  svill  be  possible  rough¬ 
ly  once  ^h  hour  during  observing.  This  should  reduce  the 
correlations  between  the  parameters,  allowing  the  determin¬ 
ation  of  absolute  declinations. 

Fii^ly,  a  more  ^histicated  star  tracker  svill  be  installed, 
allowing  obMrvations  of  objects  2-4  mag  fainter  than  our 
current  third-magnitude  limit.  This  will  ^low  more  stars  to 
be  observed  and  hmce  a  better  determination  of  the  baseline. 

Shaoera/.  (1987)  used  the  Mark  II  optical  interferometer 
with  a  3.1  m  baseline  and  achieved  a  precision  of  3  arcsec. 
They  estimated  that  the  precision  could  be  improved  by  a 
factor  of  300  by  improving  the  mechanical  and  thermal  sta¬ 
bility  of  the  system  and  1^  monitoring  the  siderostat  posi¬ 
tions  with  a  laser  metrology  system.  In  this  paper,  we  show  a 
factm*  of  1 00  improvement  over  the  Mark  II  system  and  ha ve 
not  yet  incorporated  the  laser  metrology  system.  Since  the 
noise  in  our  daU  is  still  a  factor  of  2  larger  than  that  expected 


from  the  atmosphere,  measurement  of  the  siderostat  posi¬ 
tions  should  improve  the  data  by  a  similar  factor. 

V.  SUMMARY 

Preliminary  astrometric  data  from  a  white-light  fringe¬ 
tracking  optical  interferometer  have  been  presented.  After 
fitting  a  baseline  to  the  data,  the  residual  delays  were  6.S/um 
rms.  The  noise  was  still  two  to  three  times  that  expected  from 
the  atmosphere  and  was  probably  dominated  by  mechanical 
roughness  and  thermal  drifts.  The  residual  delays  were  used 
to  determine  new  positions  for  the  observed  stars.  The  dau 
demonstrate  that  relative  astrometry  can  be  performed  with 
a  precision  of  Letter  than  22  mas  in  declination  and  S  ms  of 
time  in  right  ascension.  The  agreemoit  with  the  FK5  posi¬ 
tions  is  reasonable  considering  the  accuracy  of  the  FKS  cau- 
log  and  the  limiutions  of  our  small  daU  set. 

APPENDIX 

a}  Correction  for  the  A  tmosphere 

Since  the  index  of  refraction  of  air  is  greater  than  1,  the 
optical  path  length  differs  from  the  geometric  path  length. 
We  determined  this  difference  (1.89  m)  for  a  star  at  the 
zenith  by  integrating  the  index  of  refraction  as  a  function  of 
altitude.  We  used  a  standard  model  atmosphere  sUrting  at 
the  altitude  of  Mount  Wilson  and  the  standard  dependence 
of  the  index  of  refraction  on  temperature  and  density.  For  a 
plane-parallel  atmosphere,  possible  on  a  flat  Earth  without 
weather,  this  effect  is  identical  for  each  siderostat,  and  the 
atmosphere  has  no  systematic  effect  on  the  observed  delay. 
However,  pressure  gradients  along  the  baseline  will  result  in 
index-of-refraction  gradients  and  the  correction  to  the  delay 
will  be  different  for  the  two  siderosuts.  This  effect  is  small 
enough  to  ignore.  Making  the  sundard  assumption  that  the 
pressure  gradient  balances  the  Coriolis  force,  a  10  m/s  wind 
perpendicular  to  the  baseline  translates  into  a  change  of  opti¬ 
cal  path  length  oro.2S/im.  A  more  important  effect  is  due  to 
the  curvature  of  the  Earth.  The  zenith  direction  differs  for 
the  two  siderostats.  This  angle  is  0^4  for  a  1 2  m  baseline.  The 
difference  in  zenith  angle  to  the  star  can  be,  at  most,  this 
much.  We  assume  a  plane-parallel  atmosphere  at  each  sider¬ 
ostat  and  write  the  correction  to  the  delay  as 

Ad  =  1.89x10*  (sec(zl) -sec(z2) I,  (AI) 

where  zl  and  z2  are  the  zenith  angle  of  the  sur  as  viewed 
from  the  two  siderosuts  and  Ad  is  in  microns.  For  a  north- 
south  baseline,  this  effect  is  mainly  a  function  of  declination. 
However,  if  we  change  the  delay  constant  and  the  Z  compo¬ 
nent  of  the  baseline,  residual  delays  are  produced  which 
mimic  Eq.  (AI ).  Figure  7  shows  the  residual  delays  due  to 
the  atmosphere  (dashed  lines)  as  well  as  the  residuals  for 
Ai)o  =  3.0 /sm  and  Aff,  «  —  S.S3/sm  (solid  lines).  The  fit 
is  much  better  than  1  ^m  for  all  parts  of  the  sky. 

b)  Pivot-Point  Correction 

In  a  non-ideal  siderosut,  the  two  axes  will  not  intersect 
and  neither  will  be  in  the  surface  of  the  mirror.  Ibus  the 
mirror  will  translate  as  well  as  route  when  moving  between 
surs  in  different  parts  of  the  sky.  The  change  in  baseline  due 
to  this  motion  can  be  written  as 

SB  =  E  +  A  cos(^g),  (A2) 

where  E  is  the  disunce  of  the  center  of  the  mirror  from  the 
elevation  axis,  A  is  the  separation  of  the  two  axes  and  it 
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PIVOT  POINT  CORRECTION 


Fk.  7.  The  duhcd  bw  it  the  •tmotphcric  ooneetion,  in 
microat.  at  a  Ainelion  of  declination.  For  a  north-tooth 
baicline.  the  correction  it  afanott  independent  of  hour  an¬ 
gle.  The  tolid  hne  it  the  chanfe  in  delay  for  ACo  —  3.0^in 
and  AR,  «  —  S.S3  ^in.  Note  that  if  the  atmotpheric  cor¬ 
rection  nwre  not  included,  the  rctidual  delay  due  to  the 
atmoapherc  would  be  6t  to  within  a  micron  by  a  change  in 
the  bnwHne.  The  ttar  pMttiont  determined  from  the  dau 
would  not  be  affected. 


Fic.  8.  The  dathed  lines  show  the  liae  of  the  pivot-point 
correction  at  a  function  of  hour  angle  for  four  selected 
declinations.  The  solid  lines  show  the  change  in  delay 
for  the  same  declinationt  and  a  change  in  the  baseline  of 
Aff,  s  —  3.4/im  and  ACg— S.7^m.  Again,  the  errors 
in  the  residual  delay  caused  by  not  including  the  pivot- 
point  correction  are  lets  than  a  micron  for  all  parts  of 
the  sky  we  observed. 


the  routioaoftheelcvttion  axis.  The  Azimuth  (fixed)  axisof 
each  siderotut  is  pointed  toward  the  west  and  IS*  above  the 
horizon.  Since  the  orienutions  of  the  siderostats  are  almost 
identical,  we  can  write  the  change  in  delay  as 

Sdm(SB,~SEi}  (M-S-f  m  F) 

S.  (AS) 

where  the  subscripts  refer  to  the  two  siderostats  and  M,  S, 
and  F  are  unit  vectors  normal  to  the  siderostats  and  toward 
the  star  and  feed  mirror,  respectively.  For  identical  sidero- 
stats,  the  correction  vanishes.  After  installing  plastic  shims, 
A  and  £  were  measured  for  the  two  siderostats.  They  were  A  ^ 
»  —  9I.I/im,i4N  =  — 92.3/rm,£sB  — 3.4^m,and£M 


=  -4- 1.3  //m,  where  the  subscripts  refer  to  the  south  and 
north  siderostats.  For  these  values,  the  correction  is  shown 
as  the  dashed  lines  in  Fig.  8.  The  solid  lines  show  the  change 
in  delay  corresponding  to  a  change  in  the  baseline  (AB,. 
=  —  3.4pm  and  ADg  =  5.7f/m).  Except  for  more  than  4  hr 
to  the  east  at  low  declinations,  the  pivot-point  correction 
mimics  the  residual  delay  from  the  baMline  change.  If  either 
atmospheric  or  pivot-point  corrections  are  not  applied,  the 
baseline  parameters  will  not  correspond  to  the  geometric 
baseline,  but  the  star  positions  will  not  be  affected.  When  the 
pivot-point  corrections  were  applied  to  the  data,  the  residu¬ 
als  did  not  decrease.  Since  we  are  not  sure  of  the  stability  of 
the  offsets  at  the  micron  level,  we  decided  not  to  include  the 
corrections  in  the  final  fit. 
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